Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,032 A→C 5.3% intergenic (‑138/+46) yiiX ← / ← metJ putative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,0320AC 5.3% 224.5 / 5.4 75intergenic (‑138/+46)yiiX/metJputative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base C (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

TGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATC  >  NC_000913/4127885‑4128179
                                                                                                                                                   |                                                                                                                                                   
tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAg                                                                                                                                                    >  6:135857/1‑149 (MQ=255)
tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAg                                                                                                                                                    >  6:18085/1‑149 (MQ=255)
tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAg                                                                                                                                                    >  2:47586/1‑149 (MQ=255)
       cATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                             >  6:385017/1‑149 (MQ=255)
       cATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                              <  1:334062/148‑1 (MQ=255)
       cATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                              <  1:345093/148‑1 (MQ=255)
       cATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                              <  2:291627/148‑1 (MQ=255)
        aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                             >  6:73020/1‑148 (MQ=255)
        aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                             >  4:286016/1‑148 (MQ=255)
        aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                             >  7:395837/1‑148 (MQ=255)
        aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                             >  7:401260/1‑148 (MQ=255)
            ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCat                                                                                                                                          <  7:44635‑M2/147‑3 (MQ=255)
                aCGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCatcaggc                                                                                                                                     <  8:109843‑M2/148‑8 (MQ=255)
                     cATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTc                                                                                                                                >  6:251208/1‑148 (MQ=255)
                         tataAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTc                                                                                                                           >  4:322011/1‑149 (MQ=255)
                            aaaaaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                                                                                                                        >  4:67462/1‑149 (MQ=255)
                               aaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATGCGCTACGCTTatcaggcctacgtcatattgcaa                                                                                                                     <  3:305174‑M2/149‑24 (MQ=255)
                                    gTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGctct                                                                                                                <  1:105183/149‑1 (MQ=255)
                                    gTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGctct                                                                                                                <  1:47586/149‑1 (MQ=255)
                                         tCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                              >  7:392332/1‑114 (MQ=255)
                                         tCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                              <  8:392332/114‑1 (MQ=255)
                                         tCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtttttt                                                                                                           >  6:28409/1‑149 (MQ=255)
                                           gCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATGCGCTACGCTTatcaggcctacgtcatattgcaatttattgaattt                                                                                                         <  3:229853‑M2/149‑36 (MQ=255)
                                           gCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTg                                                                                                          <  4:178361/148‑1 (MQ=255)
                                            cAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt                                                                                                        <  2:84501/149‑1 (MQ=255)
                                            cAGACAACGTGCTCGTTGTTTATGCCGGATCCGGCGTGAACGCCATATCCGCCCGACAAGTTCGTGCAAATTCAATAAATAGCAAAATGACGTACGCCTGATAAGAGTAGCGCATCAGGCGATTc                                                                                                                                >  5:155128/1‑125 (MQ=255)
                                             aGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCtt                                                                                                       <  4:116883/149‑1 (MQ=255)
                                              gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtttt                                                                                                             >  7:43629/1‑142 (MQ=255)
                                              gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtttt                                                                                                             <  8:43629/142‑1 (MQ=255)
                                               aCAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATGCGCTACGCTTatcaggcctacgtcatattgcaatttattgaatttgcac                                                                                                     >  4:305174‑M2/1‑110 (MQ=255)
                                               aCAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa                                                                                                     <  5:385017/149‑1 (MQ=255)
                                               aCAACGTGCTCGTTGTTTATGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGACGCTCTTTTTTGCTTTa                                                                                                     <  4:35955/149‑1 (MQ=255)
                                                       cTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTccc                                                                                             >  5:164317/1‑149 (MQ=255)
                                                          gttgttTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCAc                                                                                           <  5:18085/148‑1 (MQ=255)
                                                          gttgttTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACg                                                                                          >  5:307264/1‑149 (MQ=255)
                                                                 aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCg                                                                                    <  3:286016/148‑1 (MQ=255)
                                                                  tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCggg                                                                                  <  5:251208/149‑1 (MQ=255)
                                                                   gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCCTCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGt                                                                                 >  5:114780/1‑149 (MQ=255)
                                                                    ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATGCGCTACGCTTatcaggcctacgtcatattgcaatttattgaatttgcacgaacttgtaggccggataag                                                                                 >  4:229853‑M2/1‑89 (MQ=255)
                                                                     cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTcc                                                                                              <  2:134002/134‑1 (MQ=255)
                                                                     cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTcc                                                                                              >  1:134002/1‑134 (MQ=255)
                                                                      ggATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTaa                                                                              <  8:44635/149‑1 (MQ=255)
                                                                        aTGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATc                                                                            <  4:200750/149‑1 (MQ=255)
                                                                            ggcgTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCa                                                                        <  5:28409/149‑1 (MQ=255)
                                                                               gTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa                                                                                                     <  1:291614/117‑1 (MQ=255)
                                                                               gTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa                                                                                                     >  2:291614/1‑117 (MQ=255)
                                                                                tGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCAtctc                                                                    <  5:292861/149‑1 (MQ=255)
                                                                                  aaCGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCAc                                                                  <  6:307264/149‑1 (MQ=255)
                                                                                       cTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCACATCTCACGCAt                                                              >  4:156435/1‑148 (MQ=255)
                                                                                          aTCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCCCACGCATGAt                                                           <  5:135857/148‑1 (MQ=255)
                                                                                           tCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGAtct                                                         >  8:90303/1‑149 (MQ=255)
                                                                                             cGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATctct                                                       <  5:135637/149‑1 (MQ=255)
                                                                                                ccTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCtttt                                                    <  7:109843/149‑1 (MQ=255)
                                                                                                  tACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGc                                                  >  7:383293/1‑149 (MQ=255)
                                                                                                       gTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCtt                                                                                                       >  8:224579/1‑91 (MQ=255)
                                                                                                           gTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCgg                                         >  2:57196/1‑149 (MQ=255)
                                                                                                           gTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCg                                          <  3:322011/148‑1 (MQ=255)
                                                                                                            tGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCgg                                         >  7:405310/1‑148 (MQ=255)
                                                                                                                 aTTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTtc                                   >  6:93484/1‑149 (MQ=255)
                                                                                                                 aTTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTtc                                   <  3:124488/149‑1 (MQ=255)
                                                                                                                  ttCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTtcg                                  >  6:71385/1‑149 (MQ=255)
                                                                                                                     aaTAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCCCGCATGCTCTCTTTTGCCGCTTCCGGGATTtcgtcg                               >  7:337129/1‑149 (MQ=255)
                                                                                                                       tAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCt                             <  5:350221/149‑1 (MQ=255)
                                                                                                                         aaTTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTg                            <  7:163905/148‑1 (MQ=255)
                                                                                                                                atatGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCttt                    <  7:154573/149‑1 (MQ=255)
                                                                                                                                     aCGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGc                <  4:291221/148‑1 (MQ=255)
                                                                                                                                     aCGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGGTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTACGGGATTTCGTCGCTGCGCTCTTTACg                 >  5:304766/1‑147 (MQ=255)
                                                                                                                                         aGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCCCTATTTTTTGCTTTAGTAtt                                                                                                >  2:131949/1‑64 (MQ=255)
                                                                                                                                          ggCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCg          <  8:405310/149‑1 (MQ=255)
                                                                                                                                          ggCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCACCATCTCTCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATc           <  5:311255/148‑1 (MQ=255)
                                                                                                                                           gCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCgg         >  6:348944/1‑149 (MQ=255)
                                                                                                                                               gATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCtt                                                      <  7:43379/100‑1 (MQ=255)
                                                                                                                                               gATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCtt                                                      >  8:43379/1‑100 (MQ=255)
                                                                                                                                                aTAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGcatca    >  4:162865/1‑149 (MQ=255)
                                                                                                                                                aTAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGcatc     >  8:133198/1‑148 (MQ=255)
                                                                                                                                                  aaGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCatcatc  >  3:330143/1‑149 (MQ=255)
                                                                                                                                                   aGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCatcatc  <  6:164317/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATC  >  NC_000913/4127885‑4128179

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