Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,877,558 G→C 5.3% intergenic (‑93/+147) yidB ← / ← gyrB DUF937 family protein/DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,5580GC 5.3% 594.5 / 18.8 170intergenic (‑93/+147)yidB/gyrBDUF937 family protein/DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (85/76);  new base C (6/3);  total (91/79)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 176 total aligned reads displayed.

TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATT  >  NC_000913/3877410‑3877706
                                                                                                                                                    |                                                                                                                                                    
 atTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAgatga                                                                                                                                                     <  4:1635381/149‑1 (MQ=255)
   ttCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGAtc                                                                                                                                                   <  4:308287/149‑1 (MQ=255)
    tCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGAtct                                                                                                                                                  <  3:1606742/149‑1 (MQ=255)
     ccccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGAtctc                                                                                                                                                 <  1:1632992/149‑1 (MQ=255)
       ccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCg                                                                                                                                               >  4:840030/1‑149 (MQ=255)
       ccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGCTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCg                                                                                                                                               >  3:1201749/1‑149 (MQ=255)
          cgtCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATATGCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGccc                                                                                                                                            <  2:1668052/147‑1 (MQ=255)
            tCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCtt                                                                                                                                          >  3:63344/1‑149 (MQ=255)
              gCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCtttt                                                                                                                                        <  2:1214507/149‑1 (MQ=255)
               ccTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTc                                                                                                                                       <  2:861270/149‑1 (MQ=255)
                cTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCt                                                                                                                                      >  2:83374/1‑149 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                     >  1:1117405/1‑149 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                     >  3:1351381/1‑149 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                     >  2:1138115/1‑149 (MQ=255)
                    cAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTtata                                                                                                                                  <  4:1348697/149‑1 (MQ=255)
                     aGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTc                                                                                                                                       <  3:651715/143‑1 (MQ=255)
                             aCCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGc                                                                                                                         >  4:1323163/1‑149 (MQ=255)
                               cGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                    >  1:1541442/1‑136 (MQ=255)
                               cGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCg                                                                                                                       >  1:571143/1‑149 (MQ=255)
                                gACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCgg                                                                                                                      <  1:341455/149‑1 (MQ=255)
                                 aCAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGa                                                                                                                     >  1:987445/1‑149 (MQ=255)
                                   aaCTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATg                                                                                                                   >  3:1529656/1‑149 (MQ=255)
                                    aCTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgc                                                                                                                  >  3:1660537/1‑149 (MQ=255)
                                     cTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcg                                                                                                                 >  3:698343/1‑149 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcgg                                                                                                                >  2:1465570/1‑149 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcgg                                                                                                                >  3:48396/1‑149 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcgg                                                                                                                >  3:451662/1‑149 (MQ=255)
                                         atcaAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGt                                                                                                             <  4:699582/149‑1 (MQ=255)
                                          tcaAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTAAAGCTTCTTGCCGGATGCGGCGTg                                                                                                            >  4:496784/1‑149 (MQ=255)
                                            aaaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGaa                                                                                                          >  1:898642/1‑149 (MQ=255)
                                            aaaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGGAGAGTTTTGTGGGATGATAATTGCAGATTGCAGGAATAAAAAGGGCCGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCGGGATGCGGCGTGaa                                                                                                          >  4:1641587/1‑149 (MQ=255)
                                             aaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAAc                                                                                                         <  3:464857/149‑1 (MQ=255)
                                             aaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAAc                                                                                                         <  3:981982/149‑1 (MQ=255)
                                                   cccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGAtctc                                                                                                                                                 >  2:1495112/1‑103 (MQ=255)
                                                   cccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGAtctc                                                                                                                                                 <  1:1495112/103‑1 (MQ=255)
                                                   cccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTtat                                                                                                                                   <  2:735667/117‑1 (MQ=255)
                                                   cccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                   <  4:706746/149‑1 (MQ=255)
                                                    ccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                  <  4:639247/149‑1 (MQ=255)
                                                    ccATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGCCCttt                                                                                                                                         <  2:211859/110‑1 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 >  1:325135/1‑149 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 >  3:913724/1‑149 (MQ=255)
                                                      aTGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                >  2:524285/1‑149 (MQ=255)
                                                       tGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                               >  3:479584/1‑149 (MQ=255)
                                                               aCCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                    <  1:809749/104‑1 (MQ=255)
                                                                ccTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                                      >  4:68543/1‑149 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGGTGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCTGATGCGGCGTGAACGCCTTATCCGGCCCACaaaa                                                                                     >  4:1531801/1‑149 (MQ=255)
                                                                  tGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                              >  2:821338/1‑139 (MQ=255)
                                                                  tGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAt                                                                                    >  3:411625/1‑149 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCTGAAGTAAAAAGTGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  <  2:1251579/149‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  <  3:1449946/149‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  <  3:1331576/149‑1 (MQ=255)
                                                                     aaaGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                 <  4:451662/149‑1 (MQ=255)
                                                                     aaaGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCAt                                                                                 <  1:1015151/149‑1 (MQ=255)
                                                                      aaGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATg                                                                                <  4:913724/149‑1 (MQ=255)
                                                                         cGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAAc                                                                                                         <  1:1711008/121‑1 (MQ=255)
                                                                               gtgtTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                                                       <  4:1529656/149‑1 (MQ=255)
                                                                                     aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                            <  3:1623478/122‑1 (MQ=255)
                                                                                     aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                 >  3:281477/1‑149 (MQ=255)
                                                                                       ccAAAAAGTATAGATTTGAGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGc                                                               >  4:257490/1‑149 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCaa                                                                       >  2:1683576/1‑139 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                                                             >  1:1291502/1‑149 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                                                             <  3:1592566/149‑1 (MQ=255)
                                                                                         aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCaa                                                             <  2:1604297/149‑1 (MQ=255)
                                                                                              gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAttt                                                        <  1:1137495/149‑1 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                                                       >  4:1018980/1‑149 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                                                       >  1:290922/1‑149 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                                                       <  3:819067/149‑1 (MQ=255)
                                                                                                ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                      <  1:1465570/149‑1 (MQ=255)
                                                                                                ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                      <  4:1201749/149‑1 (MQ=255)
                                                                                                ataGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTc                                                      <  2:538303/149‑1 (MQ=255)
                                                                                                   gATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATGCGGCCTACAAAATCATGAAAATTCAATACCTTGCAAGATTTTCGTa                                                   <  2:1096179/149‑1 (MQ=255)
                                                                                                   gATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                                                   <  3:688332/149‑1 (MQ=255)
                                                                                                   gATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                                                   >  2:1408129/1‑149 (MQ=255)
                                                                                                   gATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTa                                                   <  1:202743/149‑1 (MQ=255)
                                                                                                    aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAg                                                  <  3:772197/149‑1 (MQ=255)
                                                                                                          gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                            <  3:1018980/149‑1 (MQ=255)
                                                                                                          gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                            <  2:290922/149‑1 (MQ=255)
                                                                                                           ggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                     >  1:1460803/1‑107 (MQ=255)
                                                                                                              tgatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAg                                        >  2:972643/1‑149 (MQ=255)
                                                                                                               gatAATTGCGGATTGCAGAAATAAAAAGGGCGGTGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGc                                       <  1:833154/149‑1 (MQ=255)
                                                                                                                   aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAACATCATGAAAATTCAATACATTGCAACATTTCCGTAGCCACGATATGCGTAg                                   <  4:1383176/149‑1 (MQ=37)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                  <  4:1477525/149‑1 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                  >  2:47194/1‑149 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                  <  3:1244868/149‑1 (MQ=255)
                                                                                                                        ggATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAt                              >  4:1199312/1‑149 (MQ=255)
                                                                                                                           ttGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAg                           <  1:83374/149‑1 (MQ=255)
                                                                                                                            tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAgg                          <  3:215203/149‑1 (MQ=255)
                                                                                                                              cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTATGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:389370/1‑149 (MQ=255)
                                                                                                                               aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                       >  2:1147079/1‑149 (MQ=255)
                                                                                                                               aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                       >  4:540598/1‑149 (MQ=255)
                                                                                                                                gAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACg                      >  1:534292/1‑149 (MQ=255)
                                                                                                                                 aaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGc                     >  1:713398/1‑149 (MQ=255)
                                                                                                                                 aaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGc                     >  2:294095/1‑149 (MQ=255)
                                                                                                                                       aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTTCCCGGTGCCGCCTGGAACCCCTATCCCGCCTACAAAATCATTGAAATTCAATACATTGCAAGATTTTCGTAGGCCTGGTAAGCGTAGGGCATCCGGCACGCTCGCAt               >  4:767387/1‑149 (MQ=255)
                                                                                                                                       aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgcgg                                                                                                                >  2:361514/1‑52 (MQ=255)
                                                                                                                                       aaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCAt               >  3:10913/1‑149 (MQ=255)
                                                                                                                                          gggCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGt            >  1:160700/1‑149 (MQ=255)
                                                                                                                                           ggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGtt           <  4:411625/149‑1 (MQ=255)
                                                                                                                                                  atgatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCa    <  2:1291502/149‑1 (MQ=255)
                                                                                                                                                    gatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCAtt  >  3:1294245/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATT  >  NC_000913/3877410‑3877706

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