Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 641,667 Δ1 bp coding (367/453 nt) SAUSA300_RS03035 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793641,6660C.92.3% 39.9 / ‑1.8 13coding (368/453 nt)SAUSA300_RS03035membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (10/2);  total (11/2)
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.50e-01

AACAAATATTTATAATCGTTTAAAAGATTTCACGAACCCAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCGCCTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTGAAATTTAATAAGACTAAATTTTTTGTAGCATCGTATGCTAAGCCACCA  >  NC_007793/641572‑641758
                                                                                              |                                                                                            
aaCAAATATTTATAATCGTTTAAAAGATTTCACGAACCCAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATAGCGCCTGCa                                                                                         >  1:2644625/1‑100 (MQ=255)
                    tAAAAGATTTCACGAACCCAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCaat                                                                    <  1:3930664/100‑3 (MQ=255)
                              cACGAACCCAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTattgattg                                                          >  1:5889637/1‑100 (MQ=255)
                                     ccAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACt                                                   >  1:3722719/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAATCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:4364444/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:1669217/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:3049675/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:3649883/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:5400402/1‑100 (MQ=255)
                                      cAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTg                                                  >  1:5545466/1‑100 (MQ=255)
                                        gAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTGaa                                                >  1:4230058/1‑100 (MQ=255)
                                                                  gccGAATAATAAACCTAATGCGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTGAAATTTAATAAGACTAAATTTTTTGTAg                      >  1:639782/3‑100 (MQ=255)
                                                                                      cGATGGCG‑CTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTGAAATTTAATAAGACTAAATTTTTTGTAGCATCGTATGCTAAGccactg  <  1:1413308/100‑3 (MQ=255)
                                                                                              |                                                                                            
AACAAATATTTATAATCGTTTAAAAGATTTCACGAACCCAGAAACAATTAATTTGGAAATTTGGTCGGCGAATAATAAACCTAATGCGATGGCGCCTGCAATAAGTGTAACCTCTAGCATGGTATTGATTGCTGTACTGAAATTTAATAAGACTAAATTTTTTGTAGCATCGTATGCTAAGCCACCA  >  NC_007793/641572‑641758

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