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% 161.3 / NA 44V231A (GTC→GCC) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/27);  total (17/27)

AAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCT  >  NC_000913/1904453‑1904699
                                                                                                                                     |                                                                                                                 
aaaCTGGCGGAAGGGGCGTCTATCCTCGGCCTGTTGGTCATGGGGGCATTGGTTAACAAGTGGACACATGTTAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                              <  1:741430/139‑1 (MQ=255)
aaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                              <  1:505823/139‑1 (MQ=255)
 aaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGAc                                                                                                             <  2:9173/139‑1 (MQ=255)
      aCGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTc                                                                                                        <  1:862598/139‑1 (MQ=255)
       cGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCt                                                                                                       >  2:1263150/1‑139 (MQ=255)
          aaGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                    <  1:967387/139‑1 (MQ=255)
          aaGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                    <  2:803851/139‑1 (MQ=255)
                           ggCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                   <  1:1224035/139‑1 (MQ=255)
                            gCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAgg                                                                                     >  1:188865/1‑136 (MQ=255)
                            gCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAgg                                                                                     <  2:188865/136‑1 (MQ=255)
                            gCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCt                                                                                  <  1:1263150/139‑1 (MQ=255)
                                gTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                              >  2:485598/1‑139 (MQ=255)
                                           gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgct                                                                   >  2:601817/1‑139 (MQ=255)
                                            ggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgctg                                                                  <  1:485598/139‑1 (MQ=255)
                                                    ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                          <  1:298632/139‑1 (MQ=255)
                                                       aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                       <  2:574949/139‑1 (MQ=255)
                                                       aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                       <  1:571937/139‑1 (MQ=255)
                                                       aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                       >  2:1074409/1‑139 (MQ=255)
                                                              gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAgg                                                                                     <  2:1146719/102‑1 (MQ=255)
                                                              gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAgg                                                                                     >  1:1146719/1‑102 (MQ=255)
                                                                   aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACt                                           >  2:667708/1‑139 (MQ=255)
                                                                        aaCATCCCGTTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                              >  2:884071/1‑67 (MQ=255)
                                                                        aaCATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                              <  1:884071/67‑1 (MQ=255)
                                                                        aaCATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCa                                      >  1:361548/1‑139 (MQ=255)
                                                                          cATCCCGCTGGTTGTCTCTCGCACTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCCGACTATTCTGGACCAGTGAATGCCAGGCCTGGTACGACTGCTGCTGACCTTTGCTTGTATGTGGCCAGTGCGCaaa                                    <  1:859521/139‑1 (MQ=255)
                                                                                   ggTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                           <  1:667708/139‑1 (MQ=255)
                                                                                   ggTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACAACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                           <  1:1044434/139‑1 (MQ=255)
                                                                                    gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                          <  2:1178001/139‑1 (MQ=255)
                                                                                      tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCg                        <  1:1027486/139‑1 (MQ=255)
                                                                                      tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCg                        <  2:1234339/139‑1 (MQ=255)
                                                                                              gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                        >  2:261659/1‑67 (MQ=255)
                                                                                              gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                        <  1:261659/67‑1 (MQ=255)
                                                                                                      gACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                         <  1:721754/90‑1 (MQ=255)
                                                                                                      gACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                         >  2:721754/1‑90 (MQ=255)
                                                                                                       aCCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTgg       >  2:426939/1‑139 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                    <  2:1005623/41‑1 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGa                                                                                                    >  1:1005623/1‑41 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                   >  2:483678/1‑42 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                   <  1:483678/42‑1 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                   <  1:499627/58‑1 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                   >  2:499627/1‑58 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTaa                            >  1:84111/1‑113 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTaa                            <  2:84111/113‑1 (MQ=255)
                                                                                                            aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGcttct  >  2:359392/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                 
AAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCT  >  NC_000913/1904453‑1904699

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