Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,142,075 T→C 32.6% T559A (ACG→GCG)  asmA ← predicted assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,142,0750TC32.6% 48.5 / 36.9 43T559A (ACG→GCG) asmApredicted assembly protein
Reads supporting (aligned to +/- strand):  ref base T (9/20);  new base C (7/7);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTG  >  W3110S.gb/2142006‑2142141
                                                                     |                                                                  
gTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGtt                                                                   <  1:2936076/71‑1 (MQ=255)
   aCTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTAATGGCGcttct                                                                >  1:833820/1‑71 (MQ=255)
      gATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGtttcttt                                                              >  1:866069/1‑70 (MQ=255)
         ccACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAg                                                           >  1:2746850/1‑70 (MQ=255)
          cACTTGCAGGGTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAgg                                                          >  1:142095/1‑70 (MQ=255)
          cACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAgg                                                          <  1:860045/70‑1 (MQ=255)
          cACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGa                                                         <  1:800913/71‑1 (MQ=255)
          cACTTGCAGGCTGTAAATGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAgg                                                          <  1:473/70‑1 (MQ=255)
              tGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGt                                                                    >  1:1729774/1‑56 (MQ=255)
                 aGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAaatcaat                                                  <  1:3103740/71‑1 (MQ=255)
                 aGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAaatcaat                                                  <  1:262919/71‑1 (MQ=255)
                     tGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGt                                              >  1:3146406/1‑71 (MQ=255)
                          ttGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAatcaatca                                                <  1:3729123/64‑1 (MQ=255)
                          ttGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAatcaatca                                                <  1:575892/64‑1 (MQ=255)
                          ttGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAatcaatca                                                <  1:2906061/64‑1 (MQ=255)
                          ttGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAg                                               <  1:878648/65‑1 (MQ=255)
                           tGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAatcaatca                                                <  1:3481670/63‑1 (MQ=255)
                            gAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGtttct                                                                >  1:3801411/1‑46 (MQ=255)
                            gAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGtttct                                                                >  1:1067875/1‑46 (MQ=255)
                            gAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAg                                               >  1:807948/1‑63 (MQ=255)
                            gAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGt                                              >  1:443882/1‑64 (MQ=255)
                            gAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGt                                              >  1:3705147/1‑64 (MQ=255)
                                       aGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTc                               >  1:1775670/1‑68 (MQ=255)
                                         ttGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCtt                                        <  1:1680467/57‑1 (MQ=255)
                                         ttGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCtt                                        <  1:1540670/57‑1 (MQ=255)
                                          tGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGTAATCAATCAGTTTGCt                                         >  1:1532213/1‑55 (MQ=255)
                                          tGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCt                                         >  1:2480234/1‑55 (MQ=255)
                                          tGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGcca                         >  1:1624086/1‑71 (MQ=255)
                                          tGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGcca                         >  1:414019/1‑71 (MQ=255)
                                          tGCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGcca                         >  1:1539158/1‑71 (MQ=255)
                                           gCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAg                                               <  1:2309469/48‑1 (MQ=255)
                                           gCCATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGccac                        >  1:974628/1‑71 (MQ=255)
                                            ccATAAACCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGt                                 <  1:3179501/61‑1 (MQ=255)
                                                 aaCCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCtt                                        >  1:3736283/1‑49 (MQ=255)
                                                 aaCCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACgc                  >  1:3256743/1‑70 (MQ=255)
                                                 aaCCCGCAGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGa                  >  1:395392/1‑71 (MQ=255)
                                                     cGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCg              <  1:1120933/71‑1 (MQ=255)
                                                       cAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGaa            <  1:2438259/71‑1 (MQ=255)
                                                        aGCGGTACTGGCGCTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAat           <  1:1803193/71‑1 (MQ=255)
                                                            gTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCg              <  1:1969305/64‑1 (MQ=255)
                                                             tACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGc                                          <  1:1134710/35‑1 (MQ=255)
                                                             tACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCaa      <  1:3946787/71‑1 (MQ=255)
                                                               cTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCg              <  1:3185329/61‑1 (MQ=255)
                                                               cTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCg              <  1:832203/61‑1 (MQ=255)
                                                                  gCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGtt                                <  1:3653440/40‑1 (MQ=255)
                                                                  gCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGtt                                <  1:1938605/40‑1 (MQ=255)
                                                                  gCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTg  <  1:2795125/70‑1 (MQ=255)
                                                                  gCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTg  <  1:3424867/70‑1 (MQ=255)
                                                                  gCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTGCCAGCCACCCACGACCCGAATATCAAACt   <  1:2024987/69‑1 (MQ=255)
                                                                     tttctttcAGGAAATCAATCAGTTTGCTTTCCCCGtt                                <  1:1893139/37‑1 (MQ=255)
                                                                     |                                                                  
GTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTCCAGCCACCCACGACCCGAATATCAAACTG  >  W3110S.gb/2142006‑2142141

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: