Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,863,108 T→C 19.4% L73L (TTA→CTA)  ygbM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,863,1080TC19.4% 61.8 / 11.5 31L73L (TTA→CTA) ygbMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/15);  new base C (3/3);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGAATATGCAT  >  W3110S.gb/2863057‑2863174
                                                    |                                                                  
cTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGACTATCCGCCCTTCCCGGAc                                                  >  1:1736280/1‑70 (MQ=255)
  ggCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGGAATGGGGATTATCCGCCCTTCCCGGACGt                                                >  1:664935/1‑71 (MQ=255)
         ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa                                        <  1:1295417/71‑1 (MQ=255)
         ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa                                        <  1:2454021/71‑1 (MQ=255)
           cAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGaa                                       >  1:1275266/1‑70 (MQ=255)
            aaTACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg                                      <  1:474259/70‑1 (MQ=255)
             aTACCGCTCCCGGAGATATTAATGCGGGTGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGca                                    <  1:1371407/71‑1 (MQ=255)
             aTACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGca                                    <  1:1708873/71‑1 (MQ=255)
             aTACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGca                                    <  1:2372522/71‑1 (MQ=255)
               aCCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca                                  >  1:991009/1‑71 (MQ=255)
               aCCGCTCCCGGAGATATTAATGCGGGGGAATGGGGACTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca                                  >  1:2502862/1‑71 (MQ=255)
                   cTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCa                              <  1:2380494/71‑1 (MQ=255)
                          agaTATTAATGCGGGGGAATGGGGACTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCg                       <  1:1935604/71‑1 (MQ=255)
                           gaTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGa                      <  1:496105/71‑1 (MQ=255)
                              atTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGa                            >  1:1833363/1‑62 (MQ=255)
                               ttAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCa                                           >  1:1259542/1‑46 (MQ=255)
                               ttAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg                                <  1:923429/57‑1 (MQ=255)
                               ttAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGAc                           <  1:1489219/62‑1 (MQ=255)
                                     cGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTg            <  1:1318606/71‑1 (MQ=255)
                                     cGGGGGAATGGGGACTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca                                  >  1:2483491/1‑49 (MQ=255)
                                         gggAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGaa         <  1:466800/71‑1 (MQ=255)
                                         gggAATGGGGACTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGaa         <  1:1625130/71‑1 (MQ=255)
                                         gggAATGGGGACTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGaa         <  1:2599910/71‑1 (MQ=255)
                                              tGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa                                        <  1:828836/35‑1 (MQ=255)
                                              tGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg                                      <  1:1423859/37‑1 (MQ=255)
                                              tGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCa                              <  1:1750495/45‑1 (MQ=255)
                                               ggggATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGaa         >  1:1774809/1‑65 (MQ=255)
                                               ggggATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGaa         >  1:2549630/1‑65 (MQ=255)
                                                gggATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca                                  >  1:1810540/1‑39 (MQ=255)
                                                gggATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca                                  >  1:1686810/1‑39 (MQ=255)
                                                gggATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGAATATGCAt  >  1:1076532/1‑71 (MQ=255)
                                                    |                                                                  
CTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGACCTGGCACTGGAATATGCAT  >  W3110S.gb/2863057‑2863174

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: