Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 277,660 T→G intergenic (+13/+96) mmuM → / ← afuC homocysteine S‑methyltransferase/CP4‑6 prophage; ABC transporter ATP‑binding protein AfuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913277,6600TG100.0% 33.6 / NA 11intergenic (+13/+96)mmuM/afuChomocysteine S‑methyltransferase/CP4‑6 prophage; ABC transporter ATP‑binding protein AfuC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (4/7);  total (4/7)

ATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCG  >  NC_000913/277532‑277786
                                                                                                                                |                                                                                                                              
aTTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTAc                                                                                                                            <  2:384294/133‑1 (MQ=255)
 ttGCGCGCAGCTGGCCGATTATCTTCCGAAGTGGCAGGCCGATGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACa                                                                                                                           <  2:166047/133‑1 (MQ=255)
                ggATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGc                                                                                                            <  1:89108/133‑1 (MQ=255)
                                 ggCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAgcgcagcgca                                                                                           <  2:473587/133‑1 (MQ=255)
                                                                     gTCGCGCCACGCCTGAGGATATTGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATc                                                       <  1:250054/133‑1 (MQ=255)
                                                                       cGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCg                                                     <  1:380066/133‑1 (MQ=255)
                                                                                   gCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAAcc                                         >  2:408398/1‑133 (MQ=255)
                                                                                    cGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACca                                        >  1:480802/1‑133 (MQ=255)
                                                                                                             aaGCTGAGGGTTTATCGGGGCTTCATCGTTTATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAg               <  1:182766/133‑1 (MQ=255)
                                                                                                             aaGCTGAGGGTTTATCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAg               >  1:325569/1‑133 (MQ=255)
                                                                                                                          aTCGGGGCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCg  >  1:111054/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                              
ATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCGCCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCG  >  NC_000913/277532‑277786

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