Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,511,417 A→G 5.1% intergenic (+90/+51) ypeC → / ← mntH DUF2502 domain‑containing protein YpeC/Mn(2(+))/Fe(2(+)): H(+) symporter MntH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,511,4170AG 5.1% 364.3 / 6.3 117intergenic (+90/+51)ypeC/mntHDUF2502 domain‑containing protein YpeC/Mn(2(+))/Fe(2(+)): H(+) symporter MntH
Reads supporting (aligned to +/- strand):  ref base A (51/59);  new base G (2/4);  total (53/64)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 118 total aligned reads displayed.

GGCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCA  >  NC_000913/2511268‑2511565
                                                                                                                                                     |                                                                                                                                                    
ggCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATa                                                                                                                                                      >  2:861222/1‑150 (MQ=255)
   tGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                                   <  1:760654/150‑1 (MQ=255)
   tGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                                   <  1:1173622/150‑1 (MQ=255)
   cgcgatgATCGTCACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                                   >  1:383801/7‑150 (MQ=255)
    ggCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGc                                                                                                                                                  <  1:1026850/150‑1 (MQ=255)
    ggCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGc                                                                                                                                                  >  2:1050995/1‑150 (MQ=255)
       gcgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGtt                                                                                                                                               <  2:1109819/150‑1 (MQ=255)
       gcgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGtt                                                                                                                                               <  1:414948/150‑1 (MQ=255)
       gcgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGtt                                                                                                                                               <  2:270483/150‑1 (MQ=255)
        cgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGttt                                                                                                                                              <  2:40540/150‑1 (MQ=255)
         gATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTa                                                                                                                                             <  2:632033/150‑1 (MQ=255)
            cGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgc                                                                                                                                          <  2:1242003/150‑1 (MQ=255)
             gTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgcc                                                                                                                                         >  2:1300701/1‑150 (MQ=255)
              tGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgccg                                                                                                                                        <  2:1313943/150‑1 (MQ=255)
              tGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgccg                                                                                                                                        >  2:533013/1‑150 (MQ=255)
              tGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgccg                                                                                                                                        >  1:236531/1‑150 (MQ=255)
               gACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgccgc                                                                                                                                       <  2:970633/150‑1 (MQ=255)
                  gATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAAcgccttatc                                                                                                                                    <  1:1327904/150‑6 (MQ=39)
                  gATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGGGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAAcgccttatc                                                                                                                                    <  1:1333599/150‑6 (MQ=38)
                      atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAACGCCTTATCCGcc                                                                                                                                <  1:247115/150‑1 (MQ=39)
                        cgcgGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAACGCCTTATCCGccct                                                                                                                              <  1:1335342/150‑3 (MQ=38)
                         gcgGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAAt                                                                                                                             >  2:459101/1‑150 (MQ=255)
                           ggAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAg                                                                                                                                                    <  2:454275/125‑1 (MQ=255)
                           ggAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAg                                                                                                                                                    >  1:454275/1‑125 (MQ=255)
                             aaaaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTg                                                                                                                         <  1:454767/150‑1 (MQ=255)
                             aaaaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAAAGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTg                                                                                                                         >  1:1281474/1‑150 (MQ=255)
                              aaaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGc                                                                                                                        <  2:859818/150‑1 (MQ=255)
                                 ggCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAg                                                                                                                     >  2:567577/1‑150 (MQ=255)
                                 ggCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAg                                                                                                                     <  2:610067/150‑1 (MQ=255)
                                    cGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATg                                                                                                                  <  1:1633328/150‑1 (MQ=255)
                                    cGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATg                                                                                                                  <  2:639264/150‑1 (MQ=255)
                                      gggTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCg                                                                                                                <  2:1404777/150‑1 (MQ=255)
                                       ggTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGGGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGa                                                                                                               <  1:1077838/150‑1 (MQ=255)
                                          cATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGc                                                                                                            >  1:913545/1‑150 (MQ=255)
                                          cATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGc                                                                                                            <  1:653836/150‑1 (MQ=255)
                                                  ccGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCaa                                                                                                    >  2:172224/1‑150 (MQ=255)
                                                   cGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAAc                                                                                                   <  1:167645/150‑1 (MQ=255)
                                                   cGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAAc                                                                                                   <  1:459101/150‑1 (MQ=255)
                                                    gCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACt                                                                                                  <  2:481636/150‑1 (MQ=255)
                                                    gCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACt                                                                                                  <  1:144601/150‑1 (MQ=255)
                                                      cATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTAc                                                                                                <  2:587768/150‑1 (MQ=255)
                                                           aaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATcc                                                                                           >  2:159808/1‑150 (MQ=255)
                                                           aaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATcc                                                                                           >  1:1047759/1‑150 (MQ=255)
                                                                gTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAgcg                                                                                      <  1:1349595/150‑1 (MQ=255)
                                                                      cAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTcc                                                                                <  2:236531/150‑1 (MQ=255)
                                                                      cAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTcc                                                                                >  2:904285/1‑150 (MQ=255)
                                                                      cAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGGTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTcc                                                                                <  2:485543/150‑1 (MQ=255)
                                                                          ggAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccac                                                                            <  2:1047759/150‑1 (MQ=255)
                                                                           gAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccacc                                                                           <  1:171462/150‑1 (MQ=255)
                                                                             gCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaa                                                                         <  1:632038/150‑1 (MQ=255)
                                                                             gCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACaaa                                                                         >  2:230720/1‑150 (MQ=255)
                                                                              cacaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTGACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaac                                                                        <  2:1613482/150‑1 (MQ=255)
                                                                              cacaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaac                                                                        <  1:282106/150‑1 (MQ=255)
                                                                                   tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCa                                                                   <  1:533013/150‑1 (MQ=255)
                                                                                     ccTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGa                                                                 <  2:1515070/150‑1 (MQ=255)
                                                                                       tGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGata                                                               <  1:563612/150‑1 (MQ=255)
                                                                                         aTGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATAtt                                                             >  2:1368560/1‑150 (MQ=255)
                                                                                          tGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTc                                                            >  1:997479/1‑150 (MQ=255)
                                                                                          tGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTc                                                            >  1:149721/1‑150 (MQ=255)
                                                                                           gCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCa                                                           >  2:1392744/1‑150 (MQ=255)
                                                                                           gCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCa                                                           <  1:652732/150‑1 (MQ=255)
                                                                                            cGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAg                                                          <  1:861222/150‑1 (MQ=255)
                                                                                            cGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAg                                                          >  1:1070547/1‑150 (MQ=255)
                                                                                              aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                        >  1:1438320/1‑150 (MQ=255)
                                                                                              aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCAAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                        <  1:1050995/150‑1 (MQ=255)
                                                                                               cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTATTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAAAAACCAGATATTAAgcgc                                                       >  2:520663/1‑150 (MQ=255)
                                                                                                            tATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcacc                                          <  2:1301773/150‑1 (MQ=255)
                                                                                                               cAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                       >  2:1108089/1‑150 (MQ=255)
                                                                                                                aGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACaa                                      <  1:239443/150‑1 (MQ=255)
                                                                                                                   ccTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCa                                   >  1:807378/1‑150 (MQ=255)
                                                                                                                   ccTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCa                                   >  2:1002146/1‑150 (MQ=255)
                                                                                                                     tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                 >  1:546182/1‑150 (MQ=255)
                                                                                                                     tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                 <  2:173782/150‑1 (MQ=255)
                                                                                                                     tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                 >  2:1385683/1‑150 (MQ=255)
                                                                                                                     tACAAACCGGCTCTTTAACGAGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCACATATTCAGCGCCACGACCAGCACCACAATCAcc                                 >  1:1353885/1‑150 (MQ=255)
                                                                                                                      aCAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAccc                                >  2:1064938/1‑150 (MQ=255)
                                                                                                                          aCCGGCGCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACATCCAGATATTCTGCGCCACGACCAGCACGACAATCTCCCAGcc                            <  2:738000/150‑1 (MQ=39)
                                                                                                                            cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctg                          >  2:925311/1‑150 (MQ=255)
                                                                                                                              gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtc                        <  2:1281474/150‑1 (MQ=255)
                                                                                                                                tctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctg                      <  2:383801/150‑1 (MQ=255)
                                                                                                                                 ctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                     <  1:1314295/150‑1 (MQ=255)
                                                                                                                                   ttAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCAGTCAACTACAATCCCAGCGCCGTCCCCCCCCACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCAAGActgtctgtgt                   >  2:153225/1‑148 (MQ=39)
                                                                                                                                    tAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGtttt                  <  1:219449/150‑1 (MQ=255)
                                                                                                                                       cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcg               <  2:1410918/150‑1 (MQ=255)
                                                                                                                                       cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcg               >  2:597449/1‑150 (MQ=255)
                                                                                                                                        actattggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCCCCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgc              >  1:836092/7‑150 (MQ=38)
                                                                                                                                         cgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgcg             >  2:985003/1‑150 (MQ=255)
                                                                                                                                         cgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgcg             >  2:979085/1‑150 (MQ=255)
                                                                                                                                           tAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt           >  1:156020/1‑150 (MQ=255)
                                                                                                                                            aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGttt          >  1:870284/1‑150 (MQ=255)
                                                                                                                                             gggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg         >  2:145548/1‑150 (MQ=255)
                                                                                                                                              ggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc        <  1:925311/150‑1 (MQ=255)
                                                                                                                                              ggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTAAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc        >  1:460964/1‑150 (MQ=255)
                                                                                                                                                cGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTg      <  2:1516671/150‑1 (MQ=255)
                                                                                                                                                 ggATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGt     >  1:1469229/1‑150 (MQ=255)
                                                                                                                                                  gATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGtt    >  2:1185611/1‑150 (MQ=255)
                                                                                                                                                  gATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGtt    <  1:1002146/150‑1 (MQ=255)
                                                                                                                                                   aTAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTc   >  2:1086207/1‑150 (MQ=255)
                                                                                                                                                    tAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACTACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCa  <  1:567577/150‑1 (MQ=255)
                                                                                                                                                    tAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCa  <  1:1392744/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCA  >  NC_000913/2511268‑2511565

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