Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 346,362 C→T 11.5% E436K (GAA→AAA)  prpR ← DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb346,3620CT11.5% 148.1 / 9.4 52E436K (GAA→AAA) prpRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base T (3/3);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTCGTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGCGCCGCCAGAGACACTTTCA  >  W3110S.gb/346296‑346430
                                                                  |                                                                    
gTCGTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTc                                                                      <  1:1900274/67‑1 (MQ=255)
gTCGTTCCATCATATTGCGCAGTTAACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTc                                                                      <  1:2924786/67‑1 (MQ=255)
  cGTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCAc                                                                    >  1:1983412/1‑67 (MQ=255)
   gTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGcc                                                               >  1:910008/1‑71 (MQ=255)
   gTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGcc                                                               >  1:632715/1‑71 (MQ=255)
   gTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGcc                                                               >  1:1874/1‑71 (MQ=255)
         tcatATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACtt                                                                  <  1:400868/62‑1 (MQ=255)
         tcatATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACtt                                                                  <  1:3384084/62‑1 (MQ=255)
         tcatATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACtt                                                                  <  1:1689212/62‑1 (MQ=255)
           atatTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACTTGc                                                                >  1:3069146/1‑62 (MQ=255)
           atatTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACTTGc                                                                >  1:62338/1‑62 (MQ=255)
           atatTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTTACTTGc                                                                >  1:2888530/1‑62 (MQ=255)
              ttGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTg                                                    >  1:2047204/1‑71 (MQ=255)
                  gcAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGt                                                <  1:2488659/71‑1 (MQ=255)
                      ttCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTc                                                                      <  1:273326/45‑1 (MQ=255)
                      ttCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGt                                                <  1:193403/67‑1 (MQ=255)
                        cACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATg                                            >  1:1518509/1‑69 (MQ=255)
                             atatTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGa                                     <  1:3272948/71‑1 (MQ=255)
                              tatTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGt                                                <  1:1379447/59‑1 (MQ=255)
                              tatTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaa                                    <  1:2962895/71‑1 (MQ=255)
                              tatTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaa                                    <  1:3424512/71‑1 (MQ=255)
                              tatTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaa                                    <  1:2298767/71‑1 (MQ=255)
                              tatTGCCTGGCCAGTCGTAGGGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaa                                    <  1:3044306/71‑1 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACtt                                                                  >  1:3588722/1‑39 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACtt                                                                  >  1:1503469/1‑39 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTaa                                              >  1:2897870/1‑59 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTaa                                              >  1:2954945/1‑59 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTaa                                              >  1:3721258/1‑59 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATgcag                                         >  1:106749/1‑64 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATgcag                                         >  1:3104429/1‑64 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATgcag                                         >  1:3333282/1‑64 (MQ=255)
                                ttGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaaaa                                  >  1:620211/1‑71 (MQ=255)
                                   ccTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATgg                               <  1:3358839/71‑1 (MQ=255)
                                       gCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCg                           <  1:3526650/71‑1 (MQ=255)
                                        ccAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaa                                    <  1:319568/61‑1 (MQ=255)
                                              gTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAgcg                    <  1:722500/71‑1 (MQ=255)
                                              gTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAgcg                    <  1:529977/71‑1 (MQ=255)
                                                aGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATg                                            <  1:1166777/45‑1 (MQ=255)
                                                    cagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:3448451/1‑69 (MQ=255)
                                                    cagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:1663760/1‑69 (MQ=255)
                                                    cagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:1929512/1‑69 (MQ=255)
                                                    cagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:3240099/1‑69 (MQ=255)
                                                    cagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:2358632/1‑69 (MQ=255)
                                                    aagcagcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGcgccgc                >  1:4014739/2‑69 (MQ=255)
                                                         gcaCAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGAGAGAGCGCCGCCAGAGaca        <  1:1674142/71‑1 (MQ=255)
                                                          cacaGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGCGCCGCCAGAGaca        <  1:1827859/71‑1 (MQ=255)
                                                               tttCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaaa                                   <  1:3391126/39‑1 (MQ=255)
                                                               tttCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaaa                                   <  1:844945/39‑1 (MQ=255)
                                                               tttCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGaaa                                   <  1:528144/39‑1 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAAt                                 >  1:540182/1‑39 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAAt                                 >  1:2032200/1‑39 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAAt                                 >  1:1938352/1‑39 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATggg                              >  1:3654887/1‑42 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGagag                      >  1:3793490/1‑50 (MQ=255)
                                                                 tCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGCGCCGCCAGAGACACTTTCa  >  1:665408/1‑70 (MQ=255)
                                                                  |                                                                    
GTCGTTCCATCATATTGCGCAGTTCACGAATATTGCCTGGCCAGTCGTAGTGCAGCAGCACAGTTTCACTTGCCTGTAACCCCTGGCGTAATGCAGCAGAAAATGGGGCGGAGAGCGCCGCCAGAGACACTTTCA  >  W3110S.gb/346296‑346430

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