Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,963 A→T 11.4% intergenic (‑64/+97) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,9630AT11.4% 567.6 / 40.6 227intergenic (‑64/+97)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (97/87);  new base T (6/21);  total (106/121)
Fisher's exact test for biased strand distribution p-value = 3.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

GTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGG  >  NC_002947/5657816‑5658111
                                                                                                                                                   |                                                                                                                                                    
gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                                                    <  2:1920560/150‑1 (MQ=255)
gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                                                    <  1:5079928/150‑1 (MQ=255)
 tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                                   <  2:2754063/150‑1 (MQ=255)
 tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                                   <  1:2257441/150‑1 (MQ=255)
 tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                                   <  1:3033273/150‑1 (MQ=255)
   tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCa                                                                                                                                                 >  2:5345455/1‑150 (MQ=255)
    gCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCCAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACACGTACAGTGCAg                                                                                                                                                >  1:1421164/1‑150 (MQ=255)
         ggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                                                           <  1:4209244/150‑1 (MQ=255)
            cGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgtg                                                                                                                                        <  2:759970/150‑1 (MQ=255)
            cGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgtg                                                                                                                                        >  2:1271348/1‑150 (MQ=255)
             ggTTGCGCCGGCTGCAGCTGTGAGCATCCTTCCAGTGTCGCCTGGGCCGTGAGGGGCTATAGGGCTTTCATTCTGTTTCCTCGATTTGGGGCTGGTGTGGAGCCCTTCTCGGGCTCGCCCGGTGCCAGCGGTTCTGTGGAGTGTATGTgg                                                                                                                                       <  2:674861/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:3817300/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:4446335/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       >  1:3112690/1‑150 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:5291191/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:862716/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:982049/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  2:2916814/150‑1 (MQ=255)
             ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                       <  1:1489690/150‑1 (MQ=255)
               ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                     <  2:3903982/150‑1 (MQ=255)
               ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                     <  2:1807282/150‑1 (MQ=255)
               ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                     <  2:2135563/150‑1 (MQ=255)
               gtGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                     <  2:3538740/149‑1 (MQ=255)
                 gcgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGc                                                                                                                                   <  2:3914358/150‑1 (MQ=21)
                 gcgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGc                                                                                                                                   <  2:2443325/150‑1 (MQ=25)
                   gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCgg                                                                                                                                 <  2:2301753/150‑1 (MQ=255)
                   gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCGCCCACAGGTACAGTGCAGTCTCTGTGGGAGCgg                                                                                                                                 <  2:1419089/150‑1 (MQ=255)
                   gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCCGCTGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCgg                                                                                                                                 <  2:348833/150‑1 (MQ=21)
                    ccGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAGGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCggg                                                                                                                                <  1:5234755/150‑1 (MQ=255)
                     cGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGc                                                                                                                               >  2:4925513/1‑150 (MQ=255)
                           cagctgTGAGCATCCTGCCAGTGTCGCCTGGGCCGTCACGGCAAAACGCGCTTTCATTCAGTTTCCTCGTTTTGGGGCTGCTGTGCAGCCCTTCTGGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGccc                                                                                                                         <  1:1729133/150‑1 (MQ=38)
                           cagcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGccc                                                                                                                         >  2:642772/1‑150 (MQ=255)
                              cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCcgcg                                                                                                                      >  1:1254174/1‑150 (MQ=255)
                                gTGAGCATCCTGCACGTGTCGCCAGGGCCGTCACGGCAAAACGCGCTTTCATTCTGTTGCCTCGTTTTGGGGGTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCCCAGTTACAGTGCAGACTCTGTGGGAGCGGGCGAGCCCGCGaa                                                                                                                    <  2:4764806/150‑1 (MQ=37)
                                gTGAGCACCCCGCCCGTGTCGCCCGGGCCGTCACGGCCAAAAGCGCTTCCAGACCGTTACCGCGATTTGGGGCCGCTGTGCGGCCCTCCGCGGGCTCGCCCGCACCCACAGCGACTGCCCGGGCCCGGTGGGAGCGGGCGAGCCCGAGaa                                                                                                                    <  1:5487416/150‑1 (MQ=9)
                                           gCCTGTGGCGCCAGGGCCGCCACGGCCAAAAGCGCTTTCATCCAGTTCCCTCGATTGGGGGGTGCTGTGCCGCCCTTCTCGGGCTCGCCCGCGCCCACCGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                         <  2:4272950/150‑1 (MQ=37)
                                           gCCAGTGTCGCCTGGGCCGTCACGGCAAAACGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTGGGGCCCGCCCGCCCCCCCCGGAACAGGGCAGTCTGTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                         <  2:751850/150‑1 (MQ=38)
                                           gCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                         >  1:3674257/1‑150 (MQ=255)
                                           gCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTGCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                         <  2:3650412/150‑1 (MQ=255)
                                            cctgggtCGCCCGGGCCGTCCCGGCCAACCGCGCTTGCCTTCCGTGGCCGCGCTTGGGGGCTGCGGGGCAGCCCTTCTCGGGCTCGCCGGCGCCCCCCGGTCCCGGCCCGTCGCGGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGc                                                                                                        <  2:4281460/145‑1 (MQ=255)
                                                ggtCGCCCGGGCCGGCACGGCAAAACGCGCTTGCTGTCAGTTCCCGCGATTTGGGGGTGGGTTGCCGCCCTACTGGGGCTCGCCCGCGCCCACCGGTACTGGGCAGTGCCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGcc                                                                                                    <  2:1881108/149‑1 (MQ=255)
                                                 gtCGCCTGGGCCGTCACGGCAAAAAGCGGTTTCATTCCGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTGGGGCTCGCCGGCTCCCACCGGTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGccc                                                                                                   <  2:3807174/150‑1 (MQ=255)
                                                 gtCGCCTGGGCCGGCACGGCCAAAAGCGCTTACATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTCCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGccc                                                                                                   <  2:2987209/150‑1 (MQ=255)
                                                  tCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGcccc                                                                                                  >  2:716841/1‑150 (MQ=255)
                                                    gCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCctg                                                                                                >  1:1994559/1‑149 (MQ=255)
                                                     ccAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                                               >  1:1657029/1‑150 (MQ=255)
                                                      cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttct                                                                                              >  2:5063603/1‑150 (MQ=255)
                                                      cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTCTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttct                                                                                              >  1:2834895/1‑150 (MQ=255)
                                                       aGGGCCGTCACGGCCAAAAGCGCTTTCATCCCGTTCCCTCGCTTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttctt                                                                                             <  1:4685867/150‑1 (MQ=37)
                                                        gggCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCttcttc                                                                                            >  2:2823391/1‑150 (MQ=255)
                                                         ggCCGTCCGGCCCAAACGCGCTTTCATACAGTTGCCGCTATTTGGGGCTGCGGTGCAGCCCTTCGCGGGCTCGCCCGCGCCCCCCGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCAGccccaaatcga                                                                                           <  2:4230289/150‑8 (MQ=18)
                                                            cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                                                        <  2:3674257/150‑1 (MQ=255)
                                                             gTCACGGCCACAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAc                                                                                       >  2:1190118/1‑150 (MQ=255)
                                                             gTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAc                                                                                       >  2:2523496/1‑150 (MQ=255)
                                                                aCGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCGCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACcac                                                                                    >  1:784957/1‑150 (MQ=255)
                                                                    ccAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAt                                                                                >  1:4472650/1‑150 (MQ=255)
                                                                      aaaaGCGCTTTCAGGCAGTGTCCGCGAGTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGGACCGTGCCGGCGCGGGGGGAGCGGGCGGGCCCGCGCAGGGCCGCGCAGCGGCCCCTGGTTCCAGACcagcggggg                                                                              >  2:5184096/1‑143 (MQ=19)
                                                                            gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTc                                                                        >  1:4127972/1‑150 (MQ=255)
                                                                            gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTACCACAGGTACACTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAACGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAAACAACGtta                                                                        >  1:636881/1‑149 (MQ=39)
                                                                              tttCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCCCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGCAGGGCCGCACAGCGGCCCCGTCATAAAGACCACAAAACAACCTTCAg                                                                      >  1:1382739/1‑150 (MQ=37)
                                                                                  attcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCGCCCACAGGTACAGTGCAGTCGCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAACAACGTTAAGGTTg                                                                  >  1:2606242/1‑150 (MQ=255)
                                                                                   ttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                                                 <  2:1842081/150‑1 (MQ=255)
                                                                                                tttgggggCTGCTGTGCAGCCCTTCGCGGGCGCGCCCCCACCCCCCGGAACTGTGCAGTCTCTGGGGGTGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATc                                                    <  1:3089997/146‑1 (MQ=255)
                                                                                                   tGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                 <  1:1853046/150‑1 (MQ=255)
                                                                                                     gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAggg                                               >  2:1231188/1‑150 (MQ=255)
                                                                                                       gctgctGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACCAATCAACGTTCAGGTTGCCTCTTACCCGATCAAGGGCa                                             >  2:4550838/1‑150 (MQ=255)
                                                                                                        ctgctgTTCTGCCCTTCTCGGGCTGGCCCGCTCCCACCTGTACAGTGCTTTCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCaa                                            <  2:1421164/150‑1 (MQ=38)
                                                                                                         tgctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAAc                                           >  2:3706523/1‑150 (MQ=255)
                                                                                                         tgctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGGTCAGGTTGCCTATTACCCGATCAAGGGCAAc                                           >  1:4577929/1‑150 (MQ=255)
                                                                                                              tgCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgc                                      <  1:5054726/150‑1 (MQ=255)
                                                                                                                cggCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACAGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgcgc                                    >  2:370633/3‑150 (MQ=17)
                                                                                                                   cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                                 >  2:5585269/1‑150 (MQ=255)
                                                                                                                   cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAACCGATCAAGGGCAACCCGGCGCCtt                                 >  1:1252341/1‑150 (MQ=255)
                                                                                                                      ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgcg                              >  2:4149510/1‑150 (MQ=255)
                                                                                                                         tcGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtg                           <  2:2128500/150‑1 (MQ=255)
                                                                                                                         gggggCGCGCCCCCCCCCACTGGTACCGTGCTTTCTTTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtg                           <  1:3289518/148‑1 (MQ=255)
                                                                                                                          cGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtgg                          <  2:1677470/150‑1 (MQ=255)
                                                                                                                            ggCTCGCCCGCCCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                        <  2:1557123/150‑1 (MQ=255)
                                                                                                                            ggCCCGCCCGCCCCCCCCGGAACCGGGCCGCCGCTGTGGGCGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                        <  1:3450242/150‑1 (MQ=255)
                                                                                                                                cgcccgccCCCACCGGTACAGTGCCGTCCCGGGGGGCGCGGGCGAGCCCGCGCAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACACATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCt                    <  1:3370809/150‑1 (MQ=255)
                                                                                                                                  cccGCTCCCACTGGTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                  <  1:5004167/150‑1 (MQ=255)
                                                                                                                                  cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                  >  2:680254/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCCCCCGGTACAGTGCCGTCTCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                 <  2:2043236/150‑1 (MQ=39)
                                                                                                                                   ccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                 >  1:4047850/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                 >  2:3972440/1‑150 (MQ=25)
                                                                                                                                   ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                 >  1:2364974/1‑150 (MQ=21)
                                                                                                                                   ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                 >  1:2958855/1‑150 (MQ=21)
                                                                                                                                     gcccccACAGTTACCGTGCATTCTTTGTGGGTGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcgtc               <  1:4242418/147‑1 (MQ=255)
                                                                                                                                       tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTcgtcgt             >  1:1089617/1‑150 (MQ=255)
                                                                                                                                         ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc           >  1:4563798/1‑150 (MQ=255)
                                                                                                                                         ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc           >  1:133640/1‑150 (MQ=255)
                                                                                                                                         ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc           >  1:4830515/1‑150 (MQ=255)
                                                                                                                                          cacaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACt          >  1:4502341/1‑150 (MQ=255)
                                                                                                                                           acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt         >  2:3251230/1‑150 (MQ=255)
                                                                                                                                           acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt         >  2:245422/1‑150 (MQ=255)
                                                                                                                                           acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt         >  1:3279287/1‑150 (MQ=255)
                                                                                                                                           acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt         >  2:4652920/1‑150 (MQ=255)
                                                                                                                                             aGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtttt       >  2:1037258/1‑150 (MQ=255)
                                                                                                                                               gTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTcc     >  1:4610867/1‑150 (MQ=255)
                                                                                                                                                tACTGTGCAGTCTGTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCa    <  1:3516814/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGG  >  NC_002947/5657816‑5658111

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