Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 1,871 C→A intergenic (+104/+222) 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,8710CA97.3% 98.4 / ‑4.4 37intergenic (+104/+222)SAUSA300_RS14710/SAUSA300_RS14715membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (34/2);  minor base T (1/0);  total (35/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.87e-01

TGAGGGGATTTTGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGCTTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGA  >  NC_007790/1774‑1965
                                                                                                  |                                                                                                  
tGAGAGGATTTTGAAACGAGTTTCTTCTTGTTTTCGCACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAtt                                                                                                  >  1:3758026/1‑100 (MQ=255)
  ccaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAttct                                                                                                >  1:5878961/4‑98 (MQ=255)
   aaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAttgag                                                                                               >  1:4063017/3‑97 (MQ=255)
          ctGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑gagga                                                                                      >  1:1628279/2‑98 (MQ=255)
           tGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGGGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAtc                                                                                     >  1:2065016/1‑99 (MQ=255)
                       cttcttGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTCG‑GAA‑GAGATCGgaagagcacac                                                                         >  1:3980210/1‑94 (MQ=255)
                         tcttGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTAT‑GAC‑GAGATCGgaagagcacacgt                                                                       >  1:198960/1‑92 (MQ=255)
                         tcttGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTAT‑AAC‑GAGATCGgaagagcacacgt                                                                       >  1:4808686/1‑92 (MQ=255)
                          cttGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GCA‑GAGATCGgaagagcacacgtc                                                                      >  1:1294186/1‑91 (MQ=255)
                            tGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGGGC‑GAGATCGgaagagcacacgtct                                                                     >  1:174819/1‑90 (MQ=255)
                            tGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GACTGAGATCGgaagagcacacgtct                                                                     >  1:1423139/1‑90 (MQ=255)
                            tGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTa                                                                   >  1:5084974/1‑100 (MQ=255)
                               tttCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCa                                                                >  1:265122/1‑100 (MQ=255)
                               ttccacaCTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCa                                                                >  1:1278699/4‑100 (MQ=255)
                                ttCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCAt                                                               >  1:5103095/1‑100 (MQ=255)
                                  cacaCTGTTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAT‑GAGAT‑GGAAGAAAAATATGATTATCATg                                                              >  1:1684560/1‑100 (MQ=255)
                                  cacaCTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGa                                                             >  1:3140259/1‑100 (MQ=255)
                                   gcacTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGaa                                                            >  1:6064712/2‑100 (MQ=255)
                                    cacTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAg                                                           >  1:3668286/1‑100 (MQ=255)
                                                  tgccTATTGGTGTTGTTGCTTACTTTTTGTTTTTTTTTAAAGATGATGTTTTT‑GCC‑GAGAT‑GGAAGAAAATTTTAATAATCTTGAGGTTAATAAAAAATa                                              >  1:3774422/3‑100 (MQ=255)
                                                  gtgcTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGACGAAAAATATGATTATCATGAAGATAATAAAAAATa                                              >  1:6101819/4‑100 (MQ=255)
                                                       aTTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAgacga                                         >  1:3166623/1‑100 (MQ=255)
                                                             gttgttGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCAtt                                   >  1:1869299/1‑100 (MQ=255)
                                                                gttgCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTAt                                >  1:1976042/1‑100 (MQ=255)
                                                                  tgCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGc                              >  1:5959963/1‑100 (MQ=255)
                                                                   gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGGCGACGCATTTATGcc                             >  1:1183370/1‑100 (MQ=255)
                                                                    tgtACTTTTTGTTTTTCTATAAAGATGCTGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCg                            >  1:5939160/3‑100 (MQ=255)
                                                                    cTTACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCg                            >  1:4425457/1‑100 (MQ=255)
                                                                      tACTTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GCGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCgag                          >  1:1861475/1‑100 (MQ=255)
                                                                      tACCTTTTGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCgag                          >  1:2275987/1‑100 (MQ=255)
                                                                          ttttGTTTTTCTATAAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGaaaa                      >  1:2383270/1‑100 (MQ=255)
                                                                                      ataAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgt          >  1:3894987/1‑95 (MQ=255)
                                                                                      ataAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgt          >  1:4109451/1‑95 (MQ=255)
                                                                                      ataAAGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgt          >  1:4447403/1‑95 (MQ=255)
                                                                                        aaaGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:3822129/100‑8 (MQ=255)
                                                                                        aaaGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgtcg        <  1:2906812/100‑8 (MQ=255)
                                                                                         aaGATGATGATTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgtaga       >  1:652072/1‑92 (MQ=255)
                                                                                           gatgatgatTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcatagatc     >  1:4225092/1‑90 (MQ=255)
                                                                                              gatgatTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga  <  1:125353/100‑1 (MQ=255)
                                                                                              gatgatTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga  >  1:1030939/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TGAGGGGATTTTGAAACGAGTTTCTTCTTGTTTTCACACTG‑TTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGCTTTT‑GAC‑GAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGA  >  NC_007790/1774‑1965

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