Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,246,085 (C)5→4 coding (210/1071 nt) dhaK ← dihydroxyacetone kinase subunit DhaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,246,0810C.100.0% 38.9 / NA 10coding (214/1071 nt)dhaKdihydroxyacetone kinase subunit DhaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (5/5);  total (5/5)

GTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  NZ_CP009273/1245945‑1246213
                                                                                                                                        |                                                                                                                                    
gTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccc                                                                                                                                    >  1:169776/1‑139 (MQ=255)
                         cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATaa                                                                                                          >  1:53571/1‑139 (MQ=255)
                              gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACcaca                                                                                                     <  1:114798/139‑1 (MQ=255)
                                                       cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccg                                                                            >  2:404247/1‑139 (MQ=255)
                                                            tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAcc                                                                       >  2:430295/1‑139 (MQ=255)
                                                                          ttGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGa                                                         >  1:327617/1‑139 (MQ=255)
                                                                             cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACtt                                                      >  2:263923/1‑139 (MQ=255)
                                                                                                ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                    <  2:169776/138‑1 (MQ=255)
                                                                                                     gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                              <  2:53571/139‑1 (MQ=255)
                                                                                                                 ggtgaAAATTTCGCCCGGACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATaca                  <  1:430295/139‑1 (MQ=255)
                                                                                                                                 ggACAGGCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAg  <  1:263923/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAG  >  NZ_CP009273/1245945‑1246213

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