Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,036,248 G→A D232D (GAC→GAT)  SAUSA300_RS10265 ← DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,036,2480GA97.8% 119.5 / ‑5.1 46D232D (GAC→GAT) SAUSA300_RS10265DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (5/39);  total (5/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

ATCTCTTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGGTCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGTTA  >  NC_007793/2036162‑2036344
                                                                                      |                                                                                                
aTCTCTTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTg                                                                                     <  1:3785070/100‑1 (MQ=255)
  ctctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTTTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGcc                                                                                   <  1:72281/100‑1 (MQ=255)
  ctctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGtcgg                                                                                   <  1:1961233/100‑4 (MQ=255)
  ctctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGcc                                                                                   <  1:5469239/100‑1 (MQ=255)
  ctctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGcc                                                                                   <  1:3018757/100‑1 (MQ=255)
  ctctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGAACAATACCGCGTGcc                                                                                   <  1:5850312/100‑1 (MQ=255)
   ttccgatctgcatTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAGTACCGCGTGCCt                                                                                  <  1:5960394/90‑1 (MQ=255)
    ctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCGCTTCACTATGATCAATACCGCGTGCCtt                                                                                 <  1:3434735/100‑1 (MQ=255)
    ctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCCCTATGATCAATACCGCGTGCCtt                                                                                 <  1:2027099/100‑1 (MQ=255)
    ctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCtt                                                                                 <  1:2837537/100‑1 (MQ=255)
    ctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCta                                                                                 <  1:254756/100‑2 (MQ=255)
    ctTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCta                                                                                 <  1:6501882/100‑2 (MQ=255)
     ttCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCtta                                                                                <  1:1687389/100‑2 (MQ=255)
      tcgtcatcatTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGt                                                                               <  1:114413/97‑1 (MQ=255)
       catcatcatTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGtt                                                                              <  1:5498453/100‑1 (MQ=255)
       catcatcatTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGtt                                                                              <  1:3984718/100‑1 (MQ=255)
         ccgatctctaCGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTtat                                                                            <  1:4272222/91‑1 (MQ=255)
           atcttTTACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTtatac                                                                          <  1:4007998/100‑2 (MQ=255)
                   cGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCg                                                                  <  1:2934899/100‑1 (MQ=255)
                     tCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGCTCAATACCGCGTTCCTTGTTATATAAAcccctag                                                                <  1:1690060/100‑4 (MQ=255)
                     tCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGag                                                                <  1:500206/100‑2 (MQ=255)
                      cTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGAcc                                                               >  1:3019592/1‑100 (MQ=255)
                            cAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCtttt                                                         <  1:5728026/100‑1 (MQ=255)
                                tgtgCGTTCAGTCCCTACTGATTTTCTTACTCTTGATTATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCca                                                     <  1:5379667/100‑1 (MQ=255)
                                tgtgCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCca                                                     <  1:6627628/100‑1 (MQ=255)
                                cttccgatctgTCCTTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCca                                                     <  1:4848843/90‑1 (MQ=255)
                                 gtgCGTTCAGTCCCTACTGATGTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCcaa                                                    >  1:795818/1‑100 (MQ=255)
                                  tgCGTTCAGTCCCTACTGATTTTCTTACTCTTTATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCcaaa                                                   <  1:5234619/100‑1 (MQ=255)
                                         aGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGt                                            <  1:1207141/100‑1 (MQ=255)
                                            cccTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTAtt                                         <  1:4133353/100‑1 (MQ=255)
                                            cccTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTAtt                                         <  1:1393712/100‑1 (MQ=255)
                                             cttccgatctttCTTGCTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTa                                        <  1:3212056/91‑1 (MQ=255)
                                             ccTACTGATTTTCTCACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTa                                        <  1:4818107/100‑1 (MQ=255)
                                                aCTGATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAAtt                                     >  1:5502758/1‑100 (MQ=255)
                                                   gATTTTCTTACTCTTGATGATTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTGTTAATTCaa                                  >  1:4248500/1‑100 (MQ=255)
                                                          ttACTCTTGATGATCTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCgg                           <  1:4843909/100‑2 (MQ=255)
                                                                    tgatTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCatat                 <  1:6001810/100‑1 (MQ=255)
                                                                     gatTTCACTTCACTATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCatagg                <  1:5544864/100‑3 (MQ=255)
                                                                      atTTCACTTCACTATGGTCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATaa               <  1:1180760/100‑1 (MQ=255)
                                                                        ttcacttcacGATGATCAATACCGCGTGCCTTGTTATATAAACCCCGGCCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCAtatagga             <  1:6550650/100‑4 (MQ=255)
                                                                               cactATGATCAATACCGCGTGCCTTGTTATATAAACCCCGGCCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTAcc      <  1:378912/100‑1 (MQ=255)
                                                                                  tATGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCggc   <  1:5144072/100‑3 (MQ=255)
                                                                                   aTGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGTTa  <  1:6485113/100‑1 (MQ=255)
                                                                                   aTGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGTTa  <  1:1333334/100‑1 (MQ=255)
                                                                                   aTGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGTTa  >  1:1307417/1‑100 (MQ=255)
                                                                                   aTGATCAATACCGCGTGCCTTGTTATATAAACCCCGACCGCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGtgg  <  1:3080852/100‑3 (MQ=255)
                                                                                      |                                                                                                
ATCTCTTCATCATCATTCACGTCTGTTGCAAATGTGCGTTCAGTCCCTACTGATTTTCTTACTCTTGATGATTTCACTTCACTATGGTCAATACCGCGTGCCTTGTTATATAAACCCCGACCTCTTTTTCCAAACAAACGTATTAATTCAAATTCCGTTTTCTCATATAAATCTCTACCGTTA  >  NC_007793/2036162‑2036344

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