Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,632,332 2 bp→AG coding (95‑96/357 nt) SAUSA300_RS08070 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,632,3320TA86.7% 30.8 / ‑1.1 15L32F (TTA→TTT) SAUSA300_RS08070hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (4/9);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
*NC_0077931,632,3330AG86.7% 25.2 / ‑3.2 15L32S (TTA→TCA) SAUSA300_RS08070hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (4/9);  total (4/11)
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.00e+00

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCAT  >  NC_007793/1632280‑1632427
                                                    ||                                                                                              
cAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGc                                                  <  1:4721570/100‑1 (MQ=255)
                                tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                  >  1:2390683/7‑100 (MQ=255)
                                tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                  >  1:2946612/7‑100 (MQ=255)
                                 cactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATAc                 >  1:4568314/6‑100 (MQ=255)
                                  actcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACa                <  1:2683281/96‑1 (MQ=255)
                                   ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa               >  1:2131754/4‑100 (MQ=255)
                                   ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa               <  1:5736548/97‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:1147150/100‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:2204886/100‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:2720558/100‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:2887236/100‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:484466/100‑1 (MQ=255)
                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc            <  1:521497/100‑1 (MQ=255)
                                       tACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCa           <  1:4389652/100‑1 (MQ=255)
                                                tttCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCtagt  <  1:2599159/100‑4 (MQ=255)
                                                    ||                                                                                              
CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCAT  >  NC_007793/1632280‑1632427

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