Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 418,298 C→G 5.3% L26L (CTC→CTG PP_0345 → heptose kinase WapP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947418,2980CG 5.3% 375.2 / 9.5 114L26L (CTC→CTGPP_0345heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base C (54/54);  new base G (2/4);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 115 total aligned reads displayed.

TGCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAAA  >  NC_002947/418151‑418447
                                                                                                                                                   |                                                                                                                                                     
tgctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCAGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCa                                                                                                                                                     >  1:12868/1‑150 (MQ=255)
  ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCCAGCCACGCTCTGGGGCTGCACTCTAATCGCGCGGACCACACAGGCACCGGGGTACTAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   >  2:1477915/1‑150 (MQ=255)
  ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   <  2:1401962/150‑1 (MQ=255)
  cacGCCTCGCTGTATCCCAAGTCCGAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   <  2:12868/148‑1 (MQ=255)
     gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCg                                                                                                                                                >  2:1851980/1‑150 (MQ=255)
       ctggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGaa                                                                                                                                              <  2:1002620/150‑1 (MQ=255)
          gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAcac                                                                                                                                           <  2:1755489/150‑3 (MQ=255)
          gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                           <  2:555998/150‑1 (MQ=255)
                gcgcAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacc                                                                                                                                     >  2:281308/1‑150 (MQ=255)
                gcgcAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacc                                                                                                                                     >  1:2237501/1‑150 (MQ=255)
                   cAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                                                  >  2:428155/1‑150 (MQ=255)
                    aaGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcg                                                                                                                                 >  2:388493/1‑150 (MQ=255)
                      ggCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgac                                                                                                                               <  2:782703/150‑1 (MQ=255)
                       gCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacc                                                                                                                              <  1:689692/150‑1 (MQ=255)
                         ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                            >  1:1683512/1‑150 (MQ=255)
                            gAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                           <  1:1688266/132‑1 (MQ=255)
                            gAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                           >  2:1688266/1‑132 (MQ=255)
                            caaaCTCTCCGAACGCAAGCAACGCTACAGGGGGGCACTCTGATGGCGAGTTGCACACTGGCGCCGGGGTACGAACACCTGGCGGCAGACTAAGGCAGCCTCGACGCGGTCTTCGCCCTCCAGCGCGAGCGGCTCACCCGCGACCCGCTc                                                                                                                         >  2:293232/2‑150 (MQ=255)
                                 tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGa                                                                                                                    >  1:1494120/1‑150 (MQ=255)
                                  ctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                   <  2:222625/150‑1 (MQ=255)
                                  ctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                   <  2:806916/150‑1 (MQ=255)
                                      gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGACCGAATGAGCCGCGACCCGCTCAGTGAAGTgg                                                                                                               >  1:612579/1‑150 (MQ=255)
                                        tCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAc                                                                                                                                             <  1:232232/118‑1 (MQ=255)
                                        tCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAc                                                                                                                                             >  2:232232/1‑118 (MQ=255)
                                          gcaagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACg                                                                                                           >  2:1187895/1‑150 (MQ=255)
                                          gcaagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACg                                                                                                           >  2:342623/1‑150 (MQ=255)
                                           gaatcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGc                                                                                                          >  1:84739/5‑150 (MQ=255)
                                             agcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCAt                                                                                                        <  1:388493/150‑1 (MQ=255)
                                                aaCGCTATGGGGATGCACTCAGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGa                                                                                                     <  2:1444750/150‑1 (MQ=255)
                                                    cTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGt                                                                                                 <  2:84739/150‑1 (MQ=255)
                                                        ggggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                             <  2:2237501/150‑1 (MQ=255)
                                                        ggggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                             >  2:1841320/1‑150 (MQ=255)
                                                         gggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACgg                                                                                            <  1:313292/150‑1 (MQ=255)
                                                           gATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGa                                                                                          <  2:1009192/150‑1 (MQ=255)
                                                            aTGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGAt                                                                                         <  1:342623/150‑1 (MQ=255)
                                                             tGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATc                                                                                        <  1:852923/150‑1 (MQ=255)
                                                                aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATc                                                                                                       >  1:2197899/1‑132 (MQ=255)
                                                                aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                     >  2:1407865/1‑150 (MQ=255)
                                                                aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGGGAGGTGGTACGCATCGAGCGTGACGGGAtagat                                                                                     >  2:193526/1‑146 (MQ=255)
                                                                aCTCTGATGGCGAGTTGGACACTGGCGCCGGGCTACGACCACCTCGCGGCGGACTTCCGCAGCCTCGACGCGGTGTGCGCCCTCCAGGGCGAACGGCTGACACCCGACCCGCGCCGTGTAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                     <  1:293232/150‑1 (MQ=255)
                                                                  tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGTTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   <  1:1851980/150‑1 (MQ=255)
                                                                  tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   >  2:2222565/1‑150 (MQ=255)
                                                                  tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   <  1:1673552/150‑1 (MQ=255)
                                                                  tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   >  1:8691/1‑150 (MQ=255)
                                                                    tGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAactact                                                                                 >  2:1129287/1‑150 (MQ=255)
                                                                       tGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   <  1:1657359/145‑1 (MQ=255)
                                                                       aggCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAacta                                                                                   >  2:1657359/2‑145 (MQ=255)
                                                                         gCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTc                                                                            <  2:1867473/150‑1 (MQ=255)
                                                                         gCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTc                                                                            >  2:535899/1‑150 (MQ=255)
                                                                           gAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCaa                                                                          >  2:28753/1‑150 (MQ=255)
                                                                                    gcTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCGTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                 >  2:1744553/2‑150 (MQ=255)
                                                                                    acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                 >  2:2006915/1‑150 (MQ=255)
                                                                                    acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                 >  2:203516/1‑150 (MQ=255)
                                                                                    acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                 <  1:1129287/150‑1 (MQ=255)
                                                                                      tGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                                              <  2:1695729/119‑1 (MQ=255)
                                                                                      tGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCg                                                               >  1:395553/1‑150 (MQ=255)
                                                                                       tgcgcTGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgc                                                                        >  2:1213339/2‑140 (MQ=255)
                                                                                       tgcgcTGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgc                                                                        <  1:1213339/139‑1 (MQ=255)
                                                                                            cGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCg                                                         <  1:2006915/150‑1 (MQ=255)
                                                                                            cGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCg                                                         <  1:28753/150‑1 (MQ=255)
                                                                                             gcggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCgg                                                        <  2:1465880/148‑1 (MQ=255)
                                                                                                gTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCaa                                                     >  2:2409371/1‑150 (MQ=255)
                                                                                                 tACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTa                                                                                                             >  2:962920/1‑93 (MQ=255)
                                                                                                 tACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTa                                                                                                             <  1:962920/93‑1 (MQ=255)
                                                                                                  aCGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGc                                                   <  2:1683512/150‑1 (MQ=255)
                                                                                                   cGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGca                                                  <  1:1813970/150‑1 (MQ=255)
                                                                                                     aaCACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATc                                                                                                       <  2:2399354/95‑1 (MQ=255)
                                                                                                          cTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   <  2:737258/46‑1 (MQ=255)
                                                                                                          cTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   >  1:737258/1‑46 (MQ=255)
                                                                                                           tGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                        >  2:1356369/1‑72 (MQ=255)
                                                                                                           tGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                        <  1:1356369/72‑1 (MQ=255)
                                                                                                           tGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGcc                                          <  1:1383194/150‑1 (MQ=255)
                                                                                                               ggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAg                                                    <  2:73178/136‑1 (MQ=255)
                                                                                                               ggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAg                                                    >  1:73178/1‑136 (MQ=255)
                                                                                                                gcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                                <  2:1832372/75‑1 (MQ=255)
                                                                                                                gcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                                >  1:1832372/1‑75 (MQ=255)
                                                                                                                gcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgggg                                                            <  2:2251377/127‑1 (MQ=255)
                                                                                                                gcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgggg                                                            >  1:2251377/1‑127 (MQ=255)
                                                                                                                   gACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCGTCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTg                                  >  2:249283/1‑150 (MQ=255)
                                                                                                                     cTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTggg                                <  2:612579/150‑1 (MQ=255)
                                                                                                                      ttgggCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgggg                               >  1:2402313/4‑150 (MQ=255)
                                                                                                                        cGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                     <  1:606385/94‑1 (MQ=255)
                                                                                                                        cGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                     >  2:606385/1‑94 (MQ=255)
                                                                                                                         ggCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAactac                                                                                  <  2:2287951/96‑1 (MQ=255)
                                                                                                                         ggCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAactac                                                                                  >  1:2287951/1‑96 (MQ=255)
                                                                                                                          gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                           <  2:1494120/150‑1 (MQ=255)
                                                                                                                              ccTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                       >  1:1720077/1‑150 (MQ=255)
                                                                                                                               cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGt                                                                                                                 >  2:1742909/1‑59 (MQ=255)
                                                                                                                               cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGt                                                                                                                 <  1:1742909/59‑1 (MQ=255)
                                                                                                                               cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGc                      >  1:50701/1‑150 (MQ=255)
                                                                                                                                    cgcgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACGCGCTCAGTGAAGTTGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCaggaa                 <  2:1170520/150‑5 (MQ=255)
                                                                                                                                       ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTc                                                                            >  1:113466/1‑88 (MQ=255)
                                                                                                                                       ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGc              <  1:203516/150‑1 (MQ=255)
                                                                                                                                        gtcctCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                                               <  1:151407/80‑1 (MQ=255)
                                                                                                                                          gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCa                         <  2:1240001/136‑1 (MQ=255)
                                                                                                                                          gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCa                         >  1:1240001/1‑136 (MQ=255)
                                                                                                                                           ttCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGaa          <  1:1477915/150‑1 (MQ=255)
                                                                                                                                             cGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATg        <  1:1841320/150‑1 (MQ=255)
                                                                                                                                               cccTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGGTCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGAGGCAGCCACGCATCAAGCCCGAATGGc      >  2:1213323/1‑150 (MQ=255)
                                                                                                                                                   ccAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCaaaa  >  2:374664/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAAA  >  NC_002947/418151‑418447

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