Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,108,462 A→G S295S (AGT→AGC)  SAUSA300_RS10650 ← terminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,108,4620AG100.0% 19.7 / NA 10S295S (AGT→AGC) SAUSA300_RS10650terminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (2/8);  total (2/8)

AAGGTCAACTTCAGGCAAATTCATTCGCTTAGTCATGAATTCGGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTACTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTGGGTTCGCCTTTTCCCACGTCTGTCT  >  NC_007793/2108369‑2108533
                                                                                             |                                                                       
aaGGTCAACTTCAGGCAAATTCATTCGCTTAGTCATGAATTCAGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGc                                                                   <  1:3919271/100‑1 (MQ=255)
     cAACTTCAGGCAAATTCATTCGCTTAGACATGAATTCAGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAgtgt                                                              <  1:4251309/100‑1 (MQ=255)
        cTTCAGGCAAATTCATTCGCTTAGTCATGAATTCAGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGtttt                                                           <  1:66195/100‑1 (MQ=255)
             ggCAAATTCATTCGCTTAGTCATGAATTCAGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGt                                                      <  1:1893402/100‑1 (MQ=255)
                aaaTTCATTCGCTTAGTCATGAATTCAGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTAta                                                   <  1:3020257/100‑2 (MQ=255)
                                        ccaggCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTg                           >  1:1095060/4‑100 (MQ=255)
                                         caggCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTgg                          <  1:3977967/98‑1 (MQ=255)
                                                          ctGAATGGTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTGGGTTCGCCTTTTCCCACg         >  1:847033/2‑100 (MQ=255)
                                                                ggTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTGGGTTCGCCTTTTCCCACgtctgtt   <  1:5775692/100‑2 (MQ=255)
                                                                 gTAAATCGTTATATTCTTCTTCAATCGTGCTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTGGGTTCGCCTTTTCCCACgtctgtgg  <  1:4677190/100‑3 (MQ=255)
                                                                                             |                                                                       
AAGGTCAACTTCAGGCAAATTCATTCGCTTAGTCATGAATTCGGGCTTATTTGAACGGTTGAATGGTAAATCGTTATATTCTTCTTCAATCGTACTTAGCAGTGTTTTAGCGTATTCTGATAACGGTTTATGTAACATTGGGTTCGCCTTTTCCCACGTCTGTCT  >  NC_007793/2108369‑2108533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑