Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,298,035 T→C E35G (GAA→GGA)  yraQ ← putative inner membrane permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,298,0350TC100.0% 144.6 / NA 41E35G (GAA→GGA) yraQputative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/20);  total (21/20)

CGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTG  >  NC_000913/3297911‑3298169
                                                                                                                            |                                                                                                                                      
cgcgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATAc                                                                                                                          >  1:1621510/1‑139 (MQ=255)
cgcgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATAc                                                                                                                          <  1:76495/139‑1 (MQ=255)
  cgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCa                                                                                                                        <  2:709402/139‑1 (MQ=255)
   gAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCag                                                                                                                       >  2:1317659/1‑139 (MQ=255)
     ggAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCagag                                                                                                                     >  1:524229/1‑139 (MQ=255)
          tAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                    >  1:1708294/1‑135 (MQ=255)
          tAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                    <  2:1708294/135‑1 (MQ=255)
          tAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCgg                                                                                                                >  2:596474/1‑139 (MQ=255)
                   gcgTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                       <  1:1317659/139‑1 (MQ=255)
                   gcgTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                       <  1:596474/139‑1 (MQ=255)
                              cgccgcCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTca                                                                                                                               >  1:1054836/1‑104 (MQ=255)
                              cgccgcCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTca                                                                                                                               <  2:1054836/104‑1 (MQ=255)
                                   ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    <  2:784852/94‑1 (MQ=255)
                                   ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    >  1:784852/1‑94 (MQ=255)
                                            ggTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                       >  2:89691/1‑114 (MQ=255)
                                            ggTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                       <  1:89691/114‑1 (MQ=255)
                                                aGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCcaccac                                                                          >  1:60488/1‑139 (MQ=255)
                                                aGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCcaccac                                                                          >  2:1280078/1‑139 (MQ=255)
                                                            cgcAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACg                                                                                                     >  2:223299/1‑100 (MQ=255)
                                                            cgcAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACg                                                                                                     <  1:223299/100‑1 (MQ=255)
                                                             gcAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTg                                                             <  2:340771/139‑1 (MQ=255)
                                                              cAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTgc                                                            <  2:60488/139‑1 (MQ=255)
                                                              cAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTgc                                                            <  2:113214/139‑1 (MQ=255)
                                                                   aTAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCt                                                       <  2:756438/139‑1 (MQ=255)
                                                                            ccGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAg                                              >  1:1502896/1‑139 (MQ=255)
                                                                                 aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                                  >  1:709465/1‑50 (MQ=255)
                                                                                 aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                                  <  2:709465/50‑1 (MQ=255)
                                                                                 aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                         <  2:524229/139‑1 (MQ=255)
                                                                                  cTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                        <  1:1706800/139‑1 (MQ=255)
                                                                                  cTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                        <  1:1280078/139‑1 (MQ=255)
                                                                                     tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                             <  2:659136/131‑1 (MQ=255)
                                                                                     tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                             >  1:659136/1‑131 (MQ=255)
                                                                                          ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                    >  1:1503938/1‑55 (MQ=255)
                                                                                          ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                    <  2:1503938/55‑1 (MQ=255)
                                                                                          ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATa                                >  2:647434/1‑139 (MQ=255)
                                                                                          ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACG‑GAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAg                               >  2:425144/1‑139 (MQ=255)
                                                                                                       aaCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCgg                                                                                     <  1:344955/73‑1 (MQ=255)
                                                                                                       aaCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCgg                                                                                     >  2:344955/1‑73 (MQ=255)
                                                                                                             ttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGt             >  2:1229936/1‑139 (MQ=255)
                                                                                                               gCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTa           >  1:232873/1‑139 (MQ=255)
                                                                                                                        aGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTg  >  2:1328843/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                      
CGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTG  >  NC_000913/3297911‑3298169

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