Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,037 A→C 5.2% intergenic (‑143/+41) 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,0370AC 5.2% 245.5 / 5.3 77intergenic (‑143/+41)yiiX/metJputative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base C (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGC  >  NC_000913/4127892‑4128183
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                                                 >  7:395837/1‑148 (MQ=255)
 aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                                 >  7:401260/1‑148 (MQ=255)
 aTAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgcg                                                                                                                                                 >  4:286016/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:47586/149‑1 (MQ=255)
                             gTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGctct                                                                                                                    <  1:105183/149‑1 (MQ=255)
                                  tCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                                  <  8:392332/114‑1 (MQ=255)
                                  tCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAgc                                                                                                                                                  >  7:392332/1‑114 (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                                                                                                                 <  8:43629/142‑1 (MQ=255)
                                       gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtttt                                                                                                                 >  7:43629/1‑142 (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                                                                                                  >  1:134002/1‑134 (MQ=255)
                                                              cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTcc                                                                                                  <  2:134002/134‑1 (MQ=255)
                                                               ggATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTaa                                                                                  <  8:44635/149‑1 (MQ=255)
                                                                 aTGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATc                                                                                <  4:200750/149‑1 (MQ=255)
                                                                     ggcgTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCa                                                                            <  5:28409/149‑1 (MQ=255)
                                                                        gTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa                                                                                                         >  2:291614/1‑117 (MQ=255)
                                                                        gTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa                                                                                                         <  1:291614/117‑1 (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                                       <  3:124488/149‑1 (MQ=255)
                                                                                                          aTTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTtc                                       >  6:93484/1‑149 (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)
                                                                                                                                              cGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCTTCGCTGCGCTCTTTACGCAGATCGGCACCATCCgg   >  6:214062/1‑149 (MQ=255)
                                                                                                                                               gTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCgg   >  7:63269/1‑148 (MQ=255)
                                                                                                                                                tAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGc  >  5:284100/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGC  >  NC_000913/4127892‑4128183

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: