Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,176,313 C→T 5.0% W174* (TGG→TAG)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,3130CT 5.0% 250.7 / 2.1 80W174* (TGG→TAG) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (2/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

GGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAG  >  NC_000913/3176165‑3176461
                                                                                                                                                    |                                                                                                                                                    
ggTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGcc                                                                                                                                                      <  6:145274/149‑1 (MQ=255)
ggTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGcc                                                                                                                                                      <  6:278589/149‑1 (MQ=255)
  tCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCaa                                                                                                                                                    <  6:372396/149‑1 (MQ=255)
    gACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAAc                                                                                                                                                   >  1:391469/1‑148 (MQ=255)
              ccAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGc                                                                                                                                         <  6:127345/148‑1 (MQ=255)
               cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAg                                                                                                                                       <  2:388237/149‑1 (MQ=255)
               cagcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAg                                                                                                                                       <  2:17568/149‑1 (MQ=255)
                agcagGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAgg                                                                                                                                      >  4:268726/1‑149 (MQ=255)
                    ggcggcGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGc                                                                                                                                  <  3:342277/149‑1 (MQ=255)
                     gcggcgACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCg                                                                                                                                 >  7:191271/1‑149 (MQ=255)
                     gcggcgACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCg                                                                                                                                 >  2:237500/1‑149 (MQ=255)
                            ccATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatgatg                                                                                                                          <  1:176630/149‑1 (MQ=255)
                             cATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatgatga                                                                                                                         >  6:100676/1‑149 (MQ=255)
                               ttCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGatgatgatg                                                                                                                       >  1:354447/1‑149 (MQ=255)
                                tCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCTAGCAAACAACAACCCGCATATAGCGGatgatgatg                                                                                                                       >  2:217259/1‑148 (MQ=255)
                                tCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGatgatgatg                                                                                                                       >  5:353568/1‑148 (MQ=255)
                                  cAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcag                                                                                                                    >  7:293263/1‑149 (MQ=255)
                                   aGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagc                                                                                                                   <  1:217259/149‑1 (MQ=255)
                                       gagaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagca                                                                                                               <  8:191271/149‑1 (MQ=255)
                                          aTCCAGCTCCTGATGTATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcagc                                                                                                             <  4:36016/148‑1 (MQ=255)
                                           tCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatgatg                                                                                                                          <  2:220173/134‑1 (MQ=255)
                                           tCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatgatg                                                                                                                          >  1:220173/1‑134 (MQ=255)
                                           tCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCaa                                                                                                           <  1:237500/149‑1 (MQ=255)
                                            ccAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAAc                                                                                                          >  1:239240/1‑149 (MQ=255)
                                               gCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTa                                                                                                       <  6:101100/149‑1 (MQ=255)
                                               gCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTa                                                                                                       <  2:238589/149‑1 (MQ=255)
                                                cTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTAt                                                                                                      <  2:354447/149‑1 (MQ=255)
                                                    tgatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCAGCGTTACGCAGACTGTGTAGATCGAGCCAACTACAACCCGCAGGAAGCTCATGATGATGCAGCAGCAGCAACGTATGGCg                                                                                                  >  6:54381/1‑149 (MQ=255)
                                                     gatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCTAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGt                                                                                                 >  7:236257/1‑149 (MQ=255)
                                                     gatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCTAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGt                                                                                                 >  7:230638/1‑149 (MQ=255)
                                                     gatgaATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGt                                                                                                 >  1:273573/1‑149 (MQ=255)
                                                         aaTATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCt                                                                                              <  7:309955/148‑1 (MQ=255)
                                                              gACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCa                                                                                        >  4:295491/1‑149 (MQ=255)
                                                                 cGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCTAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCg                                                                                     <  8:236257/149‑1 (MQ=255)
                                                                 cGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCTAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCg                                                                                     <  8:230638/149‑1 (MQ=255)
                                                                      aGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAg                                                                                >  3:97893/1‑149 (MQ=255)
                                                                       gCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGt                                                                               <  7:323006/149‑1 (MQ=255)
                                                                        cAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGtt                                                                              >  5:123183/1‑149 (MQ=255)
                                                                         aaGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGttt                                                                             <  3:283260/149‑1 (MQ=255)
                                                                           gTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCGGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTAc                                                                           <  5:54381/149‑1 (MQ=255)
                                                                           gTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTAc                                                                           <  2:391469/149‑1 (MQ=255)
                                                                            tATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACg                                                                          <  8:293263/149‑1 (MQ=255)
                                                                              tttGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGTCGCATCGGCCAGTTTAc                                                                           <  7:173744/146‑1 (MQ=255)
                                                                              tttGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTAc                                                                           >  8:173744/1‑146 (MQ=255)
                                                                                 gACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTc                                                                       >  7:318898/1‑147 (MQ=255)
                                                                                  aCGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCa                                                                     <  3:268726/148‑1 (MQ=255)
                                                                                      gCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCGACGTATGGCg                                                                                                  <  6:59282/115‑1 (MQ=255)
                                                                                      gCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCGACGTATGGCg                                                                                                  >  5:59282/1‑115 (MQ=255)
                                                                                      gCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCAc                                                                >  8:88248/1‑149 (MQ=255)
                                                                                       cGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCAc                                                                <  2:69225/148‑1 (MQ=255)
                                                                                           aaCTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAg                                                           <  2:273573/149‑1 (MQ=255)
                                                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGttt                                                                             <  2:255748/128‑1 (MQ=255)
                                                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGttt                                                                             >  1:255748/1‑128 (MQ=255)
                                                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                         >  7:86410/1‑148 (MQ=255)
                                                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                         >  1:100593/1‑148 (MQ=255)
                                                                                              ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                         >  5:164378/1‑148 (MQ=255)
                                                                                                cGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCAAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGa                                                       <  5:207467/148‑1 (MQ=255)
                                                                                                 gCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaaa                                                     <  6:309965/149‑1 (MQ=255)
                                                                                                 gCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaaa                                                     <  5:100676/149‑1 (MQ=255)
                                                                                                  ccAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAAc                                                    <  7:91978/149‑1 (MQ=255)
                                                                                                  ccAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAAc                                                    >  7:388516/1‑149 (MQ=255)
                                                                                                         ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                         >  4:274902/1‑137 (MQ=255)
                                                                                                         ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCt                                                         <  3:274902/137‑1 (MQ=255)
                                                                                                         ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctca                                             <  1:34644/149‑1 (MQ=255)
                                                                                                           tATCCAGTTCTACCGCGTTACGCAGACTGAGTTTATCGAGCCACCTACAACCCGCAGGTAGCGGAACATCAAGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCTGCAACTTTACGTTACATCCACTCAAGCTGAAACTCgctcagc                                           >  3:282241/1‑149 (MQ=255)
                                                                                                           tATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagc                                           >  3:356110/1‑149 (MQ=255)
                                                                                                             tCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagct                                          >  3:165699/1‑148 (MQ=255)
                                                                                                                       ccGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCgg                                >  5:306265/1‑148 (MQ=255)
                                                                                                                           gTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGc                           <  7:234951/149‑1 (MQ=255)
                                                                                                                              aCGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCa                         <  3:295491/148‑1 (MQ=255)
                                                                                                                               cGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCaaa                       <  8:217380/149‑1 (MQ=255)
                                                                                                                                   gACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACAcc                   >  5:280422/1‑149 (MQ=255)
                                                                                                                                   gACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCGCACCGTGCGGCACGCCAAACAcc                   >  1:235784/1‑149 (MQ=255)
                                                                                                                                    aCTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCt                  <  6:123183/149‑1 (MQ=255)
                                                                                                                                       gtgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                <  7:28174/148‑1 (MQ=255)
                                                                                                                                       gtgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                >  8:28174/1‑148 (MQ=255)
                                                                                                                                        tgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                <  8:318898/147‑1 (MQ=255)
                                                                                                                                        tgtTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGc               <  6:353568/148‑1 (MQ=255)
                                                                                                                                          ttGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTAt            >  4:169825/1‑149 (MQ=255)
                                                                                                                                            gATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATc           >  8:70404/1‑148 (MQ=255)
                                                                                                                                             aTCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATc           >  5:234771/1‑147 (MQ=255)
                                                                                                                                                  gCCAACTACAACCCGCAGGTAGCGGATGATGCTGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAAc    <  6:375760/149‑1 (MQ=255)
                                                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAg        <  7:181243/144‑1 (MQ=255)
                                                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAg        >  8:181243/1‑144 (MQ=255)
                                                                                                                                                   ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCaa     >  6:246745/1‑147 (MQ=255)
                                                                                                                                                    cAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAg  <  2:239240/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAG  >  NC_000913/3176165‑3176461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: