Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,601,699 T→C intergenic (‑29/+234) SAUSA300_RS13370 ← / ← SAUSA300_RS13375 hypothetical protein/short‑chain dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,601,6990TC100.0% 43.6 / NA 17intergenic (‑29/+234)SAUSA300_RS13370/SAUSA300_RS13375hypothetical protein/short‑chain dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/15);  total (2/15)

CCAATTTTATCTGCAATCTCTCTTATAGCTTGATATAGCAATCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGATAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACAAT  >  NC_007793/2601603‑2601775
                                                                                                |                                                                            
ccAATTTTATCTGCAATCTCTCTTATAGCTTGATATAGCAATCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACaaa                                                                           >  1:5561636/1‑100 (MQ=255)
          cTGCAATCTCTCTTATAGCTTGATATAGCAATCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTa                                                                 <  1:1141204/100‑1 (MQ=255)
                tctctcTTGTAGCTTGATATAGCAATCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTaa                                                           <  1:437677/100‑1 (MQ=255)
                                        aTCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAGCTCg                                   <  1:1429084/100‑1 (MQ=255)
                                        aTCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCg                                   <  1:1492281/100‑1 (MQ=255)
                                        aTCTAAAACTTTCTTTTGCTTTTAACATAAAATTCTCCTTAAGTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCg                                   <  1:4379644/100‑1 (MQ=255)
                                            aaaaCTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGAt                               >  1:4752769/1‑100 (MQ=255)
                                                                aCATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACg           <  1:2065064/100‑1 (MQ=255)
                                                                aCATAAAATTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACg           <  1:5339268/100‑1 (MQ=255)
                                                                      cgatctcCTTAGTTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTAc     <  1:3418299/97‑1 (MQ=255)
                                                                      aaTTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTAc     <  1:114839/100‑1 (MQ=255)
                                                                      aaTTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTAc     <  1:2069951/100‑1 (MQ=255)
                                                                       aTTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACa    <  1:3583747/100‑1 (MQ=255)
                                                                       aTTCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACa    <  1:1816886/100‑1 (MQ=255)
                                                                        ttCTCCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACac   <  1:5147420/100‑2 (MQ=255)
                                                                         tctcCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACAAt  <  1:2843780/100‑1 (MQ=255)
                                                                         tctcCTTAATTTATTAGAATAGACAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACAAt  <  1:2781615/100‑1 (MQ=255)
                                                                                                |                                                                            
CCAATTTTATCTGCAATCTCTCTTATAGCTTGATATAGCAATCTAAAACTTTCTTTTTCTTTTAACATAAAATTCTCCTTAATTTATTAGAATAGATAAAATACAGCTTATTTTAAAGATACTCCAAAAGTTTGAACTCGAGATTCTGATAAATTAGATGCACGACCTACAAT  >  NC_007793/2601603‑2601775

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