Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 1,954 2 bp→TGC intergenic (+187/+138) 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,9531.T91.7% 61.5 / ‑0.3 24intergenic (+186/+140)SAUSA300_RS14710/SAUSA300_RS14715membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (17/5);  total (19/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
*NC_0077901,9540AG88.0% 48.8 / 0.1 25intergenic (+187/+139)SAUSA300_RS14710/SAUSA300_RS14715membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (17/5);  total (20/5)
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.57e-01
*NC_0077901,9550TC88.0% 45.5 / ‑0.1 25intergenic (+188/+138)SAUSA300_RS14710/SAUSA300_RS14715membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (17/5);  total (20/5)
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.57e-01

GATGCTTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTC  >  NC_007790/1867‑2052
                                                                                        |||                                                                                                
gatgatTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga                                                                                         <  1:125353/100‑1 (MQ=255)
gatgatTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga                                                                                         >  1:1030939/1‑100 (MQ=255)
  tgattttGACGAGATGGAAGAAAAATATGATTATCATGAGGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAgaagaag                                                                                       >  1:1094791/4‑99 (MQ=255)
      tttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   >  1:4612656/1‑100 (MQ=255)
      tttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   >  1:3154310/1‑100 (MQ=255)
      cttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   <  1:3636334/99‑1 (MQ=255)
       ttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTa                                                                                  >  1:265200/1‑100 (MQ=255)
        tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaagatcggaaga                                                                                 >  1:1403690/1‑92 (MQ=255)
        tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:1658630/1‑100 (MQ=255)
        tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:730542/1‑100 (MQ=255)
        tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:2248526/1‑100 (MQ=255)
        tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAAT‑AAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAAc                                                                                >  1:734945/1‑100 (MQ=255)
         gACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAAc                                                                                >  1:2660365/1‑100 (MQ=255)
          ccGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATTCCGAGAAAATTTATTGTGCGTTGAGAAGAAACCTTAACt                                                                               >  1:5054260/2‑100 (MQ=255)
           cGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTa                                                                              >  1:2168409/1‑100 (MQ=255)
           cGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTa                                                                              >  1:5200938/1‑100 (MQ=255)
                 ggAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAATCTTg                                                                        >  1:5982382/1‑100 (MQ=255)
                  gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:2294083/100‑2 (MQ=255)
                  gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:3322774/100‑2 (MQ=255)
                  gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:2090425/100‑2 (MQ=255)
                                                         ccgacgCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCAt                               >  1:390854/2‑100 (MQ=255)
                                                               gccccgtGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACa                         >  1:1045771/7‑100 (MQ=255)
                                                                    aTGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGtt                     >  1:2800528/1‑100 (MQ=255)
                                                                    aTGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGtt                     >  1:301323/1‑100 (MQ=255)
                                                                                      g‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTc  >  1:2049674/1‑100 (MQ=255)
                                                                                        |||                                                                                                
GATGCTTTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTC  >  NC_007790/1867‑2052

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