Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,655,934 T→C 5.6% intergenic (‑32/+32) yhiD ← / ← hdeB putative Mg(2+) transport ATPase, inner membrane protein/acid‑resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,655,9340TC 5.6% 389.4 / 8.0 125intergenic (‑32/+32)yhiD/hdeBputative Mg(2+) transport ATPase, inner membrane protein/acid‑resistance protein
Reads supporting (aligned to +/- strand):  ref base T (65/53);  new base C (3/4);  total (68/57)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 133 total aligned reads displayed.

GCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGT  >  NC_000913/3655786‑3656082
                                                                                                                                                    |                                                                                                                                                    
gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCt                                                                                                                                                      <  2:1419579/149‑1 (MQ=255)
gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCt                                                                                                                                                      >  2:1098568/1‑149 (MQ=255)
gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCt                                                                                                                                                      <  2:328987/149‑1 (MQ=255)
gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCt                                                                                                                                                      <  2:1147768/149‑1 (MQ=255)
gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCt                                                                                                                                                      <  1:1368200/149‑1 (MQ=255)
  cAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCa                                                                                                                                                    >  1:593531/1‑149 (MQ=255)
        aTACATGTGTGCGTATCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGTCAATTGCCGCGAGTATAAGACGAAAAATAAATTCTGATGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGc                                                                                                                                              <  2:206193/149‑1 (MQ=255)
        aTACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGc                                                                                                                                              <  2:79512/149‑1 (MQ=255)
            aTGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACt                                                                                                                                          >  1:44503/1‑149 (MQ=255)
               tgtgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGTTGAGCCATCTTGCTGCTCCt                                                                                                                                       <  2:618559/149‑1 (MQ=255)
                gtgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGAAATTGCGGCGAGTATAAGTCGTATAATAAATTCTGCGGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAAGTTGCTTCTCCtt                                                                                                                                      >  2:154516/1‑149 (MQ=255)
                 tgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCTGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttt                                                                                                                                     >  1:541565/1‑149 (MQ=255)
                 tgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttt                                                                                                                                     >  2:384012/1‑149 (MQ=255)
                     tAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgc                                                                                                                                 <  1:1426944/149‑3 (MQ=38)
                      aaCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATATCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCGGCTCCtttttgca                                                                                                                                <  2:1400363/149‑4 (MQ=37)
                      aaCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgca                                                                                                                                <  1:189347/149‑4 (MQ=38)
                       aCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACtt                                                                                                                               >  2:43008/1‑149 (MQ=255)
                       aCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACtt                                                                                                                               >  2:417428/1‑149 (MQ=255)
                        cccTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTACTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTg                                                                                                                              >  1:1049400/1‑149 (MQ=255)
                          cTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCa                                                                                                                            <  2:693646/149‑1 (MQ=255)
                           tGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATCGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAc                                                                                                                           <  2:62560/149‑1 (MQ=255)
                                   tGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGAGGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTaa                                                                                                                   <  2:1357470/149‑1 (MQ=255)
                                      cgcgCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTc                                                                                                                >  1:524668/1‑149 (MQ=255)
                                        cgcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCgg                                                                                                              <  2:1364683/149‑1 (MQ=255)
                                          cATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCa                                                                                                            >  1:11367/1‑149 (MQ=255)
                                            tttGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGAAAg                                                                                                          >  1:1119678/1‑149 (MQ=255)
                                             ttGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGt                                                                                                         <  1:872030/149‑1 (MQ=255)
                                               gCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCa                                                                                                       >  2:842263/1‑149 (MQ=255)
                                                 cTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCAtt                                                                                                     >  1:686333/1‑149 (MQ=255)
                                                  tttCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTa                                                                                                    <  1:1280629/149‑1 (MQ=255)
                                                   ttCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAg                                                                                                   <  1:863383/149‑1 (MQ=255)
                                                    tCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGa                                                                                                  >  2:663399/1‑149 (MQ=255)
                                                    tCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGa                                                                                                  <  1:1186210/149‑1 (MQ=255)
                                                     ccATGCCGATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATTGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGAt                                                                                                 >  2:780804/1‑149 (MQ=255)
                                                     ccATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGAt                                                                                                 >  1:784542/1‑149 (MQ=255)
                                                       aTGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGc                                                                                               >  1:290755/1‑149 (MQ=255)
                                                       aTGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATGAATTCGGCAAGTCATTAGATGc                                                                                               >  2:1045368/1‑149 (MQ=255)
                                                               aGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAtttt                                                                                       <  1:1136054/149‑1 (MQ=255)
                                                                gcgcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAttttt                                                                                      <  1:435940/149‑1 (MQ=255)
                                                                gcgcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAttttt                                                                                      <  2:648112/149‑1 (MQ=255)
                                                                  gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGa                                                                                    >  2:1035637/1‑149 (MQ=255)
                                                                   ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                   <  2:1107711/149‑1 (MQ=255)
                                                                   ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                   <  1:50208/149‑1 (MQ=255)
                                                                   ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                   <  2:1006260/149‑1 (MQ=255)
                                                                   ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                   <  2:415345/149‑1 (MQ=255)
                                                                    cacaGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                                                                                  <  1:1035637/149‑1 (MQ=255)
                                                                      caGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGt                                                                                <  2:596482/149‑1 (MQ=255)
                                                                            aTTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGAAATTCATTAGATGCTTGATTTTTGAAGTTATACaa                                                                          >  1:304233/1‑149 (MQ=255)
                                                                                ccGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttt                                                                      >  1:1388290/1‑149 (MQ=255)
                                                                                   cgAGTATAAGACGTTTTATAAATTCTTCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTc                                                                   <  1:750902/149‑1 (MQ=255)
                                                                                   cgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTc                                                                   >  1:1570420/1‑149 (MQ=255)
                                                                                    gAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGTGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCt                                                                  <  1:408269/149‑1 (MQ=255)
                                                                                    gAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCt                                                                  >  2:1264803/1‑149 (MQ=255)
                                                                                    gAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTGGCGACTCCTTTTTACTTGCAGCTCATTAATTCGGGAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttgtggt                                                                  >  2:1488241/1‑144 (MQ=255)
                                                                                        ataAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCgg                                                              >  1:1259138/1‑149 (MQ=255)
                                                                                            gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGttt                                                          >  1:1433261/1‑149 (MQ=255)
                                                                                              cGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGttttt                                                        <  1:1176235/149‑1 (MQ=255)
                                                                                               gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGttttta                                                       >  2:1458313/1‑148 (MQ=255)
                                                                                               gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTc                                                       >  1:150465/1‑149 (MQ=255)
                                                                                                tataATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCt                                                      <  2:1095431/149‑1 (MQ=255)
                                                                                                tataATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCt                                                      >  1:821418/1‑149 (MQ=255)
                                                                                                 ataataAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCtt                                                     <  2:150465/149‑1 (MQ=255)
                                                                                                  taataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTa                                                    >  1:1177353/1‑149 (MQ=255)
                                                                                                   aataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTAc                                                   >  2:483230/1‑149 (MQ=255)
                                                                                                    ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACa                                                  >  2:127050/1‑149 (MQ=255)
                                                                                                    ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACa                                                  >  1:1476630/1‑149 (MQ=255)
                                                                                                    ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACa                                                  >  1:660462/1‑149 (MQ=255)
                                                                                                    ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACa                                                  >  1:1516741/1‑149 (MQ=255)
                                                                                                    ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGAAACTCATAAATACGGCAAGACAATAGATGCATGATTTTTTAAGGTAAACAAATTTTTATTAGGGTTTATCTTACa                                                  >  2:122342/1‑149 (MQ=255)
                                                                                                     taaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAg                                                 >  2:1223147/1‑149 (MQ=255)
                                                                                                        aTTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTAt                                              >  2:32080/1‑149 (MQ=255)
                                                                                                          tCTGCTGTCATATAAAAATGCAAAAAGGATCAGCATGATTGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTAtta                                            <  1:43008/149‑2 (MQ=255)
                                                                                                          tCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTc                                            <  1:168725/149‑1 (MQ=255)
                                                                                                            tgctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGa                                          <  2:1259138/149‑1 (MQ=255)
                                                                                                            tgctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGa                                          <  2:752988/149‑1 (MQ=255)
                                                                                                             gctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGAt                                         <  1:922650/149‑1 (MQ=255)
                                                                                                               tgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCa                                       >  2:1474649/1‑149 (MQ=255)
                                                                                                                 tCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                                     >  2:911917/1‑149 (MQ=255)
                                                                                                                    tataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTAATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTAGATCACTTTa                                  >  2:69715/1‑149 (MQ=255)
                                                                                                                        aaaaTGCAAAAAGGAGTAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGaa                              >  2:1595480/1‑149 (MQ=255)
                                                                                                                          aaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTAGATCACTTTAGGAAtt                            >  1:612447/1‑149 (MQ=255)
                                                                                                                           aagtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGTGGGTTTTTCTTACAGTATTCGATCACTTTAGGAAttt                           >  1:1066777/5‑149 (MQ=37)
                                                                                                                            tGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtg                          >  1:1424789/1‑149 (MQ=255)
                                                                                                                             gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtga                         <  1:911917/149‑1 (MQ=255)
                                                                                                                             gCAAAAACGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGTCAAGTCATTAGACGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtga                         <  2:234345/149‑1 (MQ=255)
                                                                                                                              cAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag                        <  2:892187/149‑1 (MQ=255)
                                                                                                                               aaaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagt                       >  1:592914/1‑149 (MQ=255)
                                                                                                                                aaaaGGAGTAGCAAGTAGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtg                      >  2:777948/1‑149 (MQ=38)
                                                                                                                                     gAGTAGCAAGTTGCGCCAACTTGCTGCTCCATTTTACTAGCACCTCATTAATTCGGCAAGTAATTAGATGCTTGATTATAGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTATCGATGTTGAGTGAGCTc                 >  2:319496/1‑149 (MQ=255)
                                                                                                                                       gcagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCgg               <  2:592914/149‑1 (MQ=255)
                                                                                                                                       gcagcaAGATGGCTCAACTTGCTACTCCTTTTTACTCGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCgg               <  1:71885/149‑1 (MQ=255)
                                                                                                                                         agcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGGTGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt             <  2:24539/149‑1 (MQ=255)
                                                                                                                                         agcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt             <  1:756899/149‑1 (MQ=255)
                                                                                                                                          gcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGttt            >  2:1574443/1‑149 (MQ=255)
                                                                                                                                              gATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCAtt        <  2:1384139/149‑1 (MQ=255)
                                                                                                                                                tGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTa      <  1:1264803/149‑1 (MQ=255)
                                                                                                                                                  gCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa    <  2:496269/149‑1 (MQ=255)
                                                                                                                                                  gCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTCCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa    <  1:842263/149‑1 (MQ=255)
                                                                                                                                                   cTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAg   >  1:1354614/1‑149 (MQ=255)
                                                                                                                                                    tCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGt  >  1:1021353/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGT  >  NC_000913/3655786‑3656082

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