Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,632,339 2 bp→AA coding (88‑89/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,3390GA81.2% 24.2 / ‑0.3 16T30I (ACA→ATA) SAUSA300_RS08070hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (4/9);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0077931,632,3400TA87.5% 35.9 / ‑1.3 16T30S (ACA→TCA) SAUSA300_RS08070hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (4/10);  total (4/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 = 9.82e-01

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAA  >  NC_007793/1632280‑1632432
                                                           ||                                                                                            
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)
                                                     gaatttGAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaa  <  1:3280461/97‑1 (MQ=255)
                                                           ||                                                                                            
CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAA  >  NC_007793/1632280‑1632432

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