Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,462,102 A→G 10.8% F53L (TTC→CTC)  sixA ← phosphohistidine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,462,1020AG10.8% 166.4 / 12.8 65F53L (TTC→CTC) sixAphosphohistidine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base G (2/5);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

AACTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTCG  >  W3110S.gb/2462034‑2462172
                                                                    |                                                                      
aaCTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCg                                                                      <  1:1540475/71‑1 (MQ=255)
     cTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGcg                                                                  >  1:3065751/1‑69 (MQ=255)
     cTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGc                                                                   <  1:2207368/69‑1 (MQ=255)
     cTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCt                                                                  >  1:2009435/1‑70 (MQ=255)
           tGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAg                                                            >  1:29199/1‑70 (MQ=255)
             gAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGaa                                                          >  1:112186/1‑70 (MQ=255)
             gAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGaa                                                          >  1:4882251/1‑70 (MQ=255)
               aGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAAcgc                                                       <  1:2248573/71‑1 (MQ=255)
                   aGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACTCGTTc                                                   <  1:1291985/71‑1 (MQ=255)
                   aGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:3211192/71‑1 (MQ=255)
                           aCAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCa                                           >  1:3875117/1‑71 (MQ=255)
                           aCAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:998194/63‑1 (MQ=255)
                           aCAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:2191408/63‑1 (MQ=255)
                           aCAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:4893146/63‑1 (MQ=255)
                            cAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   >  1:3592982/1‑62 (MQ=255)
                            cAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   >  1:3184606/1‑62 (MQ=255)
                            cAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   >  1:4778209/1‑62 (MQ=255)
                             aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACt                                         <  1:1502586/71‑1 (MQ=255)
                             aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACt                                         <  1:4908335/71‑1 (MQ=255)
                             aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACt                                         <  1:2764399/71‑1 (MQ=255)
                             aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCa                                                                >  1:3189249/1‑48 (MQ=255)
                             aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCa                                                                >  1:5123284/1‑48 (MQ=255)
                                ccccTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACtttt                                      >  1:4067070/1‑71 (MQ=255)
                                   cTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGa                                                           >  1:4302797/1‑47 (MQ=255)
                                       ttcttcCAGTGTTTGCTCGGCTCTCAGGAACGGGCTCACCAGAACGCGTTCg                                                  <  1:1198735/52‑1 (MQ=255)
                                       ttcttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTcagcag                                                            <  1:2971471/42‑1 (MQ=255)
                                       ttcttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAg                                                            <  1:1297332/42‑1 (MQ=255)
                                         cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGaa                                                          >  1:435764/1‑42 (MQ=255)
                                         cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGaa                                                          >  1:1969354/1‑42 (MQ=255)
                                         cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:3280542/49‑1 (MQ=255)
                                         cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGAc                                  >  1:2842153/1‑66 (MQ=255)
                                         cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                              <  1:1818203/70‑1 (MQ=255)
                                            ccAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGc                          >  1:4277926/1‑71 (MQ=255)
                                            ccAGTGATTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGc                          >  1:658565/1‑71 (MQ=255)
                                                 gtTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGtt                     >  1:1687200/1‑71 (MQ=255)
                                                  tttGCTCGGCTCGCAGGAACGGGCTCACCAGAACTCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGt                      <  1:294322/69‑1 (MQ=255)
                                                  tttGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:5077910/40‑1 (MQ=255)
                                                  tttGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                                   <  1:150462/40‑1 (MQ=255)
                                                  tttGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGt                      <  1:1714229/69‑1 (MQ=255)
                                                    tGCTCGGCTCGCAGGAACGGGCTCCCCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGc                  >  1:2574617/1‑71 (MQ=255)
                                                    tGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACttttt                                     >  1:1044952/1‑52 (MQ=255)
                                                    tGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACttttt                                     >  1:91130/1‑52 (MQ=255)
                                                     gctcggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGc                  <  1:1329056/70‑1 (MQ=255)
                                                      ctcggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGcc                 <  1:355641/70‑1 (MQ=255)
                                                      ctcggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGcc                 <  1:2548618/70‑1 (MQ=255)
                                                      ctcggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCa                <  1:2217549/71‑1 (MQ=255)
                                                       tcggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCAt               >  1:2977192/1‑71 (MQ=255)
                                                        cggctcgCAGGAACGGGCTCACCTGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGcctt               >  1:1584289/1‑68 (MQ=255)
                                                         ggctcgCAGGAACGGGCTCCCCAGAACGCGTTCGATTTCCACTTTTTGa                                   >  1:2639388/1‑49 (MQ=255)
                                                         ggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTc                    >  1:4388284/1‑64 (MQ=255)
                                                          gctcgCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAg            <  1:1192828/71‑1 (MQ=255)
                                                          gctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAg            >  1:3035706/1‑71 (MQ=255)
                                                           ctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAgg           >  1:1050222/1‑71 (MQ=255)
                                                              gCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCga        <  1:4437549/71‑1 (MQ=255)
                                                              gCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCga        <  1:4501313/71‑1 (MQ=255)
                                                                aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATc              >  1:2701899/1‑63 (MQ=255)
                                                                aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATc              >  1:2943856/1‑63 (MQ=255)
                                                                aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATc              >  1:5104278/1‑63 (MQ=255)
                                                                aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCgaga      <  1:1329409/71‑1 (MQ=255)
                                                                 ggAACGGGCTCACCAGAACTCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     <  1:1202632/71‑1 (MQ=255)
                                                                 ggAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     >  1:4675196/1‑71 (MQ=255)
                                                                 ggAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     <  1:4977131/71‑1 (MQ=255)
                                                                 ggAACGGGCTCACAAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     >  1:2760003/1‑71 (MQ=255)
                                                                  gAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     <  1:4047496/70‑1 (MQ=255)
                                                                  gAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGAt     <  1:3430042/70‑1 (MQ=255)
                                                                   aaCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTc   <  1:1466486/71‑1 (MQ=255)
                                                                   aaCGGGCGCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTc   <  1:3104991/71‑1 (MQ=255)
                                                                    aCGGGCTCACCAGAACGCGTTCGATTTCCACttttt                                     >  1:3185888/1‑36 (MQ=255)
                                                                    aCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTCg  <  1:818436/71‑1 (MQ=255)
                                                                    aCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTCg  <  1:2699372/71‑1 (MQ=255)
                                                                    |                                                                      
AACTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTCG  >  W3110S.gb/2462034‑2462172

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: