Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,088,483 C→T 73.9% Q599Q (CAG→CAA ycdR ← predicted enzyme associated with biofilm formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,088,4830CT73.9% 34.0 / 12.7 23Q599Q (CAG→CAAycdRpredicted enzyme associated with biofilm formation
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (10/7);  total (14/9)
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.70e-01

GCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGC  >  W3110S.gb/1088418‑1088545
                                                                 |                                                              
gCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAgg                                                           >  1:1147548/1‑71 (MQ=255)
   tGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGAt                                                        >  1:1341743/1‑71 (MQ=255)
    gATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATg                                                       >  1:20400/1‑71 (MQ=255)
         tGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGtttt                                                   <  1:757231/70‑1 (MQ=255)
              aTTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTg                                                              <  1:32685/54‑1 (MQ=255)
               tttttCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAgg                                                           <  1:1437014/56‑1 (MQ=255)
                 tttCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGAt                                          >  1:1677853/1‑71 (MQ=255)
                       ccAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTc                                     <  1:1732802/70‑1 (MQ=255)
                       ccAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTc                                     <  1:284713/70‑1 (MQ=255)
                        cAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGAt                                                        >  1:934596/1‑50 (MQ=255)
                        cAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGAt                                                        >  1:284746/1‑50 (MQ=255)
                                      tAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAAttt                                             >  1:829800/1‑47 (MQ=255)
                                      tAATTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTg                                        >  1:470308/1‑52 (MQ=255)
                                       aaTTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTg                    >  1:131255/1‑71 (MQ=255)
                                       aaTTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTg                    >  1:1212029/1‑71 (MQ=255)
                                       aaTTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTg                    >  1:458353/1‑71 (MQ=255)
                                       aaTTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTg                    >  1:677339/1‑71 (MQ=255)
                                       aaTTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTg                    >  1:1084682/1‑71 (MQ=255)
                                        aTTCTAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTgg                   >  1:112532/1‑71 (MQ=255)
                                         ttCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGt                  <  1:192368/71‑1 (MQ=255)
                                            tAAAATAGATTTGTCTTTAGCTTGAGGGATGTTTTTAAttt                                             <  1:478157/41‑1 (MQ=255)
                                              aaaTAGATTTGTCTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGcc             <  1:17449/71‑1 (MQ=255)
                                                          cTTTAGCTTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGc  <  1:783448/70‑1 (MQ=255)
                                                                 |                                                              
GCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGC  >  W3110S.gb/1088418‑1088545

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: