Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 50,527 C→T 8.5% E232E (GAG→GAA apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb50,5270CT 8.5% 132.6 / 3.1 48E232E (GAG→GAAapaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base T (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCAC  >  W3110S.gb/50459‑50596
                                                                    |                                                                     
caggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCt                                                                      <  1:3966419/70‑1 (MQ=255)
caggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTc                                                                     <  1:3087548/71‑1 (MQ=255)
 aggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTc                                                                     >  1:1802319/1‑70 (MQ=255)
    ttAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGc                                                                 <  1:946669/71‑1 (MQ=255)
     tAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGc                                                                 >  1:3221188/1‑70 (MQ=255)
      aaTGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCg                                                                <  1:198773/70‑1 (MQ=255)
            ccaccCCAGCAGCAGACGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGa                                                               >  1:1662666/1‑65 (MQ=255)
                    gcagcagcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCa                                                  >  1:633231/1‑70 (MQ=255)
                        cagcCGTTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:3365430/71‑1 (MQ=255)
                        cagcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:340133/71‑1 (MQ=255)
                         agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:4374201/70‑1 (MQ=255)
                         agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             >  1:375547/1‑70 (MQ=255)
                         agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             >  1:2410122/1‑70 (MQ=255)
                            cGGTATCCAGCGCGTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATg                                         >  1:657565/1‑71 (MQ=255)
                              gTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGc                                        >  1:1456756/1‑70 (MQ=255)
                              gTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCt                                       >  1:1677056/1‑71 (MQ=255)
                              gTATCCAGCGCGTATATACCTTCAGGCGTACCTTTGCCCTCCAGCGa                                                               >  1:405343/1‑47 (MQ=255)
                               tATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                      <  1:717163/71‑1 (MQ=255)
                                  ccAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAt                                   >  1:3211087/1‑71 (MQ=255)
                                   cAGCGCGTAAATACCTTCCGGCGTACCTTTGCCTTcc                                                                    >  1:740868/1‑37 (MQ=37)
                                    aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTg                                                      <  1:449598/50‑1 (MQ=255)
                                    aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCaa                                                 <  1:3650893/55‑1 (MQ=255)
                                    aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCaa                                                 <  1:1482152/55‑1 (MQ=255)
                                    aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAttc                                 >  1:4336703/1‑71 (MQ=255)
                                      cgcgTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAttctt                               <  1:4341370/71‑1 (MQ=255)
                                         gTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:19926/54‑1 (MQ=39)
                                         gTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:762887/54‑1 (MQ=39)
                                          tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:527912/53‑1 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCGTCCCGCGATGCCCAGTGACCAAAGGCGa                                           >  1:408695/1‑53 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGtcc                                                          <  1:3203038/38‑1 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCa                                                  <  1:467169/46‑1 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa                                    <  1:714358/60‑1 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa                                    <  1:4418000/60‑1 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa                         >  1:309818/1‑71 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa                         >  1:2920827/1‑71 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa                         >  1:2265887/1‑71 (MQ=255)
                                            tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa                         >  1:4415008/1‑71 (MQ=255)
                                              taCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa                         >  1:121667/1‑69 (MQ=255)
                                               aCCTGCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAt                                   <  1:2806392/58‑1 (MQ=255)
                                                  ttCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                                             <  1:3062820/45‑1 (MQ=255)
                                                    ccGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGc                   >  1:3486854/1‑69 (MQ=255)
                                                      ggCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAgg               <  1:840020/71‑1 (MQ=255)
                                                          tACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAgg               <  1:482816/67‑1 (MQ=255)
                                                             cTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCa        <  1:4092305/71‑1 (MQ=255)
                                                               ttGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGaa             >  1:2353053/1‑64 (MQ=255)
                                                               ttGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGaa             >  1:4055750/1‑64 (MQ=255)
                                                               ttGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCaaa      >  1:1403434/1‑71 (MQ=255)
                                                                tGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAAc     >  1:366495/1‑71 (MQ=255)
                                                                   ccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCa        <  1:3159028/65‑1 (MQ=255)
                                                                   ccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCAc  <  1:583440/71‑1 (MQ=255)
                                                                   ccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCAc  <  1:2416484/71‑1 (MQ=255)
                                                                    |                                                                     
CAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCAC  >  W3110S.gb/50459‑50596

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: