Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,904,586 T→C V231A (GTC→GCC)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,5860TC100.0% 77.8 / NA 21V231A (GTC→GCC) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/12);  total (9/12)

AAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGT  >  NC_000913/1904453‑1904719
                                                                                                                                     |                                                                                                                                     
aaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                                                  <  2:247832/139‑1 (MQ=255)
aaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                                                  <  2:412004/139‑1 (MQ=255)
          aaGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                                        <  1:588850/139‑1 (MQ=255)
          aaGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                                        <  2:173582/139‑1 (MQ=255)
                 gTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTa                                                                                                                 <  2:163161/139‑1 (MQ=255)
                          cGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGc                                                                                                        >  2:608581/1‑139 (MQ=255)
                                 tttGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                       <  1:49293/117‑1 (MQ=255)
                                 tttGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                       >  2:49293/1‑117 (MQ=255)
                                           gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTc                                                                                                                            >  1:508486/1‑102 (MQ=255)
                                           gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTc                                                                                                                            <  2:508486/102‑1 (MQ=255)
                                                    ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                                              <  1:200619/139‑1 (MQ=255)
                                                    ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                                              <  2:460362/139‑1 (MQ=255)
                                                       aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                                           <  1:34444/139‑1 (MQ=255)
                                                                             cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaaa                                                     >  2:441838/1‑139 (MQ=255)
                                                                               cgcgGGTTGTCTCTCGCATTACTGACCAGACGGGCAAGGAACACGTTACTACTGCCCAGACTATTCTGGGCAAGGTAATGGCAGGCCTGGTACCACTGCTGCTGACCTTTGCTGGTATGTGGCTACTGCGCAAAAAAGt                                                   >  1:271539/1‑139 (MQ=255)
                                                                               cGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGt                                                   >  2:170673/1‑139 (MQ=255)
                                                                                   ggTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                                               <  1:84515/139‑1 (MQ=255)
                                                                                    gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                                                  >  1:7922/1‑55 (MQ=255)
                                                                                    gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                                                  <  2:7922/55‑1 (MQ=255)
                                                                                            tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTgg                                      >  2:359821/1‑139 (MQ=255)
                                                                                                                                tactGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGt  >  2:117173/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
AAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGT  >  NC_000913/1904453‑1904719

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