Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,017,703 T→C 45.5% N1175D (AAC→GAC)  rhsB ← rhsB element core protein RshB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,017,7030TC45.5% 4.6 / 43.5 33N1175D (AAC→GAC) rhsBrhsB element core protein RshB
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base C (3/12);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 7.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAG  >  W3110S.gb/4017635‑4017757
                                                                    |                                                      
atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAACTGATGCGGGTTCTCTTCATTCAGCAGGTCGc                                                      >  1:1099927/1‑71 (MQ=255)
   atACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGccc                                                    >  1:3907729/1‑70 (MQ=35)
         ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtc                                               <  1:831073/69‑1 (MQ=35)
              ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCata                                         >  1:928128/1‑70 (MQ=40)
              ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCata                                         >  1:1380843/1‑70 (MQ=40)
                 gCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt                                     >  1:377114/1‑71 (MQ=39)
                   aGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTgc                                   >  1:3588097/1‑71 (MQ=255)
                      cGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca                                           <  1:189636/60‑1 (MQ=31)
                        gATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcac                               <  1:3652700/70‑1 (MQ=35)
                        gATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTTCGcacc                              <  1:2101466/71‑1 (MQ=255)
                             gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt                         <  1:3004132/71‑1 (MQ=40)
                              ctgctgCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTg                        <  1:5140168/71‑1 (MQ=255)
                                 ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCata                                         <  1:2391709/51‑1 (MQ=35)
                                    cAGCTGATGCGGCTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc                    <  1:2542781/69‑1 (MQ=255)
                                     aGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcccc                 <  1:5062214/71‑1 (MQ=255)
                                         gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTcc             >  1:3159162/1‑71 (MQ=39)
                                          aTGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCg            <  1:5132986/71‑1 (MQ=39)
                                          aTGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcc                   <  1:609819/64‑1 (MQ=255)
                                           tGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTc                                      <  1:4061411/44‑1 (MQ=255)
                                              gggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCt        <  1:1076901/71‑1 (MQ=255)
                                              gggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCt        <  1:2632357/71‑1 (MQ=255)
                                                 ttCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTc                                      <  1:3621687/38‑1 (MQ=40)
                                                  tctcTTCATTCAGCAGGTCGCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgat     >  1:4463963/1‑69 (MQ=14)
                                                   ctctTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatg    <  1:295420/70‑1 (MQ=255)
                                                   ctctTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatg    <  1:4956163/70‑1 (MQ=255)
                                                      ttcattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTg          <  1:2364806/61‑1 (MQ=255)
                                                      ttcattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCt        <  1:1873844/63‑1 (MQ=255)
                                                        cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTg          >  1:4488858/1‑59 (MQ=35)
                                                        cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTg          >  1:3633454/1‑59 (MQ=35)
                                                        cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTg          >  1:1288828/1‑59 (MQ=35)
                                                        cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCt        >  1:4755283/1‑61 (MQ=35)
                                                        cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCt        >  1:1620390/1‑61 (MQ=35)
                                                          ttcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAg  <  1:2266019/65‑1 (MQ=14)
                                                                    |                                                      
ATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAG  >  W3110S.gb/4017635‑4017757

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