Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 1,967 C→G intergenic (+200/+126) SAUSA300_RS14710 → / ← SAUSA300_RS14715 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077901,9670CG88.0% 54.3 / 2.7 25intergenic (+200/+126)SAUSA300_RS14710/SAUSA300_RS14715membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base G (20/2);  total (23/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 = 1.30e-01

TTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑A‑TGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTT  >  NC_007790/1872‑2034
                                                                                                 |                                                                   
ttttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCt                                                                  >  1:3571407/1‑100 (MQ=255)
  ttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTa                                                                >  1:1424125/1‑100 (MQ=255)
  ttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTa                                                                >  1:4506837/1‑100 (MQ=255)
  ttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTa                                                                >  1:5018753/1‑100 (MQ=255)
   tGACGAGATGGAAGAAAAATATGATTATCATGAGGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTaa                                                               >  1:4822400/1‑100 (MQ=255)
   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTaa                                                               >  1:1327063/1‑100 (MQ=255)
   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTaa                                                               >  1:3358931/1‑100 (MQ=255)
   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTaa                                                               >  1:3143780/1‑100 (MQ=255)
   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTaa                                                               >  1:1157575/1‑100 (MQ=255)
    gccGAGATGGAAGAAAAATATGATTATCAGGAGGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAgcccccag                                                              >  1:2158383/3‑95 (MQ=255)
    gACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAAc                                                              >  1:1590955/1‑100 (MQ=255)
      atacctGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCGTTAACTa                                                            >  1:2201106/6‑100 (MQ=255)
            ggAAGAAAAATATGATTATCATGAAGTTAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTg                                                      >  1:4763828/1‑100 (MQ=255)
            ggAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGCTGAGAAGAAGCCTTAACTAAACTTg                                                      >  1:1242376/1‑100 (MQ=255)
             gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACtaag                                                     <  1:3953809/100‑4 (MQ=255)
              aagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGGGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTgga                                                    >  1:3555115/1‑98 (MQ=255)
              aagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTgga                                                    >  1:2341119/1‑98 (MQ=255)
               agaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTGGAg                                                   <  1:1465436/100‑1 (MQ=255)
                                   ctgATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑TGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTg                               >  1:2416009/3‑100 (MQ=255)
                                               caTAGACGACGCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGcc                   >  1:4129629/2‑100 (MQ=255)
                                                   gacgacGCATTTATGCCGAGAAAATTTATTG‑A‑TGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCa              >  1:1814264/1‑100 (MQ=255)
                                                    ccgacgCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCa              >  1:1545002/2‑100 (MQ=255)
                                                    ccgacgCATTTATGCCGAGAAAATTTATTG‑TGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCa              >  1:2804317/2‑100 (MQ=255)
                                                    acgacgCATTTATGCCGAGAAAATTTATTG‑A‑TGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCAt             >  1:686590/1‑100 (MQ=255)
                                                       acgCATTTATGCCGAGAAAATTTATTGTG‑CGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTc           >  1:3051170/1‑100 (MQ=255)
                                                               aTGCCGAGAAAATTTATTG‑A‑TGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGGGCTATTAAGCCGACCATTCGACAAGttt  >  1:1793070/1‑100 (MQ=255)
                                                                                                 |                                                                   
TTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑A‑TGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTT  >  NC_007790/1872‑2034

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