Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 521,750 A→G R141R (AGA→AGG)  SAUSA300_RS02475 → DNA polymerase III subunit delta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793521,7500AG97.1% 83.6 / ‑5.7 34R141R (AGA→AGG) SAUSA300_RS02475DNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (33/0);  total (34/0)
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.00e+00

TGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGATGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAG  >  NC_007793/521662‑521845
                                                                                        |                                                                                               
tGAAATTTCTTGAAGAACCACCGGACAGTACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGt                                                                                      >  1:4740987/1‑100 (MQ=255)
tGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAgatcgg                                                                                      >  1:3817407/1‑95 (MQ=255)
          tGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAg                                                                            >  1:6035370/1‑100 (MQ=255)
          tGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAg                                                                            >  1:4577875/1‑100 (MQ=255)
             agaaCCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGccg                                                                         >  1:4416640/1‑99 (MQ=255)
               taccaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAt                                                                       >  1:2916846/2‑100 (MQ=255)
               caccaccGGACAATACGATTGCTATTTTATTGTCTACCAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAt                                                                       >  1:5798732/2‑100 (MQ=255)
               aaCCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAACCTAt                                                                       >  1:2916860/1‑100 (MQ=255)
                accaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGCCACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTAta                                                                      >  1:1759416/1‑99 (MQ=255)
                 ccaccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTg                                                                     >  1:5666075/1‑100 (MQ=255)
                   gccGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGAt                                                                   >  1:1957899/2‑100 (MQ=255)
                   accGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGAt                                                                   >  1:320581/1‑100 (MQ=255)
                                      aTTTTATTGTCT‑‑AAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGa                                              >  1:4371927/1‑100 (MQ=255)
                                               ccTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGGGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGACTAGtca                                       >  1:6374227/2‑98 (MQ=255)
                                                cTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGa                                      >  1:1234924/1‑100 (MQ=255)
                                                  gcAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:3390413/2‑100 (MQ=255)
                                                  aCAAAACTTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:2204299/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACTCAATCCATTCAAGATGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGACTCGCTGAAc                                    >  1:6508314/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAGTAGATTAGTTGAAc                                    >  1:4823821/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGTTTAgagcgga                                    >  1:3855344/1‑94 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:4482247/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:6389247/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:122117/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1162681/1‑100 (MQ=255)
                                                  aCAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGGGGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAAc                                    >  1:1210266/1‑100 (MQ=255)
                                                           gAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGt                           >  1:5142919/1‑100 (MQ=255)
                                                           gAGCAAATTTTAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGt                           >  1:98491/1‑100 (MQ=255)
                                                                     tAGACACAATCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGt                 >  1:1936187/1‑100 (MQ=255)
                                                                           gcctCCATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCCGa           >  1:1902451/4‑100 (MQ=255)
                                                                               acATTCAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGGAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATg       >  1:6107891/2‑100 (MQ=255)
                                                                                   tCAAGGCGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTa   >  1:4967019/1‑100 (MQ=255)
                                                                                    cAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAg  >  1:5200073/1‑100 (MQ=255)
                                                                                    cAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAg  >  1:6259297/1‑100 (MQ=255)
                                                                                    cAAGGTGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTCGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAg  >  1:3188207/1‑100 (MQ=255)
                                                                                        |                                                                                               
TGAAATTTCTTGAAGAACCACCGGACAATACGATTGCTATTTTATTGTCTACAAAACCTGAGCAAATTTTAGACACAATCCATTCAAGATGTCAGCATGTATATTTCAAGCCTATTGATAAAGAAAAGTTTATAAATAGATTAGTTGAACAAAACATGTCTAAGCCAGTAGCTGAAATGATTAG  >  NC_007793/521662‑521845

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: ‑