Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,217,050 A→G 11.1% A75A (GCA→GCG yjaA → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,217,0500AG11.1% 155.6 / 10.5 63A75A (GCA→GCGyjaAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/19);  new base G (2/5);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 9.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

CCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTG  >  W3110S.gb/4216987‑4217114
                                                               |                                                                
ccacccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGGGCACTGGaaa                                                           >  1:1518844/1‑71 (MQ=255)
ccacccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaaa                                                           >  1:1032229/1‑71 (MQ=255)
ccacccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaa                                                            >  1:2309265/1‑70 (MQ=255)
  acccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatg                                                         <  1:107489/71‑1 (MQ=255)
   cccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatga                                                        <  1:3115539/71‑1 (MQ=255)
   cccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatga                                                        <  1:2066797/71‑1 (MQ=255)
    ccacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaa                                                       <  1:1537737/71‑1 (MQ=255)
     cacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCa                                                                  >  1:2896727/1‑59 (MQ=255)
     cacGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaat                                                      >  1:1908281/1‑71 (MQ=255)
      acGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaaa                                                           >  1:760468/1‑65 (MQ=255)
      acGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaaa                                                           >  1:2975031/1‑65 (MQ=255)
             aCCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGa                                                             <  1:2279372/56‑1 (MQ=255)
                tGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaaa                                                           >  1:1542589/1‑55 (MQ=255)
                tGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaatg                                                     <  1:3789326/61‑1 (MQ=255)
                tGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaa                                                       >  1:2450974/1‑59 (MQ=255)
                 ggCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaa                                                       >  1:2484536/1‑58 (MQ=255)
                  gCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaat                                                      >  1:2655252/1‑58 (MQ=255)
                  gCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGt                                         <  1:2631288/71‑1 (MQ=255)
                    aTTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGttt                                           <  1:3503385/67‑1 (MQ=255)
                     ttCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaa                                                            <  1:1568556/49‑1 (MQ=255)
                     ttCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGaa                                                            <  1:1161490/49‑1 (MQ=255)
                     ttCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaat                                                      <  1:1838401/55‑1 (MQ=255)
                     ttCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTca                                       >  1:3150493/1‑70 (MQ=255)
                     ttCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTc                                        >  1:1110452/1‑69 (MQ=255)
                      tCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGt                                         <  1:3112942/67‑1 (MQ=255)
                        ttttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTcaca                                     >  1:665957/1‑69 (MQ=255)
                        ttttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTcaca                                     >  1:1977235/1‑69 (MQ=255)
                        ttttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTcaca                                     >  1:656640/1‑69 (MQ=255)
                        ttttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAgg                                   >  1:2931679/1‑71 (MQ=255)
                        ttttttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAgg                                   >  1:160897/1‑71 (MQ=255)
                           tttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGtt                                 >  1:863669/1‑70 (MQ=255)
                           tttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGtt                                 >  1:856998/1‑70 (MQ=255)
                           tttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGtt                                 >  1:2336960/1‑70 (MQ=255)
                           tttGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGtt                                 >  1:224665/1‑70 (MQ=255)
                             tGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGa                               <  1:2365978/70‑1 (MQ=255)
                              gCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAg                              <  1:1089286/70‑1 (MQ=255)
                               ccAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAgg                             <  1:1151174/70‑1 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACg                         >  1:3365963/1‑71 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACg                         >  1:1330018/1‑71 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaatgaa                                                       >  1:2264166/1‑41 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaat                                                          <  1:3712361/38‑1 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCACTGGAaat                                                          <  1:1061632/38‑1 (MQ=255)
                                  gCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACg                         <  1:487468/71‑1 (MQ=255)
                                   cAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACg                         >  1:2393539/1‑70 (MQ=255)
                                    aaaaCGTATGGATATTGTTGTGAGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCa                       <  1:1071180/71‑1 (MQ=255)
                                    aaaaCGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGt                                         >  1:439002/1‑53 (MQ=255)
                                     aaaCGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAg                                                   >  1:2647741/1‑42 (MQ=255)
                                       aCGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATtc                    >  1:984507/1‑71 (MQ=255)
                                         gTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATtctc                  >  1:1697872/1‑71 (MQ=255)
                                              gATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCAtt                     >  1:457999/1‑63 (MQ=255)
                                              gATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATg              >  1:3425991/1‑70 (MQ=255)
                                              gATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATg              >  1:2289786/1‑70 (MQ=255)
                                                 atTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGtt                                 <  1:291468/48‑1 (MQ=255)
                                                 atTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTGTCCATGaa            >  1:1398015/1‑69 (MQ=255)
                                                   tgttgtGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGTACGCAt                      >  1:2595232/1‑57 (MQ=255)
                                                   tgttgtGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTgg        >  1:961702/1‑71 (MQ=255)
                                                    gttgtGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGa                               >  1:2670213/1‑47 (MQ=255)
                                                    gttgtGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGa                               >  1:1459744/1‑47 (MQ=255)
                                                    gttgtGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGaa                           >  1:466015/1‑51 (MQ=255)
                                                       gtgAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGc    >  1:2811549/1‑71 (MQ=255)
                                                          aGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTTCATGAAGTGGTCGcag  <  1:2582986/70‑3 (MQ=255)
                                                          aGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTTCATGAAGTGGTCGcag  <  1:1790848/70‑3 (MQ=255)
                                                          aGCGCGCTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTTCATGAAGTGGTCGcag  <  1:1717435/70‑3 (MQ=255)
                                                          aGCGCGCTGGAAATGAATGAGGGAGGTTTGTCACAGGTTGAGGAACGCATTCTTCATGAAGTGGTCGcag  <  1:2508033/70‑3 (MQ=255)
                                                               |                                                                
CCACCCACGTTCAACCTGGCATTCTTTTTTGCCAGCAAAACGTATGGATATTGTTGTGAGCGCACTGGAAATGAATGAGGGCGGTTTGTCACAGGTTGAGGAACGCATTCTCCATGAAGTGGTCGCTG  >  W3110S.gb/4216987‑4217114

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: