Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 273,644 G→A A266A (GCG→GCA ykfC → CP4‑6 prophage; PF00078 domain‑containing protein YkfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913273,6440GA100.0% 102.3 / NA 29A266A (GCG→GCAykfCCP4‑6 prophage; PF00078 domain‑containing protein YkfC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)

GAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGG  >  NC_000913/273567‑273721
                                                                             |                                                                             
gAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGt                                                                             >  1:63456/1‑80 (MQ=255)
 aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTc                                                                            <  2:63420/80‑1 (MQ=255)
 aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTc                                                                            <  2:63421/80‑1 (MQ=255)
        ccTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGaaaa                                                                    <  1:63423/81‑1 (MQ=255)
        ccTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGaaaa                                                                    <  2:63422/81‑1 (MQ=255)
         cTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAAc                                                                   >  2:63467/1‑81 (MQ=255)
          tGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACt                                                                  >  2:63484/1‑81 (MQ=255)
             gCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGc                                                               <  1:63424/81‑1 (MQ=255)
                   aaGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGa                                                         <  1:63425/81‑1 (MQ=255)
                    aGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGaa                                                        <  1:63426/81‑1 (MQ=255)
                    aGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGa                                                         >  1:63464/1‑80 (MQ=255)
                         gAAAAGATCGGTGGTAATGGAATAACAGTAACCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCcg                                                   <  1:63427/81‑1 (MQ=255)
                           aaaGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCcgcg                                                 >  2:63491/1‑81 (MQ=255)
                                tCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGc                                            <  1:63428/81‑1 (MQ=255)
                                  ggtggtACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCg                                           <  1:63429/80‑1 (MQ=255)
                                     ggtACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACt                                       >  1:63497/1‑81 (MQ=255)
                                      gtACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTg                                      >  1:63496/1‑81 (MQ=255)
                                            gAATAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCt                                 <  1:63430/80‑1 (MQ=255)
                                              aTAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATg                              <  2:63432/81‑1 (MQ=255)
                                              aTAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATg                              <  1:63431/81‑1 (MQ=255)
                                               tAACAGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGc                             <  1:63433/81‑1 (MQ=255)
                                                   aGTATCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCgat                         >  2:63471/1‑81 (MQ=255)
                                                      aTCCAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCgatgat                      >  2:63482/1‑81 (MQ=255)
                                                         cAACGGGGCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGAtttt                   <  2:63434/81‑1 (MQ=255)
                                                              gggCCGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCt              >  2:63492/1‑81 (MQ=255)
                                                                 ccGAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCAt           <  1:63435/81‑1 (MQ=255)
                                                                   gAAGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCg         >  1:63476/1‑81 (MQ=255)
                                                                    aaGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGt        >  2:63474/1‑81 (MQ=255)
                                                                     aGTACGGCAGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGt        >  2:63481/1‑80 (MQ=255)
                                                                          ggcagTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAgg  >  1:63483/1‑81 (MQ=255)
                                                                             |                                                                             
GAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGG  >  NC_000913/273567‑273721

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