Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 495,541 C→T A226T (GCA→ACA)  SAUSA300_RS02370 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793495,5410CT95.0% 50.4 / ‑3.6 20A226T (GCA→ACA) SAUSA300_RS02370hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (1/18);  minor base G (0/1);  total (1/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

GCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGCCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGCCAGCCTGTAAC  >  NC_007793/495444‑495637
                                                                                                 |                                                                                                
gcgcCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGcat                                                                                                <  1:903470/100‑3 (MQ=255)
 cgcCAAGTGAACCTATTAACATTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCat                                                                                               <  1:900063/100‑1 (MQ=255)
 cgcCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCat                                                                                               <  1:1364485/100‑1 (MQ=255)
 cgcCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCat                                                                                               <  1:129311/100‑1 (MQ=255)
 cgcCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCac                                                                                               <  1:344506/100‑2 (MQ=255)
 cgcCAAGTGAACCTATTAACAATGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCat                                                                                               <  1:141020/100‑1 (MQ=255)
  gcCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCata                                                                                              <  1:5543725/100‑1 (MQ=255)
    cAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGGCacaca                                                                                            <  1:1737247/100‑5 (MQ=255)
     aaGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCatatat                                                                                           <  1:35063/100‑1 (MQ=255)
                           ccTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCta                                                                     >  1:3671808/2‑100 (MQ=255)
                            cTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCtag                                                                    <  1:5078333/100‑2 (MQ=255)
                                                       aTAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACt                                         <  1:752691/100‑1 (MQ=255)
                                                         aaaaaGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTgg                                       <  1:1500094/100‑1 (MQ=255)
                                                             aGTTCATCGTTTCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTAt                                   <  1:5427617/100‑1 (MQ=255)
                                                                       ttCATATTTAATGCCTTTACCATCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAgg                         <  1:2146971/100‑2 (MQ=255)
                                                                                  tgctcttccgaTCTGTCATATATATGATTAATGCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATg              <  1:5169485/90‑1 (MQ=255)
                                                                                  gtgcTTTCCGATCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATg              <  1:5176990/97‑1 (MQ=255)
                                                                                    ctcttccgaTCTGTTATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGcc            <  1:1427141/92‑1 (MQ=255)
                                                                                            aTCTGTCCTATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGCCAGCCTGTa    <  1:1958098/100‑1 (MQ=255)
                                                                                            aTCTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGCCAGCCTGTa    <  1:2235246/100‑1 (MQ=255)
                                                                                              cTGTCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGCCAGCCTGTAAc  <  1:5624611/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTTTACCATCTGCCATATATATGATTAATTCTGTTAAACCTATACTCAAAAAAGTTCCGATAATTGTACTGGCTATGGTAATGAGTGTACGCATATGCCAGCCTGTAAC  >  NC_007793/495444‑495637

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