Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,704 G→C 6.5% intergenic (+23/‑205) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7040GC 6.5% 345.1 / 8.9 113intergenic (+23/‑205)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (41/63);  new base C (6/2);  total (48/65)
Fisher's exact test for biased strand distribution p-value = 6.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

ATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGA  >  NC_000913/3334555‑3334852
                                                                                                                                                     |                                                                                                                                                    
aTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGCCACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCTCTGTTCCACGCGCTCGTTAATCCCCTCCCCTCAACCCGAGCAg                                                                                                                                                      >  1:1077992/1‑150 (MQ=255)
aTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                      >  2:1587858/1‑150 (MQ=255)
    aCGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCcgcgcg                                                                                                                                                  >  2:107270/1‑150 (MQ=255)
          cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGATGCACTCATCGGCCTCGCGCACATCGCTGTTAACGCGATCGTTAATCGCCTCCCCTCATCCGGCGCGCTGCGCGGGATg                                                                                                                                           <  2:464734/150‑1 (MQ=255)
             cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGt                                                                                                                                         <  2:1443086/150‑1 (MQ=255)
                gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                      <  1:229394/150‑1 (MQ=255)
                gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                      <  1:1071473/150‑1 (MQ=255)
                gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                      <  1:570094/150‑1 (MQ=255)
                gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                      <  1:1616820/150‑1 (MQ=255)
                gCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                      <  2:17753/150‑1 (MQ=255)
                  cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAg                                                                                                                                    <  2:1453110/150‑1 (MQ=255)
                  cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAg                                                                                                                                    <  2:1510304/150‑1 (MQ=255)
                  cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAg                                                                                                                                    <  1:795029/150‑1 (MQ=255)
                   gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGATGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGt                                                                                                                                   <  1:1261996/150‑1 (MQ=255)
                    aGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGCTGGTTCCAgga                                                                                                                                  >  1:1316864/1‑148 (MQ=255)
                     ggAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                 >  1:1301953/1‑150 (MQ=255)
                     ggAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                 >  2:940807/1‑150 (MQ=255)
                         gaagCAGACAAAGCCATCGCAGCGTTACAGGGGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTc                                                                                                                             <  2:345308/150‑1 (MQ=255)
                           agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                           >  1:1054536/1‑150 (MQ=255)
                             cAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACGCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCACCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGGTCCAGTAAGTTCCATa                                                                                                                         >  2:25239/1‑150 (MQ=255)
                              aGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGGCGCGCGGGGTGGTTCCAGTTAGTTCCATaa                                                                                                                        >  1:1436561/1‑150 (MQ=255)
                                ccAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                      <  1:1657726/149‑1 (MQ=255)
                                 cAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                     >  1:714649/1‑150 (MQ=255)
                                 cAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                     >  2:1289084/1‑150 (MQ=255)
                                  aaaGCCATCGCAGCGTTACATGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAac                                                                                                                    <  1:1130984/150‑1 (MQ=255)
                                  aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAac                                                                                                                    <  1:1299510/150‑1 (MQ=255)
                                    aGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAacac                                                                                                                  >  2:1695080/1‑150 (MQ=255)
                                      ccATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTACTCCCCTCCCCTCATCCCGCGCCGCCCGCGGGATGGTGCCAGTAAGTTCCAGAAAAACACtt                                                                                                                >  2:271754/1‑150 (MQ=255)
                                          cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGGTCCAGTAAGGTCCATAAAAACGCTTATTc                                                                                                            >  1:1594443/1‑150 (MQ=255)
                                           gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCa                                                                                                           <  2:1523104/150‑1 (MQ=255)
                                                  tACAGGTGCTCCCGGACACCCCTTGGCTAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                                    <  1:899080/150‑1 (MQ=255)
                                                  tACAGGTGCTCCCGGAACACCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATAGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                                    <  1:44976/150‑1 (MQ=255)
                                                    cAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                  >  1:47558/1‑150 (MQ=255)
                                                      ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAaca                                                                                                >  2:12099/1‑150 (MQ=255)
                                                          cTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATACCGCGAATCGCTCGGGCTGGGTCCAGTTAATTCCATAAAAACACATAATAACCTCTATAacaacag                                                                                            >  1:1318009/1‑146 (MQ=255)
                                                           tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGa                                                                                           >  1:276909/1‑150 (MQ=255)
                                                            cccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCCCGCACCGCGCGGGATTGTTCCAGTAACTTCCAAAAAAAAACTTAATAACCTCAATAACACCTGaa                                                                                          >  2:802966/1‑150 (MQ=255)
                                                                    ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAggg                                                                                  <  2:966671/150‑1 (MQ=255)
                                                                     cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                 <  2:735507/150‑1 (MQ=255)
                                                                       cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAATAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                               >  1:1268526/1‑150 (MQ=255)
                                                                       cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGGTTATCCGCTCCCCTCATCCCGCGCAGCGCGCGGGAAGGTTCCAGTAAGGTCCATTAACACCCTTATTCAGCTCTAGCACCCCTGAAAAATAGGGagg                                                                               >  1:1220564/1‑148 (MQ=255)
                                                                         tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                             <  2:714649/150‑1 (MQ=255)
                                                                         tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                             <  1:1212797/150‑1 (MQ=255)
                                                                          ggCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGc                                                                            >  2:26101/1‑150 (MQ=255)
                                                                           gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCa                                                                           <  2:479930/150‑1 (MQ=255)
                                                                                aaGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCttt                                                                      <  2:88195/150‑1 (MQ=255)
                                                                                 aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                     >  2:1571718/1‑150 (MQ=255)
                                                                                 aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                     <  1:1411894/150‑1 (MQ=255)
                                                                                 aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                     <  2:166754/150‑1 (MQ=255)
                                                                                 aGCACTCATCGGCCTCGCGCACATCGCGGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCAGCTTta                                                                     <  1:107270/150‑1 (MQ=255)
                                                                                   cACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTGCAGTAAGTTCCATAAAAACACTTATTCAGCGCTTTAACACCTGAAAAATAGGGATGTGGATCTTtata                                                                   >  2:854239/1‑150 (MQ=255)
                                                                                   cACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCACCACTATAACACCTGAAAAATAGGGATGTGCATCTTtata                                                                   >  1:286501/1‑150 (MQ=255)
                                                                                   cACTCATCGGCCTCGAGCACATCGCTGTTCAACGCGATAGGTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTAATTAAGCTCTATAACAACTGAAAAATAGGGATGTGCATCTTtata                                                                   >  1:54713/1‑150 (MQ=255)
                                                                                       cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGGTGGTTCCAGTAAGTTGCATAAAAACACGTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatttttt                                                               >  1:1471955/1‑145 (MQ=255)
                                                                                         tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCt                                                             <  2:14386/150‑1 (MQ=255)
                                                                                         tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCt                                                             <  2:869090/150‑1 (MQ=255)
                                                                                          cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATtgtg                                                            >  1:133042/1‑150 (MQ=255)
                                                                                           ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTATGa                                                           >  2:1706389/1‑150 (MQ=255)
                                                                                           ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTAAAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                           >  2:391373/1‑150 (MQ=255)
                                                                                           ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTACAGTAAGTTCCATAAAAACACTTATTCATCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATTTATTCTGa                                                           >  2:1826828/1‑150 (MQ=255)
                                                                                            gCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                          <  1:908783/150‑1 (MQ=255)
                                                                                                 gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                     <  1:1008812/150‑1 (MQ=255)
                                                                                                  cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATAttaa                                                    >  1:770353/1‑148 (MQ=255)
                                                                                                    cacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                  >  2:50526/1‑150 (MQ=255)
                                                                                                         cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCtt                                             <  1:445386/150‑1 (MQ=255)
                                                                                                          gCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTCCAGGAAGTTCCATAAAAACACTTATTCAGCGGTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                            >  2:907901/1‑150 (MQ=255)
                                                                                                           cTGTTCAACGCGATCGTTAATCCGCTCCCCTCATCCCGCGCAGCGCTCTGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACGTGATAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                           <  2:486274/150‑1 (MQ=255)
                                                                                                            tGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCTCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTAc                                          <  1:1624772/150‑1 (MQ=255)
                                                                                                             gTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                         <  1:1340813/150‑1 (MQ=255)
                                                                                                             gTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                         >  2:659639/1‑150 (MQ=255)
                                                                                                              ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                        <  1:1121404/150‑1 (MQ=255)
                                                                                                              ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                        <  1:406775/150‑1 (MQ=255)
                                                                                                               tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAgg                                       >  1:1086654/1‑150 (MQ=255)
                                                                                                               tCAACGCGATCGTTAATCCCCTCCCCTCATAGCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAgg                                       <  2:125490/150‑1 (MQ=255)
                                                                                                                cAACGCGATCGTTAATCCCCTCCACTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGa                                      <  1:1155558/150‑1 (MQ=255)
                                                                                                                 aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                     >  1:47845/1‑150 (MQ=255)
                                                                                                                 aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATACTCTGAATATTCACACTCTTTACAGGaa                                     <  1:804591/150‑1 (MQ=255)
                                                                                                                       ttCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCATTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                               <  2:1597648/149‑1 (MQ=255)
                                                                                                                       aTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                               <  1:1587858/150‑1 (MQ=255)
                                                                                                                        tCGTTAATCCCCTCCCCTCATCCCGCGCAGCGAGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACCCTTTACAGGAACTTTTTa                              <  1:271754/150‑1 (MQ=255)
                                                                                                                         cggtAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGCGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                             <  2:260402/147‑1 (MQ=255)
                                                                                                                         cGTTAATCCCCTCACCTCATACAGAGCAGCGCGCGGGATTGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttttg                             <  1:409351/150‑3 (MQ=255)
                                                                                                                           ttAATCCCCTCCCTTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATATTTATATATTCTGAATATTCACACTCTTAACAGGAACTTTTTagag                           >  1:215393/1‑150 (MQ=255)
                                                                                                                            tAATCCCCTCCCCTCATGCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGc                          <  1:908047/150‑1 (MQ=255)
                                                                                                                            tAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGc                          <  2:1054536/150‑1 (MQ=255)
                                                                                                                                 ccctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                     >  1:306700/1‑150 (MQ=255)
                                                                                                                                  cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                    <  2:1301953/150‑1 (MQ=255)
                                                                                                                                   ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGc                   <  1:1485146/150‑1 (MQ=255)
                                                                                                                                      ccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTAACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCAt                >  2:1319130/1‑150 (MQ=255)
                                                                                                                                         tcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg             <  1:940807/150‑1 (MQ=255)
                                                                                                                                         ccATCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATATGCCATCAg             >  2:1015710/2‑150 (MQ=255)
                                                                                                                                         ccATCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCGCTCAg             >  2:1199070/2‑150 (MQ=255)
                                                                                                                                           aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATTTATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCTTCAggg           >  1:1103361/1‑149 (MQ=255)
                                                                                                                                           aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTATGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa           >  2:74221/1‑150 (MQ=255)
                                                                                                                                           aTCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCGCGATAACACCTGAAAAATAGGGATGTGCATCTTTATTTATTCTGATTATTCACAATGTTTACAGGACCTTTTTATCGCAATTGGCCATCAGGa           >  2:217335/1‑150 (MQ=255)
                                                                                                                                             cccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGt         >  1:1571812/1‑150 (MQ=255)
                                                                                                                                               cgcgCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGtat       >  2:174856/1‑150 (MQ=255)
                                                                                                                                                 cgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAg     <  2:726183/150‑1 (MQ=255)
                                                                                                                                                   cAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTAGATAACACCAGAAAAATAGGGATGTGCATCTTTATAAATTCTGAATTTTCACACTCTTTACAGGATCATTTAAGCGCCATAGGCCATCAGGAGTAGAGTg   >  1:1052209/1‑150 (MQ=255)
                                                                                                                                                    aGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGa  <  1:1695080/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGA  >  NC_000913/3334555‑3334852

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