Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,222,206 C→T 78.1% N306N (AAC→AAT ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,222,2060CT78.1% 126.0 / 34.9 64N306N (AAC→AATebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (4/10);  new base T (29/21);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 7.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGG  >  W3110S.gb/3222142‑3222272
                                                                |                                                                  
aaTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTg                                                              >  1:1245572/1‑71 (MQ=255)
aaTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTg                                                              >  1:1516366/1‑71 (MQ=255)
 aTGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTg                                                              >  1:183686/1‑70 (MQ=255)
   ggTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTgg                                                             >  1:134553/1‑69 (MQ=255)
    gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAg                                                          >  1:1298357/1‑71 (MQ=255)
    gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAg                                                          >  1:2102518/1‑71 (MQ=255)
     tCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgt                                                         <  1:2371451/71‑1 (MQ=255)
     tCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgt                                                         <  1:766156/71‑1 (MQ=255)
     tCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgt                                                         <  1:2571552/71‑1 (MQ=255)
     tCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAgt                                                         <  1:794661/71‑1 (MQ=255)
             aTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCcac                                                 <  1:260729/71‑1 (MQ=255)
             aTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCcac                                                 <  1:1918791/71‑1 (MQ=255)
                cccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACaa                                               >  1:2283279/1‑70 (MQ=255)
                cccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACaa                                               >  1:2662374/1‑70 (MQ=255)
                  cTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAAcgc                                            <  1:1979534/71‑1 (MQ=255)
                   ttATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCcac                                                 <  1:1538745/65‑1 (MQ=255)
                    tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGt                                          >  1:421427/1‑71 (MQ=255)
                    tATCTTTACCATCTGGTCATGACGCTGAAAGACGCAAACGGCAATGTTCTGGAAgtggt                                                      >  1:1166943/1‑59 (MQ=255)
                     aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGtt                                         >  1:1977524/1‑71 (MQ=255)
                     aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGtt                                         >  1:2099967/1‑71 (MQ=255)
                     aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGtt                                         >  1:897725/1‑71 (MQ=255)
                      tCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTg                                        >  1:851158/1‑71 (MQ=255)
                      tCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTg                                        >  1:865148/1‑71 (MQ=255)
                       cTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgtggt                                                      <  1:582414/56‑1 (MQ=255)
                       cTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgtgg                                                       >  1:356487/1‑55 (MQ=255)
                       cTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTgg                                       >  1:2495174/1‑71 (MQ=255)
                            ccATCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAAcgcg                                           <  1:164489/62‑1 (MQ=255)
                              aTCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:699609/71‑1 (MQ=255)
                              aTCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:2789981/71‑1 (MQ=255)
                              aTCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:2751126/71‑1 (MQ=255)
                              aTCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:889956/71‑1 (MQ=255)
                              aTCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:2497897/71‑1 (MQ=255)
                               tCTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGt                                <  1:2461778/70‑1 (MQ=255)
                               tCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                  <  1:1618834/36‑1 (MQ=255)
                               tCTGGTCATGACGCTGAAAGACGCCAACGGCAACGt                                                                  <  1:766045/36‑1 (MQ=255)
                               tCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAAcgcg                                           <  1:864375/59‑1 (MQ=255)
                               tCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAAcgcg                                           <  1:2597663/59‑1 (MQ=255)
                                cTGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGc                                                    <  1:273590/49‑1 (MQ=255)
                                 tGGTCATGACGCTGATAGACGCCAACGGCAACGTTc                                                                >  1:399004/1‑36 (MQ=255)
                                 tGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGa                              >  1:2395619/1‑70 (MQ=255)
                                 tGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGa                              >  1:1547657/1‑70 (MQ=255)
                                 tGGTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGa                              >  1:2079899/1‑70 (MQ=255)
                                   gTCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGatat                           >  1:2101167/1‑71 (MQ=255)
                                    tCATGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAgt                                                         <  1:1785399/40‑1 (MQ=255)
                                     cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGt                                          >  1:2359291/1‑54 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         >  1:418714/1‑70 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         >  1:2723452/1‑70 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         >  1:52506/1‑70 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         <  1:1658140/70‑1 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         >  1:710876/1‑70 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCa                         >  1:232591/1‑70 (MQ=255)
                                      aTGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGATGGCTTCCGTGATATCaa                        <  1:680475/71‑1 (MQ=255)
                                       tGACGCTGAAAGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCaa                        <  1:1280738/70‑1 (MQ=255)
                                       tGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAAcgcg                                           <  1:1233277/51‑1 (MQ=255)
                                         aCGCTGAAAGACGCCAACGGCAATGTTCTGGAAgtggt                                                      >  1:973642/1‑38 (MQ=255)
                                         aCGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                      <  1:2201221/70‑1 (MQ=255)
                                                 aGACGCCAACGGCAATGTTCTGGAAGTTGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACg             <  1:1072093/71‑1 (MQ=255)
                                                 aGACGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACg             <  1:1284244/71‑1 (MQ=255)
                                                    cGCCAACGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGTCTTCCGTGATATCAAAGTGCGCGACGGTc          >  1:1274283/1‑71 (MQ=255)
                                                         aCGGCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTc     >  1:810584/1‑71 (MQ=255)
                                                           ggCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTg   >  1:2332969/1‑71 (MQ=255)
                                                           ggCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTg   >  1:2786217/1‑71 (MQ=255)
                                                           ggCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTg   <  1:2323031/71‑1 (MQ=255)
                                                            gCAATGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTgg  <  1:370100/71‑1 (MQ=255)
                                                                |                                                                  
AATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGG  >  W3110S.gb/3222142‑3222272

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