Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,481,523 G→A 5.0% A74A (GCC→GCT rhlB ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,481,5230GA 5.0% 627.6 / 12.8 197A74A (GCC→GCTrhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base G (99/88);  new base A (4/6);  total (103/94)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

TAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGT  >  NC_002947/1481374‑1481669
                                                                                                                                                     |                                                                                                                                                  
tAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTg                                                                                                                                                    >  1:4536731/1‑150 (MQ=255)
  aCCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                                                  <  1:452166/150‑1 (MQ=255)
   cccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGcc                                                                                                                                                 >  1:584045/1‑150 (MQ=255)
    ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCt                                                                                                                                                <  1:1913002/150‑1 (MQ=255)
      aGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCtgt                                                                                                                                              <  2:2466858/150‑1 (MQ=255)
        gACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTgg                                                                                                                                            <  2:3160663/150‑1 (MQ=255)
          ccTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggct                                                                                                                                          >  2:865763/1‑150 (MQ=255)
           cTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggctt                                                                                                                                         >  1:5258781/1‑150 (MQ=255)
             gcgcCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgg                                                                                                                                       >  1:93482/1‑150 (MQ=255)
                 cTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggctt                                                                                                                                   <  1:3079233/150‑1 (MQ=255)
                 cTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggctt                                                                                                                                   <  1:349942/150‑1 (MQ=255)
                 cTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggctt                                                                                                                                   <  2:2568686/150‑1 (MQ=255)
                 cTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggctt                                                                                                                                   <  1:2217410/150‑1 (MQ=255)
                   ggATCGGCGTGCAGGACGGGAAGCCCAGGTCGTGGATGGCGTGCCTCAGCTCGTTGGACAGCTTGAAGTCGTGGACGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGc                                                                                                                                 <  2:1276690/150‑1 (MQ=255)
                    gATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCg                                                                                                                                >  2:3663157/1‑150 (MQ=255)
                    gATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCg                                                                                                                                >  1:4287408/1‑150 (MQ=255)
                               aGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTcgcggcgcg                                                                                                                     >  2:226147/1‑150 (MQ=255)
                                gTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTcgcggcgcgg                                                                                                                    >  1:2696191/1‑150 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGt                                                                                                                   >  1:2169579/1‑150 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGt                                                                                                                   >  2:3758251/1‑150 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTGTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGt                                                                                                                   <  2:4501415/150‑1 (MQ=255)
                                    gggAAGCCCAGGTCGTGGATGGCGTGCATAAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCATGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCGTGGGCTTGCGTTCGCGGCGCGGTTTg                                                                                                                >  1:2813602/1‑150 (MQ=255)
                                      gAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTc                                                                                                              >  2:2977450/1‑150 (MQ=255)
                                         gCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCtt                                                                                                           <  2:1601261/150‑1 (MQ=255)
                                          cccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                                                          >  1:1293132/1‑150 (MQ=255)
                                            cAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGc                                                                                                        >  1:804393/1‑150 (MQ=255)
                                            cAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGc                                                                                                        >  1:4450505/1‑150 (MQ=255)
                                             aGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGcc                                                                                                       >  1:2049168/1‑150 (MQ=255)
                                                    ggATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGAt                                                                                                >  2:5437791/1‑150 (MQ=255)
                                                    ggATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGACTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGAt                                                                                                <  2:2949001/150‑1 (MQ=255)
                                                     gATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATc                                                                                               <  2:3102957/150‑1 (MQ=255)
                                                       gggCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGACGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCgg                                                                                             <  1:2443339/149‑1 (MQ=255)
                                                        ggCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGc                                                                                            >  1:4588334/1‑150 (MQ=255)
                                                        ggCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGc                                                                                            >  1:4588625/1‑150 (MQ=255)
                                                           gTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAgg                                                                                         >  2:3362106/1‑150 (MQ=255)
                                                                tCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCg                                                                                    >  2:1760357/1‑150 (MQ=255)
                                                                   gCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgg                                                                                                                                       <  2:2779012/96‑1 (MQ=255)
                                                                   gCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCa                                                                                           <  2:4853313/140‑1 (MQ=255)
                                                                   gCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCa                                                                                           >  1:4853313/1‑140 (MQ=255)
                                                                   gCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTcg                                                                                 >  1:1377854/1‑150 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                             <  1:5340302/150‑1 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                             <  2:4599647/150‑1 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                             <  2:93482/150‑1 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGGTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                             >  2:724872/1‑150 (MQ=255)
                                                                            aCAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAg                                                                        >  1:4214111/1‑150 (MQ=255)
                                                                                cTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTg                                                                    <  1:558307/150‑1 (MQ=255)
                                                                                 ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGc                                                                   >  2:391848/1‑150 (MQ=255)
                                                                                    aaGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCa                                                                <  1:4405263/150‑1 (MQ=255)
                                                                                       tCGTGGAAGCGGGTTTTGCCTTCCCGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGcc                                                             >  1:3728714/1‑150 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                           >  1:2866698/1‑150 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                           <  2:2870230/150‑1 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                           >  1:1198033/1‑150 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                           >  2:505185/1‑150 (MQ=255)
                                                                                           ggAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTgg                                                         <  1:4765329/150‑1 (MQ=255)
                                                                                                 gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                   <  2:122762/150‑1 (MQ=255)
                                                                                                 gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                   <  1:1767869/150‑1 (MQ=255)
                                                                                                 gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGGTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                   >  2:4059622/1‑150 (MQ=255)
                                                                                                  ggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCggg                                                  <  2:432857/150‑1 (MQ=255)
                                                                                                    ttttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGt                                                <  2:2366812/150‑1 (MQ=255)
                                                                                                         ccTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCggcggc                                                                          >  1:2596228/1‑119 (MQ=255)
                                                                                                          cTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAg                                                                                          >  2:4371560/1‑102 (MQ=255)
                                                                                                           ttCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATg                                         >  1:2293352/1‑150 (MQ=255)
                                                                                                            tCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg                                        <  1:3810354/150‑1 (MQ=255)
                                                                                                            tCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg                                        <  2:251331/150‑1 (MQ=255)
                                                                                                             ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                       <  2:1377883/150‑1 (MQ=255)
                                                                                                             ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCACAGCCTTGGGGGCGGGGTTGGCATggg                                       <  2:2478295/150‑1 (MQ=255)
                                                                                                                gCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTg                                    >  2:4269496/1‑150 (MQ=255)
                                                                                                                     tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGcc                               >  1:401773/1‑150 (MQ=255)
                                                                                                                     tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGcc                               >  2:1018709/1‑150 (MQ=255)
                                                                                                                     tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGcc                               >  2:5304187/1‑150 (MQ=255)
                                                                                                                          cACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGttt                          <  1:1985705/150‑1 (MQ=255)
                                                                                                                           aCGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtttt                         >  1:1035724/1‑150 (MQ=255)
                                                                                                                           aCGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtttt                         >  1:293635/1‑150 (MQ=255)
                                                                                                                            cGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTg                        <  2:2881221/150‑1 (MQ=255)
                                                                                                                               aGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGCTCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcg                     <  2:2118930/150‑1 (MQ=255)
                                                                                                                                 tCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcgcg                   <  1:209008/150‑1 (MQ=255)
                                                                                                                                   ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcggg                 >  2:317079/1‑150 (MQ=255)
                                                                                                                                       ggCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccgg             <  1:3362106/150‑1 (MQ=255)
                                                                                                                                       ggCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccgg             <  2:4633744/150‑1 (MQ=255)
                                                                                                                                        gCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtttt                         >  1:3554550/1‑137 (MQ=255)
                                                                                                                                        gCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtttt                         <  2:3554550/137‑1 (MQ=255)
                                                                                                                                         aagccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccg              >  1:4098908/4‑147 (MQ=255)
                                                                                                                                         aagccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccg              <  2:4098908/144‑1 (MQ=255)
                                                                                                                                          ttCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCt          >  2:2346770/1‑150 (MQ=255)
                                                                                                                                           gccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCGTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCtt         >  2:3471383/2‑150 (MQ=255)
                                                                                                                                            ccACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTc        <  1:3103384/150‑1 (MQ=255)
                                                                                                                                            ccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTc        >  1:1680374/1‑150 (MQ=255)
                                                                                                                                             cacaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTcc       <  1:3691365/150‑1 (MQ=255)
                                                                                                                                              acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcgcg                   <  2:687308/137‑1 (MQ=255)
                                                                                                                                              acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg      <  1:2738312/150‑1 (MQ=255)
                                                                                                                                              acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg      >  2:1717887/1‑150 (MQ=255)
                                                                                                                                              acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg      >  1:3879919/1‑150 (MQ=255)
                                                                                                                                              acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg      >  1:3880372/1‑150 (MQ=255)
                                                                                                                                               caGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                       <  2:3287999/116‑1 (MQ=255)
                                                                                                                                               caGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                       >  1:3287999/1‑116 (MQ=255)
                                                                                                                                               caGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGACATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCgg     <  2:5243927/150‑1 (MQ=255)
                                                                                                                                                aGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAgg                                                                       >  1:2389069/1‑83 (MQ=255)
                                                                                                                                                aGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAgg                                                                       <  2:2389069/83‑1 (MQ=255)
                                                                                                                                                  gCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGt  >  1:668714/1‑150 (MQ=255)
                                                                                                                                                     ggCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccg              >  1:5194757/1‑135 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGT  >  NC_002947/1481374‑1481669

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