Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 964,634 T→C R105R (CGT→CGC ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,6340TC100.0% 594.4 / NA 156R105R (CGT→CGCycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (77/79);  total (77/79)

Only 100 of 159 total aligned reads displayed.

TTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGTATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCATT  >  NC_000913/964487‑964781
                                                                                                                                                   |                                                                                                                                                   
tttCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGcat                                                                                                                                                   >  1:1832318/1‑150 (MQ=255)
   cGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATAtt                                                                                                                                                <  1:901434/150‑1 (MQ=255)
     tttGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTccc                                                                                                                                             >  1:654718/1‑151 (MQ=255)
       tGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCt                                                                                                                                            <  4:567158/150‑1 (MQ=255)
        ggggCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTg                                                                                                                                           <  2:61836/150‑1 (MQ=255)
         gggCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCg                                                                                                                                         >  4:1721982/1‑151 (MQ=255)
           gCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTc                                                                                                                                       >  1:1126614/1‑151 (MQ=255)
             atatTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTc                                                                                                                                       >  4:1561147/1‑149 (MQ=255)
                ttatCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGtc                                                                                                                                  >  4:1466966/1‑151 (MQ=255)
                ttatCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGtc                                                                                                                                  >  4:113619/1‑151 (MQ=255)
                  atCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGtctc                                                                                                                                >  1:965198/1‑151 (MQ=255)
                   tCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCg                                                                                                                               >  3:70042/1‑151 (MQ=255)
                    cAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCg                                                                                                                               >  1:1195302/1‑150 (MQ=255)
                     agcaaagcaaATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTg                                                                                                                             >  4:312949/1‑151 (MQ=255)
                      gcaaagcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGa                                                                                                                            <  1:1620610/151‑1 (MQ=255)
                       caaagcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCGCTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGAt                                                                                                                           <  4:1347716/151‑1 (MQ=255)
                       caaagcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGAt                                                                                                                           <  4:609385/151‑1 (MQ=255)
                          agcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGt                                                                                                                         >  3:254019/1‑150 (MQ=255)
                          agcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTa                                                                                                                        <  1:1437693/151‑1 (MQ=255)
                           gcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTAt                                                                                                                       >  1:459535/1‑151 (MQ=255)
                             aaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATgg                                                                                                                     <  3:18204/151‑1 (MQ=255)
                                 tttGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGa                                                                                                                  >  1:1236934/1‑150 (MQ=255)
                                  ttGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAat                                                                                                                >  4:390154/1‑151 (MQ=255)
                                   tgtgGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAat                                                                                                                >  4:1869981/1‑150 (MQ=255)
                                   tgtgGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAat                                                                                                                <  3:60675/150‑1 (MQ=255)
                                    gtgGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAata                                                                                                               <  3:176024/150‑1 (MQ=255)
                                      gggCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCt                                                                                                            >  2:404553/1‑151 (MQ=255)
                                          aGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCg                                                                                                        <  1:1750977/151‑1 (MQ=255)
                                          aGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCg                                                                                                        >  1:569673/1‑151 (MQ=255)
                                            gagaAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCa                                                                                                      <  1:1347377/151‑1 (MQ=255)
                                              gaAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAg                                                                                                    <  4:1111780/151‑1 (MQ=255)
                                              gaAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAg                                                                                                    <  4:1790374/151‑1 (MQ=255)
                                              gaAACCTTAACTGGCGCGACGCAGGATGCAATTGATGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTGCCGCAAg                                                                                                    <  1:131848/151‑1 (MQ=255)
                                               aaaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATTACTGCAACTGACGGCATGCCCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGCCAGGTCTCGTGATGTATGGCGAATATCTTCAGCAAg                                                                                                    >  4:1785690/1‑150 (MQ=255)
                                               aaaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGc                                                                                                   <  1:939858/151‑1 (MQ=255)
                                               aaaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGc                                                                                                   <  1:1060876/151‑1 (MQ=255)
                                               aaaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGc                                                                                                   <  4:791424/151‑1 (MQ=255)
                                                    ttAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGTGAATATCTTCCGCAAGCGGttt                                                                                              >  3:236608/1‑151 (MQ=255)
                                                     tAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTtg                                                                                             >  4:1267081/1‑151 (MQ=255)
                                                      aaCTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTtgt                                                                                            >  1:62294/1‑151 (MQ=255)
                                                      aaCTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTtgt                                                                                            >  1:1176559/1‑151 (MQ=255)
                                                        cTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGc                                                                                          >  1:1536684/1‑151 (MQ=255)
                                                             gcgACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGac                                                                                     <  1:577177/151‑1 (MQ=255)
                                                               gACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGaca                                                                                    >  3:421804/1‑150 (MQ=255)
                                                                 cGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGacaaca                                                                                 <  4:1237625/151‑1 (MQ=255)
                                                                  gCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGAcaacaa                                                                                <  4:1380035/151‑1 (MQ=255)
                                                                   cAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGAcaacaa                                                                                <  1:1617266/150‑1 (MQ=255)
                                                                   cAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAAt                                                                               >  4:831318/1‑151 (MQ=255)
                                                                   cAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAAt                                                                               <  1:66757/151‑1 (MQ=255)
                                                                    tggATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATg                                                                              >  3:207902/2‑151 (MQ=255)
                                                                    aGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATg                                                                              >  4:341042/1‑151 (MQ=255)
                                                                       aTGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTc                                                                           <  1:503685/151‑1 (MQ=255)
                                                                       aTGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTc                                                                           <  1:290271/151‑1 (MQ=255)
                                                                        tGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTc                                                                           >  1:396666/1‑150 (MQ=255)
                                                                        tGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTc                                                                           >  1:393617/1‑150 (MQ=255)
                                                                         gCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCa                                                                          >  3:314096/1‑150 (MQ=255)
                                                                          cAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAAt                                                                        >  1:819785/1‑151 (MQ=255)
                                                                          cAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAAt                                                                        <  2:372541/151‑1 (MQ=255)
                                                                              tgttgAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGaa                                                                     <  4:1119733/150‑1 (MQ=255)
                                                                              tgttgAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGaa                                                                     <  1:1210347/150‑1 (MQ=255)
                                                                                 tgAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACt                                                                  <  1:113587/150‑1 (MQ=255)
                                                                                  gagaTCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTc                                                                 <  4:1920496/150‑1 (MQ=255)
                                                                                   agaTCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCaa                                                               <  2:480428/151‑1 (MQ=255)
                                                                                       cACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGtt                                                           >  2:439253/1‑151 (MQ=255)
                                                                                       cACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGtt                                                           <  4:502566/151‑1 (MQ=255)
                                                                                       cACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGt                                                            <  4:1588987/150‑1 (MQ=255)
                                                                                        aCTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTc                                                          <  1:1861889/151‑1 (MQ=255)
                                                                                        aCTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTc                                                          <  1:444911/151‑1 (MQ=255)
                                                                                            cAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTg                                                       >  1:55786/1‑150 (MQ=255)
                                                                                                 gACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTc                                                 <  2:252947/151‑1 (MQ=255)
                                                                                                     gCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATgg                                             <  4:140489/151‑1 (MQ=255)
                                                                                                      cATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTcc                                                                                                         <  4:1756717/90‑1 (MQ=255)
                                                                                                      cATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTcc                                                                                                         >  1:1756026/1‑90 (MQ=255)
                                                                                                         gACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCaa                                         <  4:569907/151‑1 (MQ=255)
                                                                                                         gACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCa                                          >  4:1549800/1‑150 (MQ=255)
                                                                                                          aCCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAAc                                        >  4:1885184/1‑151 (MQ=255)
                                                                                                            cACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACtt                                      >  4:1478739/1‑151 (MQ=255)
                                                                                                            cACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACt                                       <  4:1215142/150‑1 (MQ=255)
                                                                                                             aCTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTa                                     <  4:714194/151‑1 (MQ=255)
                                                                                                              cTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTa                                     >  4:1523931/1‑150 (MQ=255)
                                                                                                                cATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTaa                                    <  4:824993/149‑1 (MQ=255)
                                                                                                                cATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatg                                  >  1:864852/1‑151 (MQ=255)
                                                                                                                 aTTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatg                                  >  1:995790/1‑150 (MQ=255)
                                                                                                                  ttACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTaa                                    <  3:10406/147‑1 (MQ=255)
                                                                                                                  ttACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatgat                                >  1:1730912/1‑151 (MQ=255)
                                                                                                                  ttACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatgat                                >  4:914419/1‑151 (MQ=255)
                                                                                                                   tACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatgatg                               >  1:181392/1‑151 (MQ=255)
                                                                                                                          aaaTTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCtt                        <  2:421818/151‑1 (MQ=255)
                                                                                                                              tACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGa                     <  4:1417203/150‑1 (MQ=255)
                                                                                                                                  cATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGcc                <  1:1555310/151‑1 (MQ=255)
                                                                                                                                       aCAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAg              <  4:459704/148‑1 (MQ=255)
                                                                                                                                       aCAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGt           >  1:1016052/1‑151 (MQ=255)
                                                                                                                                        cAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGtt          >  4:1246899/1‑151 (MQ=255)
                                                                                                                                        cAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGtt          >  2:15358/1‑151 (MQ=255)
                                                                                                                                        cAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGtt          <  4:792476/151‑1 (MQ=255)
                                                                                                                                         aaGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTa         <  1:1399489/151‑1 (MQ=255)
                                                                                                                                               gACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCAt   >  1:576425/1‑151 (MQ=255)
                                                                                                                                               gACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCAt   <  4:603098/151‑1 (MQ=255)
                                                                                                                                               gACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCAt   >  1:326289/1‑151 (MQ=255)
                                                                                                                                                aCGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCAtt  >  4:1533521/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGTATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTATGCCATT  >  NC_000913/964487‑964781

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