Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 780,563 G→T 5.7% intergenic (‑90/+32) PP_0670 ← / ← glyA‑II bile acid/Na+ symporter family transporter/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947780,5630GT 5.7% 318.5 / 3.8 102intergenic (‑90/+32)PP_0670/glyA‑IIbile acid/Na+ symporter family transporter/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (46/48);  new base T (3/3);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 102 total aligned reads displayed.

CAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGC  >  NC_002947/780417‑780705
                                                                                                                                                    |                                                                                                                                              
cAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGcc                                                                                                                                             >  2:1748658/1‑150 (MQ=255)
 aGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCg                                                                                                                                            <  1:2265219/150‑1 (MQ=255)
 aGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCg                                                                                                                                            >  2:536300/1‑150 (MQ=255)
   tACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTCTAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCggg                                                                                                                                          <  1:1483806/150‑1 (MQ=255)
     cggcggCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGct                                                                                                                                        <  2:2385554/150‑1 (MQ=255)
      ggcggcGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctc                                                                                                                                       <  1:350909/150‑1 (MQ=255)
       gcggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctct                                                                                                                                      >  1:1428526/1‑150 (MQ=255)
        cggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GGGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCtt                                                                                                                                     <  1:775056/150‑1 (MQ=255)
        cggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCtt                                                                                                                                     <  2:767112/150‑1 (MQ=255)
         ggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGATGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCttt                                                                                                                                    <  1:240184/150‑1 (MQ=255)
             aaGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACt                                                                                                                                >  1:1449251/1‑150 (MQ=255)
                tAGATGGGGATCTGCCGGGTTTCGAGTTGCTATCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                             <  2:787656/150‑1 (MQ=255)
                 aGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                            <  2:1428526/150‑1 (MQ=255)
                  gATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGc                                                                                                                                              <  1:405878/131‑1 (MQ=255)
                  gATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGc                                                                                                                                              >  2:405878/1‑131 (MQ=255)
                  gATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAgg                                                                                                                           >  1:791746/1‑150 (MQ=255)
                     ggggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCaa                                                                                                                        >  2:747383/1‑150 (MQ=255)
                      gggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAg                                                                                                                       >  2:1461842/1‑150 (MQ=255)
                         aTCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgcgc                                                                                                                    <  1:1258880/150‑1 (MQ=255)
                            tGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                                 >  2:2261300/1‑150 (MQ=255)
                            tGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                                 <  1:652026/150‑1 (MQ=255)
                            tGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                                 >  2:1823973/1‑150 (MQ=255)
                             gCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                             >  1:275217/1‑137 (MQ=255)
                             gCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                             <  2:275217/137‑1 (MQ=255)
                               cGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCa                                                                                                              <  1:2329830/150‑1 (MQ=255)
                                gggTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAg                                                                                                             >  1:193222/1‑150 (MQ=255)
                                gggTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAg                                                                                                             <  1:409612/150‑1 (MQ=255)
                                  gTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGtt                                                                                                           <  1:536300/150‑1 (MQ=255)
                                   tttCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg                                                                                                          >  1:2369115/1‑150 (MQ=255)
                                   tttCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg                                                                                                          <  2:1619655/150‑1 (MQ=255)
                                    ttCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGc                                                                                                         >  2:654682/1‑150 (MQ=255)
                                     tCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGCGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGcc                                                                                                        <  1:1260673/150‑1 (MQ=255)
                                     tCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGcc                                                                                                        <  1:654682/150‑1 (MQ=255)
                                         gTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                             >  2:1696159/1‑125 (MQ=255)
                                         gTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                             <  1:1696159/125‑1 (MQ=255)
                                            gCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAAcc                                                                                                 >  1:557982/1‑150 (MQ=255)
                                              tctctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCgg                                                                                               <  1:2094931/150‑1 (MQ=255)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTcgtt                                                                                                                                <  1:1981318/115‑4 (MQ=38)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTcgtt                                                                                                                                >  2:1981318/1‑112 (MQ=38)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCgg                                                                                               >  1:2273847/1‑148 (MQ=255)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCgg                                                                                               <  2:2273847/148‑1 (MQ=255)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                             >  1:1555635/1‑150 (MQ=255)
                                                tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                             >  2:2176730/1‑150 (MQ=255)
                                                 ctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAg                                                                                            <  2:791746/150‑1 (MQ=255)
                                                    gTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGc                                                                                                                                         >  1:1089267/1‑102 (MQ=255)
                                                       aaGGTGCGGACGGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTCCTAGGGCAGCTATTCGGAGTCATT‑‑CCGGAAGCGCCCGGCACCGTCGCGGGCCCTTTCACTTCAGGCAAGCGTCAATCAGTTGCCGTAAACCGGCAGcttctt                                                                                      <  1:544204/150‑1 (MQ=255)
                                                         ggTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGc                                                                                    <  2:1057900/150‑1 (MQ=255)
                                                          gTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctc                                                                                                                                       >  1:504280/1‑98 (MQ=255)
                                                          gTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctc                                                                                                                                       <  2:504280/98‑1 (MQ=255)
                                                          gTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCa                                                                                   <  2:207979/150‑1 (MQ=255)
                                                               gCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc                                                                                                                      <  1:583553/110‑1 (MQ=255)
                                                               gCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc                                                                                                                      >  2:583553/1‑110 (MQ=255)
                                                                 cTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAAt                                                                                                                >  1:1318385/1‑114 (MQ=255)
                                                                 cTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAAt                                                                                                                <  2:1318385/114‑1 (MQ=255)
                                                                 cTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGcc                                                                            >  2:180365/1‑150 (MQ=255)
                                                                 cTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGcc                                                                            >  2:180273/1‑150 (MQ=255)
                                                                  tgtgCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGc                                                                                                         <  2:2350011/120‑1 (MQ=255)
                                                                  tgtgCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGc                                                                                                         >  1:2350011/1‑120 (MQ=255)
                                                                     gCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGa                                                                        <  1:2261300/150‑1 (MQ=255)
                                                                      cAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGAc                                                                       <  1:1461842/150‑1 (MQ=255)
                                                                           tGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtct                                                                  >  2:2086586/1‑150 (MQ=255)
                                                                            ggAGGTATTGGGGTCTCTGTCTGTGGACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCGTTTCGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtctc                                                                 <  1:747383/150‑1 (MQ=255)
                                                                              aGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCac                                                               >  2:129025/1‑150 (MQ=255)
                                                                              aGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCac                                                               >  2:139257/1‑150 (MQ=255)
                                                                               ggTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCggg                                                                                                                                          >  2:1639981/1‑74 (MQ=255)
                                                                               ggTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCggg                                                                                                                                          <  1:1639981/74‑1 (MQ=255)
                                                                               ggTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacg                                                              >  1:2303684/1‑150 (MQ=255)
                                                                                gTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgcgc                                                                                                                    >  2:247084/1‑95 (MQ=255)
                                                                                gTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgcgc                                                                                                                    <  1:247084/95‑1 (MQ=255)
                                                                                  aTTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggac                                                           <  1:2176730/150‑1 (MQ=255)
                                                                                   ttGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggacg                                                          >  2:1667420/1‑150 (MQ=255)
                                                                                   ttGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCAATTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggacg                                                          <  1:1748658/150‑1 (MQ=255)
                                                                                     gggggCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGc                                                        <  2:677553/150‑1 (MQ=255)
                                                                                       gggCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAAGCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGt                                                      <  1:425047/150‑1 (MQ=255)
                                                                                                 tttGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCac                                                               >  2:278219/1‑131 (MQ=255)
                                                                                                 tttGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCac                                                               <  1:278219/131‑1 (MQ=255)
                                                                                                        gAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCt                                           <  1:1444325/144‑1 (MQ=255)
                                                                                                        gAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCt                                           >  2:1444325/1‑144 (MQ=255)
                                                                                                        gAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgttg                                     >  2:1163573/1‑150 (MQ=255)
                                                                                                               gCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg                                >  1:2309590/1‑150 (MQ=255)
                                                                                                               gCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg                                >  1:1804948/1‑150 (MQ=255)
                                                                                                                cGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCatg                               <  2:2369115/150‑3 (MQ=255)
                                                                                                                  gTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTc                             >  1:204280/1‑150 (MQ=255)
                                                                                                                     cgcCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGc                          <  2:193222/150‑1 (MQ=255)
                                                                                                                        cTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAg                       >  2:2103193/1‑150 (MQ=255)
                                                                                                                         tGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAgg                      <  2:1449251/150‑1 (MQ=255)
                                                                                                                             ttAAACGAAAGAGCCCGGCACTGAGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGt                  >  2:1185594/1‑150 (MQ=255)
                                                                                                                               aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCg                <  2:1804948/150‑1 (MQ=255)
                                                                                                                                aaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGc               >  2:1128131/1‑150 (MQ=255)
                                                                                                                                 aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCa              >  1:308527/1‑150 (MQ=255)
                                                                                                                                   gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGa            <  2:2309590/150‑1 (MQ=255)
                                                                                                                                   gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGa            <  2:204280/150‑1 (MQ=255)
                                                                                                                                   gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCAATTTAGGAAAGCGAAGATCTCTTGCCGTAAACCGGCAGCATCGGGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGa            >  1:1863615/1‑150 (MQ=255)
                                                                                                                                    aaaGAGCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGAt           >  2:57404/1‑150 (MQ=255)
                                                                                                                                      agagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTc                                                     <  2:2372748/106‑1 (MQ=39)
                                                                                                                                      agagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTc                                                     >  1:2372748/1‑106 (MQ=39)
                                                                                                                                      agagCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATcc         >  2:2028679/1‑150 (MQ=255)
                                                                                                                                      agagCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATcc         >  2:2028678/1‑150 (MQ=255)
                                                                                                                                          cccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGcc     >  2:1884911/1‑150 (MQ=255)
                                                                                                                                             ggCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGc  <  1:1884911/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGG‑‑GCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGC  >  NC_002947/780417‑780705

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