Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 278,741 A→G V21A (GTA→GCA)  afuC ← 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_000913278,7410AG100.0% 44.8 / NA 14V21A (GTA→GCA) afuCCP4‑6 prophage; ABC transporter ATP‑binding protein AfuC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/10);  total (4/10)

GAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCA  >  NC_000913/278611‑278867
                                                                                                                                  |                                                                                                                              
gAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGcc                                                                                                                              <  1:79636/133‑1 (MQ=255)
      ttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTGGATATGCTCGATTCCCGTATTa                                                                                                                        >  1:250998/1‑133 (MQ=255)
       tGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTAc                                                                                                                       <  2:328984/133‑1 (MQ=255)
             tCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCaa                                                                                                                 >  1:510930/1‑133 (MQ=255)
              cGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCaaa                                                                                                                <  2:34261/133‑1 (MQ=255)
              cGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCaaa                                                                                                                <  2:486617/133‑1 (MQ=255)
                        ttCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTg                                                                                                      <  2:386992/133‑1 (MQ=255)
                                                       gctttGCCGGTGCCGGACGGGCCGTGCACCGTCACCATTTGGCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAg                                                                       <  2:250998/131‑1 (MQ=255)
                                                              cGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGt                                                                <  1:505839/133‑1 (MQ=255)
                                                                 aGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTgcg                                                             <  2:409958/133‑1 (MQ=255)
                                                                                                ccccTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAAc                              >  2:319216/1‑133 (MQ=255)
                                                                                                   cTGCGGGATGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACggg                           <  1:52179/133‑1 (MQ=255)
                                                                                                          aTGGTGAGGTTGATATTGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACg                    >  1:54937/1‑133 (MQ=255)
                                                                                                                          tGTCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGc    <  2:510930/133‑1 (MQ=255)
                                                                                                                            tCGATTGCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCa  >  1:260263/1‑133 (MQ=255)
                                                                                                                                  |                                                                                                                              
GAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCA  >  NC_000913/278611‑278867

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: