Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,511,424 T→A 5.1% intergenic (+97/+44) 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,4240TA 5.1% 308.2 / 3.1 95intergenic (+97/+44)ypeC/mntHDUF2502 domain‑containing protein YpeC/Mn(2(+))/Fe(2(+)): H(+) symporter MntH
Reads supporting (aligned to +/- strand):  ref base T (47/43);  new base A (3/2);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCG  >  NC_000913/2511276‑2511573
                                                                                                                                                    |                                                                                                                                                     
ggATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGttt                                                                                                                                                      <  2:683020/149‑1 (MQ=255)
cgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGttt                                                                                                                                                      <  1:765507/150‑1 (MQ=255)
     gTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgcc                                                                                                                                                 <  2:978746/150‑1 (MQ=255)
     gTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATAAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgcc                                                                                                                                                 >  2:1481552/1‑150 (MQ=255)
         cGATCATCGCGGAAAAGGCCCGGGTCATGGGCACCGCCATTAAACGCCTCCGCCAATCCAGCACAATTCCTGAAGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTATTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCAt                                                                                                                                             >  2:1373314/1‑150 (MQ=38)
             catcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAACGCCTTATCCGc                                                                                                                                         <  1:1472501/150‑1 (MQ=39)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAACGCCTTATCCGcc                                                                                                                                        <  2:1480400/150‑1 (MQ=39)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCgtc                                                                                                                                        <  2:1079330/150‑3 (MQ=255)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                                        <  1:1466535/150‑1 (MQ=255)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                                        <  1:133156/150‑1 (MQ=255)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                                        <  1:871924/150‑1 (MQ=255)
              atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                                        <  2:635771/150‑1 (MQ=255)
                cgcgGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCaa                                                                                                                                      >  2:1457372/1‑150 (MQ=255)
                    gAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGt                                                                                                                                  >  1:738834/1‑150 (MQ=255)
                    gAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGt                                                                                                                                  >  2:944640/1‑150 (MQ=255)
                     aaaaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTg                                                                                                                                 <  1:1387014/150‑1 (MQ=255)
                         ggCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAg                                                                                                                             <  1:963123/150‑1 (MQ=255)
                         ggCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAg                                                                                                                             <  2:634121/150‑1 (MQ=255)
                           ccGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGAt                                                                                                                           >  1:556873/1‑150 (MQ=255)
                             ggggTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGc                                                                                                                         <  1:834307/150‑1 (MQ=255)
                             ggggGCATGGGCACCGCCATTAACCGCGTCCGCCACTGGAGCACCATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCGACAACCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGc                                                                                                                         <  1:1524226/150‑1 (MQ=38)
                                gTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGAc                                                                                                                      <  1:1033879/150‑1 (MQ=255)
                                  cATGGGCACCGCCATTTAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTTCCGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCgcagc                                                                                                                    >  2:887648/1‑146 (MQ=255)
                                      ggCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCAt                                                                                                                <  1:946050/150‑1 (MQ=255)
                                       gCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGc                                                                                                                         >  2:316796/1‑140 (MQ=255)
                                       gCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGc                                                                                                                         <  1:316796/140‑1 (MQ=255)
                                       gCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCAtt                                                                                                               <  2:491192/150‑1 (MQ=255)
                                        cACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTc                                                                                                              >  2:205796/1‑150 (MQ=255)
                                         aCCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCa                                                                                                             <  1:1212855/150‑1 (MQ=255)
                                              cATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTAc                                                                                                        >  2:665484/1‑150 (MQ=255)
                                              cATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTAc                                                                                                        <  1:1042650/150‑1 (MQ=255)
                                                 tAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAAt                                                                                                     >  2:581566/1‑150 (MQ=255)
                                                  aaaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATc                                                                                                    <  1:180113/150‑1 (MQ=255)
                                                   aaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATcc                                                                                                   >  2:1170480/1‑150 (MQ=255)
                                                    aCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATccc                                                                                                  <  2:708838/150‑1 (MQ=255)
                                                     cgcgTCCGCCAATTGAGCACAATGCCTGATGCGACGCTAGTGCTTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGGGCCAGATGCGGCGCTCATTCAACTACAATCCCa                                                                                                 <  1:522338/150‑1 (MQ=255)
                                                        ggccGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAgcg                                                                                              <  2:1173081/148‑1 (MQ=255)
                                                        gTCCGCCAAAGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAgcg                                                                                              <  1:587933/150‑1 (MQ=255)
                                                           cGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCg                                                                                           <  2:1198287/150‑1 (MQ=255)
                                                           cGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCg                                                                                           >  1:1293193/1‑150 (MQ=255)
                                                            gCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATTAGGCGTTTACGCCGCATCCGCCAATAGAGCCAGATGCGACGCTCATTCAACGACAATCCCAGCGCCgg                                                                                          >  2:1406356/1‑149 (MQ=255)
                                                             ccAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTc                                                                                         <  1:1481552/150‑1 (MQ=255)
                                                             ccAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTc                                                                                         <  2:693104/150‑1 (MQ=255)
                                                                 tGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCcca                                                                                     >  1:224994/1‑150 (MQ=255)
                                                                   gAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAGGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccacc                                                                                   >  2:6181/1‑150 (MQ=255)
                                                                       acaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaaca                                                                               <  2:1357314/150‑1 (MQ=255)
                                                                       acaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaaca                                                                               <  1:107571/150‑1 (MQ=255)
                                                                       acaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCATAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaaca                                                                               >  2:697003/1‑150 (MQ=255)
                                                                        caATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaacaa                                                                              >  1:910208/1‑150 (MQ=255)
                                                                          aTGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                            <  1:952778/150‑1 (MQ=255)
                                                                            gCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCATAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAg                                                                          <  1:697003/150‑1 (MQ=255)
                                                                                  tGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTc                                                                    >  1:1453627/1‑150 (MQ=255)
                                                                                  tGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTc                                                                    >  2:1033012/1‑150 (MQ=255)
                                                                                   gCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCa                                                                   >  2:825348/1‑150 (MQ=255)
                                                                                    cGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAg                                                                  >  1:316151/1‑150 (MQ=255)
                                                                                     gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgc                                                                 >  1:988700/1‑150 (MQ=255)
                                                                                     gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgc                                                                 >  2:280411/1‑150 (MQ=255)
                                                                                           aGTGCGTCTTATGATGTCTAGATACCGGCTCTTTATCGGGTTGGGGGGTTAAGGCTTTTAGGCCTCTTCCGCCAATAGTGCCAGATGCTACGCTGATTCAAGTACAATCCCAGCGCGGTCCTTACTAACAACCTGATATTTAGCGGCagg                                                           <  2:910208/150‑3 (MQ=255)
                                                                                              gCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                        <  1:944640/150‑1 (MQ=255)
                                                                                               cGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCa                                                       <  1:1033012/150‑1 (MQ=255)
                                                                                                 tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGc                                                     >  1:322683/1‑150 (MQ=255)
                                                                                                  cTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCCGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcc                                                    >  1:1209360/1‑149 (MQ=255)
                                                                                                    tATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcacc                                                  >  1:1222149/1‑150 (MQ=255)
                                                                                                          gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATc                                            >  1:1089393/1‑150 (MQ=255)
                                                                                                           ccTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCa                                           >  1:73116/1‑150 (MQ=255)
                                                                                                            cTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAc                                          <  1:825348/150‑1 (MQ=255)
                                                                                                             tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                         >  1:334878/1‑150 (MQ=255)
                                                                                                               cAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCaacaac                                                                             <  1:956801/112‑1 (MQ=255)
                                                                                                               cAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCaacaac                                                                             >  2:956801/1‑112 (MQ=255)
                                                                                                               cAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCa                                       >  1:108620/1‑150 (MQ=255)
                                                                                                                aaaCCGGCTCTTTAACGCGTAGGGCGGATAAGGCTTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                      <  2:1298012/150‑1 (MQ=255)
                                                                                                                 aaCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGc                                     >  2:102293/1‑150 (MQ=255)
                                                                                                                 aaCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGc                                     >  1:1358515/1‑150 (MQ=255)
                                                                                                                  aCCGGCTCTTTAACGCGTAGGTCGGATAAGGCTTTTACGCCGAATCCGCCAATAGTGCCAGATGCGACGCTCATTCAAATACAATGCGAGCGCCGTCCCCACTAACTACCAGATATTCAGCGCCACGACCAGCACCACAATCACGCAGcc                                    <  2:556873/150‑1 (MQ=255)
                                                                                                                    cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctg                                  <  2:1032568/150‑1 (MQ=255)
                                                                                                                     ggCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCCCCACAATCACCCAGActgt                                 >  2:1167652/1‑150 (MQ=255)
                                                                                                                     ggCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgt                                 >  2:507897/1‑150 (MQ=255)
                                                                                                                     ggCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgt                                 >  2:777678/1‑150 (MQ=255)
                                                                                                                      gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtc                                <  2:988700/150‑1 (MQ=255)
                                                                                                                      gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGAATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtc                                >  2:133323/1‑150 (MQ=255)
                                                                                                                       ctctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtct                               >  2:9175/1‑150 (MQ=255)
                                                                                                                          tttAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGtt                            >  1:811784/1‑150 (MQ=255)
                                                                                                                            tAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGtttt                          <  1:1245536/150‑1 (MQ=255)
                                                                                                                               cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCAGAATCACCCAGCCTGTCTGTTTTAcg                       <  1:1373314/150‑1 (MQ=255)
                                                                                                                               cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcg                       >  1:130116/1‑150 (MQ=255)
                                                                                                                               cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcg                       >  2:1370886/1‑150 (MQ=255)
                                                                                                                               cactattggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcg                       >  1:1193853/8‑150 (MQ=38)
                                                                                                                                gcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgc                      >  1:1416335/1‑150 (MQ=255)
                                                                                                                                gcctAGGGCCGATAAGGCGTTTACGCCGCATCCCCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCACAGCCGTCACCACCAACAACCAGATATTCAGCGCCACCACCAGCACCACAATCAGCCAGCCTGTCTGTTTTAcgc                      >  2:137013/4‑150 (MQ=255)
                                                                                                                                 cgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgcg                     <  2:1246314/150‑1 (MQ=255)
                                                                                                                                  gTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCCGCACCACAATCACCCAGCCTGTCTGTTTTACGCGt                    >  1:1360965/1‑150 (MQ=255)
                                                                                                                                   tAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTGAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                   <  2:1414072/150‑1 (MQ=255)
                                                                                                                                   attggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGCCCAGCACCACCATCACCCAGCCTGTCTGTTTTACGCGtt                   >  2:930575/4‑150 (MQ=255)
                                                                                                                                      ggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc                >  1:885631/1‑150 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCAcc        >  2:732463/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCTACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGa     >  2:16154/1‑150 (MQ=255)
                                                                                                                                                  gTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGAt    <  1:581566/150‑1 (MQ=255)
                                                                                                                                                    ttACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCg  <  2:88866/150‑1 (MQ=255)
                                                                                                                                                    ttACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCg  >  1:142219/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCG  >  NC_000913/2511276‑2511573

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