Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,489,179 T→C 26.1% F77F (TTT→TTC ydcF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,489,1790TC26.1% 108.6 / 44.3 65F77F (TTT→TTCydcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/18);  new base C (6/11);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 8.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAC  >  W3110S.gb/1489111‑1489247
                                                                    |                                                                    
tGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACtttctt                                                                    <  1:3431195/71‑1 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:2689795/70‑1 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:1143130/1‑70 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:2949491/70‑1 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:652602/1‑70 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:1658524/1‑70 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:1687337/1‑70 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:3645999/1‑70 (MQ=255)
  tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTg                                                                   <  1:3478154/70‑1 (MQ=255)
   aaGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:232408/69‑1 (MQ=255)
    aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 <  1:2427592/70‑1 (MQ=255)
    aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 >  1:2131250/1‑70 (MQ=255)
    aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 <  1:3583887/70‑1 (MQ=255)
      aTTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:1865964/1‑68 (MQ=255)
        tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 >  1:805062/1‑66 (MQ=255)
        tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 >  1:1242079/1‑66 (MQ=255)
             gcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATc                                                       >  1:2138132/1‑71 (MQ=255)
             gcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATc                                                       >  1:1253376/1‑71 (MQ=255)
             gcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATc                                                       >  1:2950373/1‑71 (MQ=255)
             gcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:596492/1‑61 (MQ=255)
              cgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCg                                                      <  1:2038703/71‑1 (MQ=255)
              cgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCCCCATCg                                                      <  1:685586/71‑1 (MQ=255)
                aTCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTg                                                                   <  1:3342378/56‑1 (MQ=255)
                  cAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 >  1:3472826/1‑56 (MQ=255)
                  cAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 >  1:973918/1‑56 (MQ=255)
                   aGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAgcgc                                                           <  1:194338/61‑1 (MQ=255)
                     cAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCa                                               >  1:3101191/1‑71 (MQ=255)
                         ttCCTTTACTGATTAGTGGTGATATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCg                                           >  1:2090606/1‑71 (MQ=255)
                                      tAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:2890957/1‑36 (MQ=255)
                                      tAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:2393120/1‑36 (MQ=255)
                                       aGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTa                              >  1:774910/1‑70 (MQ=255)
                                           gtggtATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGc                         <  1:347287/71‑1 (MQ=255)
                                           gtggtATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGc                         <  1:3542045/71‑1 (MQ=255)
                                                 tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:250482/64‑1 (MQ=255)
                                                 tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:3953262/64‑1 (MQ=255)
                                                 tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:1687575/64‑1 (MQ=255)
                                                 tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:1595403/64‑1 (MQ=255)
                                                 tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACt                   >  1:1866500/1‑71 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                             >  1:3757362/1‑43 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                             >  1:3222120/1‑43 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                             >  1:3253978/1‑43 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:3877931/1‑71 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 <  1:3556830/71‑1 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:2485231/1‑71 (MQ=255)
                                                   ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:3927350/1‑71 (MQ=255)
                                                      cACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTac                                     <  1:3256053/48‑1 (MQ=255)
                                                       aCTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACt                               <  1:3140485/53‑1 (MQ=255)
                                                       aCTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcag              <  1:3211297/70‑1 (MQ=255)
                                                       aCTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTAt                             >  1:3314805/1‑55 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACt                                       <  1:3695178/44‑1 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACt                                       <  1:2779186/44‑1 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACt                                       <  1:3205438/44‑1 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaa                                   >  1:1324102/1‑48 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaa                                   >  1:4074206/1‑48 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaa                                   >  1:3043874/1‑48 (MQ=255)
                                                        cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaa                                   >  1:312957/1‑48 (MQ=255)
                                                          cGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagag            >  1:3063864/1‑69 (MQ=255)
                                                           gACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagag            >  1:1116668/1‑68 (MQ=255)
                                                           gACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagag            >  1:1706078/1‑68 (MQ=255)
                                                           gACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagag            >  1:3068936/1‑68 (MQ=255)
                                                            aCAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACt                   <  1:3079251/60‑1 (MQ=255)
                                                            aCAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcag         <  1:465958/70‑1 (MQ=255)
                                                            aCAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcag         <  1:2233719/70‑1 (MQ=255)
                                                            aCAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcag         <  1:1353022/70‑1 (MQ=255)
                                                            aCAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcag         <  1:1114893/70‑1 (MQ=255)
                                                              aaCTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcaga        <  1:1851413/69‑1 (MQ=255)
                                                                  tttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  >  1:3517541/1‑71 (MQ=255)
                                                                  tttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  >  1:425561/1‑71 (MQ=255)
                                                                   ttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  <  1:550100/70‑1 (MQ=255)
                                                                   ttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  <  1:3340535/70‑1 (MQ=255)
                                                                    |                                                                    
TGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAC  >  W3110S.gb/1489111‑1489247

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