Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,499,109 A→T N122I (AAT→ATT)  SAUSA300_RS07270 → sulfite reductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,499,1090AT93.3% 31.2 / ‑4.6 15N122I (AAT→ATT) SAUSA300_RS07270sulfite reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (12/2);  minor base C (1/0);  total (13/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 = 7.10e-01

ACACTGAGTGATATGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTAATGTGATGCGAAGTATACATAGTTCTTTGTTTACTGAAAATACTATTATTAATTATTTAGTAGAGTTAACCGCATGTTGGATGATA  >  NC_007793/1499015‑1499194
                                                                                              |                                                                                     
acacTGAGTGATATGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGa                                                                                  >  1:6411639/1‑100 (MQ=255)
 tacTGAGTGATATGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGAt                                                                                 >  1:398226/2‑100 (MQ=255)
      aGTGATATGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGaa                                                                            >  1:504496/1‑100 (MQ=255)
          gtatGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGtat                                                                        >  1:1456323/2‑100 (MQ=255)
            atGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATAc                                                                      >  1:3373141/1‑100 (MQ=255)
            atGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATAc                                                                      >  1:5263916/1‑100 (MQ=255)
            atGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTACTGGGATCGGAAgagcac                                                                      >  1:4946957/1‑95 (MQ=255)
              ggCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATACAt                                                                    >  1:923441/1‑100 (MQ=255)
                     caTCAACGTTACGTCCTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCtt                                                             >  1:4253853/2‑100 (MQ=255)
                     aaTCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCtt                                                             >  1:5051268/1‑100 (MQ=255)
                       aacaCGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCTTTg                                                           >  1:3462149/4‑100 (MQ=255)
                               aCGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCTTTGTTTACTGa                                                   >  1:21938/1‑100 (MQ=255)
                                                                        aaTGCGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCTTTGTTTACTGAAAATACTATTATTAGTTATTTAGTAGAGTTAACCGCATGtg          <  1:3771793/100‑2 (MQ=255)
                                                                            cGAAACATACGTATTCTATTGTGATGCGAAGTATACATAGTTCTTTGTTTACTGAAAATACTATTATTAATTATTTAGTAGAGTTAACCGCATGTTGgag      <  1:1950630/100‑2 (MQ=255)
                                                                                aCATACGTATTCTATTGTGATGCGAAGTATACATAGTTCTTTGTTTACTGAAAATACTATTATTAATTATTTAGTAGAGTTAACCGCATGTTGGATGata  >  1:1812728/1‑100 (MQ=255)
                                                                                              |                                                                                     
ACACTGAGTGATATGGCAGGAAATCAACGTTACGTCTTTCTAGATCATAACAATCAATTTGTTGCAGATCAAAATGCGAAACATACGTATTCTAATGTGATGCGAAGTATACATAGTTCTTTGTTTACTGAAAATACTATTATTAATTATTTAGTAGAGTTAACCGCATGTTGGATGATA  >  NC_007793/1499015‑1499194

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