Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,640,202 G→A A55V (GCA→GTA)  ydfR ← Qin prophage; protein YdfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,640,2020GA100.0% 59.7 / NA 18A55V (GCA→GTA) ydfRQin prophage; protein YdfR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/7);  total (11/7)

CAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGC  >  NC_000913/1640126‑1640277
                                                                            |                                                                           
cAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTAcc                                                                           <  2:285075/79‑1 (MQ=255)
       gACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAAc                                                                  <  2:285076/81‑1 (MQ=255)
          tCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTc                                                                >  1:285097/1‑80 (MQ=255)
            aaCTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAg                                                             >  1:285104/1‑81 (MQ=255)
             aCTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAgc                                                            >  2:285103/1‑81 (MQ=255)
                ttttCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAgcg                                                           >  2:285096/1‑79 (MQ=255)
                  ttCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAg                                                        <  1:285077/80‑1 (MQ=255)
                      gACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCt                                                   >  2:285101/1‑81 (MQ=255)
                          aTTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCaa                                               >  1:285094/1‑81 (MQ=255)
                                cTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGt                                         >  1:285093/1‑81 (MQ=255)
                                  gCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTc                                        >  2:285098/1‑80 (MQ=255)
                                      cTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATg                                    <  2:285078/80‑1 (MQ=255)
                                         ccAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCg                                 <  2:285079/80‑1 (MQ=255)
                                         ccAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGa                                >  1:285109/1‑81 (MQ=255)
                                                gACTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAg                         <  1:285080/81‑1 (MQ=255)
                                                 aCTTGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAgg                        <  2:285081/81‑1 (MQ=255)
                                                   ttGCCGGACGCTCCTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGt                       >  2:285106/1‑80 (MQ=255)
                                                               ccTTCTGGCTGCTACCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCAcc          >  2:285105/1‑81 (MQ=255)
                                                                       ctgctaCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGc  <  2:285082/81‑1 (MQ=255)
                                                                            |                                                                           
CAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGC  >  NC_000913/1640126‑1640277

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