Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,452,377 C→T 72.0% P181P (CCG→CCA aroC ← chorismate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,452,3770CT72.0% 33.8 / 15.5 25P181P (CCG→CCAaroCchorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (9/9);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGG  >  W3110S.gb/2452316‑2452447
                                                             |                                                                      
aGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCaaaaa                                                               <  1:2592039/71‑1 (MQ=255)
          ttttttCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGAtttt                                                     >  1:188290/1‑71 (MQ=255)
             tttCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTg                                                    <  1:2390508/69‑1 (MQ=255)
             tttCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGc                                                   <  1:1134119/70‑1 (MQ=255)
                     cacgcaTCCACCCGTCTAACGCGGCGATTTTGTCGGGGGCTGGGCAAAAAAACGGATTTTGCTCg                                                <  1:1569943/65‑1 (MQ=255)
                           tCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGc                                                   <  1:66779/56‑1 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGAtt                                                       >  1:1942381/1‑50 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTc                                  >  1:2339601/1‑71 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTc                                  >  1:2534306/1‑71 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTc                                  >  1:1821416/1‑71 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTc                                  >  1:1818075/1‑71 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGa                                               <  1:2186671/58‑1 (MQ=255)
                             aaCTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGt                                   >  1:1651055/1‑70 (MQ=255)
                                 cGTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTg                                                    <  1:2111556/49‑1 (MQ=255)
                                  gTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCttt                               <  1:2680227/69‑1 (MQ=255)
                                  gTCTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGa                             >  1:2641320/1‑71 (MQ=255)
                                    cTAACGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGata                           >  1:1043379/1‑71 (MQ=255)
                                    cTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGAttt                                                      >  1:2600639/1‑44 (MQ=255)
                                     tAACGCGTCGATTTTGTCGGGGTCCGGGCaaaaaa                                                              <  1:820742/35‑1 (MQ=255)
                                       aCGCGTCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCttt                               <  1:2185611/64‑1 (MQ=255)
                                            tCGATTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCgg                   <  1:1433858/71‑1 (MQ=255)
                                            tCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCgg                   <  1:494968/71‑1 (MQ=255)
                                               aTTTTGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCgg                   <  1:2271927/68‑1 (MQ=255)
                                                 tttGTCGGGGTCTGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGt              >  1:815709/1‑71 (MQ=255)
                                                   tGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTc             <  1:2644845/70‑1 (MQ=255)
                                                             cGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTggg  >  1:2484417/1‑71 (MQ=255)
                                                             |                                                                      
AGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGG  >  W3110S.gb/2452316‑2452447

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