Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,579,617 G→A D70D (GAC→GAT)  SAUSA300_RS07730 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,579,6170GA92.9% 28.5 / ‑4.5 14D70D (GAC→GAT) SAUSA300_RS07730hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (2/11);  total (2/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

TCACTATTTTTAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAATCCAAACGCCAACATGTCCGCTAATTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAA  >  NC_007793/1579525‑1579709
                                                                                            |                                                                                            
tctgtATCTTTAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGcca                                                                                       <  1:3202019/96‑3 (MQ=255)
tctgtATCTTTAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGcca                                                                                       <  1:4716030/96‑3 (MQ=255)
   gtATCTTTAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATccaccaact                                                                                    >  1:4207578/2‑93 (MQ=255)
        tttAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCt                                                                               <  1:4850127/100‑1 (MQ=255)
                            cttcttTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCtt                                                           <  1:4410432/100‑1 (MQ=255)
                            cttcttTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCtt                                                           <  1:4985131/100‑1 (MQ=255)
                             agacgtgtgctctTCCGATCTCACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCttg                                                          <  1:2953700/91‑2 (MQ=255)
                                  cgaTCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCtag                                                     <  1:3168907/98‑3 (MQ=255)
                                   gaTCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAATCCAAACGCCAACATGTCAGCTAATTCATCAAGTTGTACGTCTAACGGCTTACCTGGt                                                    <  1:5005589/99‑1 (MQ=255)
                                      ttccgatcttgaCACTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCTGGTttc                                                 <  1:4686840/91‑1 (MQ=255)
                                                   acTCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCt                                    <  1:2263180/100‑1 (MQ=255)
                                                     tCCTACTTGATTCGCAATACTCAAACCAAACGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCTTa                                  <  1:2102296/100‑1 (MQ=255)
                                                                                 aaCGCCAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCaa      <  1:674031/100‑1 (MQ=255)
                                                                                     gcAACATATCCGCCAACTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCTTAAAAGTTTCCAATGTATTAAACCATTCAAAGaa  >  1:5587756/2‑100 (MQ=255)
                                                                                            |                                                                                            
TCACTATTTTTAAAACTTGATTCAATCGCTTCTTTTATCTCTTCTGATGACACTCCTACTTGATTCGCAATACTCAATCCAAACGCCAACATGTCCGCTAATTCATCTAACTGAACATCTAACGGCTTACCTGGTTTCTTCTTCCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAA  >  NC_007793/1579525‑1579709

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: ‑