Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,334,215 C→G 6.2% intergenic (+28/‑21) dksA → / → gluQ RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,334,2150CG 6.2% 530.8 / 19.1 160intergenic (+28/‑21)dksA/gluQRNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (80/70);  new base G (3/7);  total (83/77)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 167 total aligned reads displayed.

GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGG  >  NC_002947/5334066‑5334363
                                                                                                                                                     |                                                                                                                                                    
ggTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGcatc                                                                                                                                                      >  1:3470633/1‑150 (MQ=255)
    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                                  >  2:3257124/1‑150 (MQ=255)
         cGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGcccc                                                                                                                                             <  1:2514789/150‑1 (MQ=255)
         cGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGcccc                                                                                                                                             <  1:3579963/150‑1 (MQ=255)
          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                            <  2:2872731/150‑1 (MQ=255)
          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                            <  1:1585501/150‑1 (MQ=255)
             cGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGCCCCGttc                                                                                                                                         <  2:3918600/150‑2 (MQ=255)
             cGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGCCCCGttc                                                                                                                                         <  2:3003859/150‑2 (MQ=255)
              ggCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTc                                                                                                                                        >  2:1715182/1‑150 (MQ=255)
               gCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCa                                                                                                                                       <  2:4717860/150‑1 (MQ=255)
               gCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTAAACGGGGCGCATGATGCGCCCCGttcag                                                                                                                                       <  1:1884768/150‑4 (MQ=255)
                cATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                      >  2:3050177/1‑150 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                     >  1:2690157/1‑150 (MQ=255)
                 aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                     >  2:4647330/1‑150 (MQ=255)
                  tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                                    <  2:1846170/150‑1 (MQ=255)
                  tCGAGATCGGCCGGCGCCGCCTGGAAGCCCGCCCTCCCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAACAGGAAAAAACGGTCGGCAAAGGCGGATCTTTCGTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                                    <  2:3663975/150‑1 (MQ=255)
                    gagaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtat                                                                                                                                  >  1:1782191/1‑150 (MQ=255)
                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTCCCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                                 <  1:1847250/150‑1 (MQ=255)
                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                                 <  1:1847265/150‑1 (MQ=255)
                      gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTatat                                                                                                                                <  2:2419191/150‑1 (MQ=255)
                      gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTatat                                                                                                                                <  1:277311/150‑1 (MQ=255)
                          ggCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGcatc                                                                                                                                                      <  1:3263953/124‑1 (MQ=255)
                           gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                                           >  2:572308/1‑150 (MQ=255)
                             cTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGc                                                                                                                         <  2:4977656/150‑1 (MQ=255)
                                  ccgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                                    <  1:572308/150‑1 (MQ=255)
                                  ccgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                                    <  1:1263679/150‑1 (MQ=255)
                                   cgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAcc                                                                                                                   <  2:3838771/150‑1 (MQ=255)
                                    accTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCa                                                                                                                  >  1:3407259/2‑150 (MQ=255)
                                     ccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                                 >  1:1263633/1‑150 (MQ=255)
                                        ggAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCTAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAc                                                                                                              <  1:4712035/150‑1 (MQ=255)
                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCACGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCCGCCTGACCATGAcc                                                                                                             >  2:4896220/1‑150 (MQ=255)
                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                             >  1:2191612/1‑150 (MQ=255)
                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                             >  2:3267247/1‑150 (MQ=255)
                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                             >  1:4419033/1‑150 (MQ=255)
                                          aaGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCg                                                                                                            >  2:3312578/1‑150 (MQ=255)
                                            gcccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGAc                                                                                                          >  1:2812682/1‑150 (MQ=255)
                                                gcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCCTCATGCGCCCCGTTTCAGTTATATGCCCAGCCTGACCATGACCGACTCCa                                                                                                      >  1:1363326/1‑150 (MQ=255)
                                                 cccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAg                                                                                                     >  2:4354229/1‑150 (MQ=255)
                                                   cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                                                   <  1:1793469/150‑1 (MQ=255)
                                                   cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                                                   <  2:697531/150‑1 (MQ=255)
                                                    tACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTa                                                                                                  >  1:3731011/1‑150 (MQ=255)
                                                            aCCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgc                                                                                          >  2:2660020/1‑150 (MQ=255)
                                                                gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCtt                                                                                      >  1:415695/1‑150 (MQ=255)
                                                                 tgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                     <  2:198060/150‑1 (MQ=255)
                                                                  gCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCg                                                                                    >  2:5113273/1‑150 (MQ=255)
                                                                  gCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCg                                                                                    >  2:3645315/1‑150 (MQ=255)
                                                                    ttCGACTGCAAGTAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCATAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCGGACCATGCCCGACTCCAGCTACATCGGGCGCTTCGcc                                                                                  <  1:1967171/150‑1 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                <  2:2590461/150‑1 (MQ=255)
                                                                      cGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccc                                                                                <  1:5113273/150‑1 (MQ=255)
                                                                       gACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccccc                                                                               >  1:1323137/1‑150 (MQ=255)
                                                                       gACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccccc                                                                               >  1:1323160/1‑150 (MQ=255)
                                                                          tGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                                 >  2:1334520/1‑113 (MQ=255)
                                                                          tGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                                 <  1:1334520/113‑1 (MQ=255)
                                                                          tGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccacc                                                                            >  1:1188565/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                         >  1:625462/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                         >  2:2512136/1‑150 (MQ=255)
                                                                             aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccc                                                                         >  1:4656561/1‑150 (MQ=255)
                                                                               ggAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAg                                                                       >  2:2306511/1‑150 (MQ=255)
                                                                               ggAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAg                                                                       <  1:1644288/150‑1 (MQ=255)
                                                                                gAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGc                                                                      <  2:1501931/150‑1 (MQ=255)
                                                                                gAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGc                                                                      <  2:1501943/150‑1 (MQ=255)
                                                                                 aaaTTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCg                                                                     <  1:3257124/150‑1 (MQ=255)
                                                                                   aTTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGc                                                                   >  1:984494/1‑150 (MQ=255)
                                                                                     tGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCtt                                                                 <  2:1611698/150‑1 (MQ=255)
                                                                                       ccGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTcc                                                               >  1:4818373/1‑150 (MQ=255)
                                                                                       ccGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTcc                                                               >  1:4961071/1‑150 (MQ=255)
                                                                                       ccGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              >  2:3280091/1‑150 (MQ=255)
                                                                                        cGATAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              <  2:2811038/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              <  1:3498211/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              <  1:1080382/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              <  2:1991768/150‑1 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              >  1:2993305/1‑150 (MQ=255)
                                                                                        cGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGGTTCATATATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                              >  1:5055038/1‑150 (MQ=255)
                                                                                         gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                             >  2:2722896/1‑150 (MQ=255)
                                                                                              aaaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACtt                                                        >  1:1945040/1‑150 (MQ=255)
                                                                                               aaGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTc                                                       <  1:4260801/150‑1 (MQ=255)
                                                                                                 ggAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCgg                                                     <  1:4031946/150‑1 (MQ=255)
                                                                                                           gTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTc                                           <  2:1237628/150‑1 (MQ=255)
                                                                                                              ggCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgcc                                        <  2:383247/150‑1 (MQ=255)
                                                                                                                cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGGGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc                                      >  1:1020254/1‑150 (MQ=255)
                                                                                                                cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc                                      <  2:2666236/150‑1 (MQ=255)
                                                                                                                     gCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCg                                 >  2:2287735/1‑150 (MQ=255)
                                                                                                                      cTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGc                                <  2:1734605/150‑1 (MQ=255)
                                                                                                                      cTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGc                                <  2:1782191/150‑1 (MQ=255)
                                                                                                                       tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                               >  1:1941112/1‑150 (MQ=255)
                                                                                                                         aTCTTTCTTCCACCTGAACGGGGCGCGTCATGCGCCCCGTTTCATTTATTTGCCCAGCCTTACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                             <  1:4572618/150‑1 (MQ=255)
                                                                                                                         aTCTTTCTTCCACCTGAACGGGGCGCGTCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                             <  1:4572990/150‑1 (MQ=255)
                                                                                                                         aTCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                             <  1:4354229/150‑1 (MQ=255)
                                                                                                                           ctttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctcct                           >  2:115869/1‑150 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTg                          >  1:3842728/1‑150 (MQ=255)
                                                                                                                            tttcttCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTg                          >  2:1918619/1‑150 (MQ=255)
                                                                                                                             ttcttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAAGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTgg                         >  1:1724956/1‑150 (MQ=255)
                                                                                                                                ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc                      >  2:2904179/1‑150 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                   <  1:2566696/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                   <  1:2817362/150‑1 (MQ=255)
                                                                                                                                   cACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                   >  2:4983823/1‑150 (MQ=255)
                                                                                                                                        gAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcg              >  1:3085121/1‑150 (MQ=255)
                                                                                                                                           cGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCgcgc           >  1:1801864/1‑150 (MQ=255)
                                                                                                                                              ggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGc                                                    >  2:4404158/1‑106 (MQ=255)
                                                                                                                                                    tcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGACCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACgg  <  2:3731011/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAACGG  >  NC_002947/5334066‑5334363

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: