Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,586,045 2 bp→CG coding (62‑63/303 nt) SAUSA300_RS07790 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,586,0450GC100.0% 30.8 / NA 14S21R (AGC→AGG) SAUSA300_RS07790hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (1/13);  total (1/13)
*NC_0077931,586,0460CG100.0% 30.7 / NA 14S21T (AGC→ACC) SAUSA300_RS07790hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (1/13);  total (1/13)

AATCTCCGTCGATAACTCCTATTGATGAGAAAAGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTCCGCTTTTTTCAGATAATAAGTCTTTTACTTCTTCCTGAGTCATGTCTTTATATTTTTTAGTCATTGTTGACTTCCTCCTCGTTT  >  NC_007793/1585954‑1586127
                                                                                           ||                                                                                 
aatCTCCGTCGATAACTCCTATTGATGAGAAAAGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTc                                                                            <  1:1968931/100‑1 (MQ=255)
aatCTCCGTCGATAACTCCTATTGATGAGAAAAGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTc                                                                            <  1:4658337/100‑1 (MQ=255)
                             aaaaGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTCCAGATAATAAGTCTTTTATTTCTTCCTgct                                               <  1:4812023/100‑3 (MQ=255)
                              aaaGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGt                                              <  1:5413007/100‑1 (MQ=255)
                                aGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGtgt                                            <  1:2897194/100‑3 (MQ=255)
                                  cAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATg                                          <  1:1888957/100‑1 (MQ=255)
                                   aaaaTATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGt                                         <  1:1343670/100‑1 (MQ=255)
                                   aaaaTATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGt                                         <  1:4747142/100‑1 (MQ=255)
                                      atatCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCtt                                      <  1:3536375/100‑1 (MQ=255)
                                       tatCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCttt                                     <  1:1153941/100‑1 (MQ=255)
                                                        cccTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCTTTATAATTTTTAGTCAtgg                    >  1:1552225/1‑98 (MQ=25)
                                                         ccTTAATTTCTTTCGCTAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCTTTATAATTTTTAGTCAtgga                   <  1:2831839/100‑4 (MQ=25)
                                                                         tAATTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCTTTATAATTTTTAGTCATGGTTGACTtcctccgagag   <  1:2200744/100‑6 (MQ=18)
                                                                          aaTTCATACAATTCTGCCGTTTTTTCAGATAATAAGTCTTTTATTTCTTCCTGCGTCATGTCTTTATAATTTTTAGTCATGGTTGACTtcctccgttttt  <  1:447537/100‑7 (MQ=18)
                                                                                           ||                                                                                 
AATCTCCGTCGATAACTCCTATTGATGAGAAAAGCAAAATATCAAATTTACTTTCTCCCTTAATTTCTTTCGCTAATTCATACAATTCTCCGCTTTTTTCAGATAATAAGTCTTTTACTTCTTCCTGAGTCATGTCTTTATATTTTTTAGTCATTGTTGACTTCCTCCTCGTTT  >  NC_007793/1585954‑1586127

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