Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,591,880 C→T 10.4% M201I (ATG→ATA yihF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,591,8800CT10.4% 194.5 / 10.8 67M201I (ATG→ATAyihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base T (2/5);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGGAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCA  >  W3110S.gb/3591810‑3591949
                                                                      |                                                                     
cAACCGGAAATTTCTTTTCTTATAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTc                                                                       <  1:1774128/71‑1 (MQ=255)
     ggAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTgg                                                                  >  1:616021/1‑71 (MQ=255)
        aaTTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGa                                                               >  1:757371/1‑71 (MQ=255)
         aTTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGAt                                                              <  1:1338173/71‑1 (MQ=255)
         aTTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGAt                                                              <  1:2540154/71‑1 (MQ=255)
          tttcttttctTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATa                                                             <  1:1829215/71‑1 (MQ=255)
           ttcttttcttCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCa                                                                      >  1:2375024/1‑61 (MQ=255)
           ttcttttcttCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATaa                                                            >  1:29476/1‑71 (MQ=255)
           ttcttttcttCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCAGTTACTTTGTTCATTGGTGATaa                                                            >  1:1179380/1‑71 (MQ=255)
                ttcttcTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGc                                                        <  1:1072121/70‑1 (MQ=255)
                        aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACaa                                                 <  1:1958038/69‑1 (MQ=255)
                        aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAat                                                <  1:612719/70‑1 (MQ=255)
                        aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAat                                                <  1:1671155/70‑1 (MQ=255)
                        aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAat                                                <  1:737990/70‑1 (MQ=255)
                        aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAat                                                <  1:234640/70‑1 (MQ=255)
                           gcgcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCt                                            >  1:2567156/1‑71 (MQ=255)
                           gcgcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCt                                            >  1:130605/1‑71 (MQ=255)
                                 gagaATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCa                                                                      >  1:516453/1‑39 (MQ=255)
                                    aaTCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATaa                                                            >  1:1067783/1‑46 (MQ=255)
                                    aaTCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTg                                   >  1:2311968/1‑71 (MQ=255)
                                      tCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAatat                                              >  1:919687/1‑58 (MQ=255)
                                      tCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGcc                                 >  1:2418917/1‑71 (MQ=255)
                                       cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAatat                                              >  1:2114984/1‑57 (MQ=255)
                                        aaTGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAata                                               >  1:2388593/1‑55 (MQ=255)
                                        aaTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACaa                                                 >  1:1598810/1‑53 (MQ=255)
                                         aTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              <  1:985537/71‑1 (MQ=255)
                                         aTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              <  1:596547/71‑1 (MQ=255)
                                         aTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              <  1:1474393/71‑1 (MQ=255)
                                         aTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              <  1:1662184/71‑1 (MQ=255)
                                         aTGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATATTCGACTTTGCCTTc                              <  1:2209825/71‑1 (MQ=255)
                                          tGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGAt                                                              <  1:1354331/38‑1 (MQ=255)
                                          tGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                   <  1:105897/49‑1 (MQ=255)
                                          tGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                   <  1:884425/49‑1 (MQ=255)
                                          tGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                   <  1:644850/49‑1 (MQ=255)
                                          tGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              >  1:2808021/1‑70 (MQ=255)
                                            gTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCa                                                       <  1:1778215/43‑1 (MQ=255)
                                             tAAAACTTTTTGCCGTTACTTTGTTCATTGGTGAt                                                              <  1:1850916/35‑1 (MQ=255)
                                             tAAAACTTTTTGCCGTTACTTTGTTCATTGGTGAt                                                              <  1:1227587/35‑1 (MQ=255)
                                             tAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                   <  1:420994/46‑1 (MQ=255)
                                             tAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                   <  1:2639478/46‑1 (MQ=255)
                                              aaaaCTTTTTGCCGTTACTTTGTTCATTGGTGATaa                                                            >  1:1194277/1‑36 (MQ=255)
                                              aaaaCTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGa                                        >  1:1175457/1‑56 (MQ=255)
                                              aaaaCTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAatat                         >  1:468628/1‑71 (MQ=255)
                                                 aCTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATc                        >  1:927207/1‑69 (MQ=255)
                                                  cTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAg                     >  1:560817/1‑71 (MQ=255)
                                                   tttttGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGAc                                       <  1:284076/52‑1 (MQ=255)
                                                     tttGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGt                   <  1:2688908/70‑1 (MQ=255)
                                                     tttGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTCCCTTCAATATCCAGGt                   <  1:623154/70‑1 (MQ=255)
                                                         ccGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGAc              <  1:1835127/71‑1 (MQ=255)
                                                            ttACTTTGTTTATTGGTGATAATTGCAGAACAATATc                                             <  1:355587/37‑1 (MQ=255)
                                                            ttACTTTGTTTATTGGTGATAATTGCAGAACAATATc                                             <  1:2696655/37‑1 (MQ=255)
                                                            ttACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGcc                                 <  1:763472/49‑1 (MQ=255)
                                                            ttACTTTGTTTATTGGGGATAATTGCAGAACAATATCTTCGACTTTGcc                                 <  1:2222597/49‑1 (MQ=255)
                                                            ttACTTTGTTCATTGGTGATAATTGCAGAACAATATc                                             <  1:2551974/37‑1 (MQ=255)
                                                            ttACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCa                             <  1:2240689/53‑1 (MQ=255)
                                                              aCTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg          <  1:589148/70‑1 (MQ=255)
                                                              aCTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg          <  1:1432486/70‑1 (MQ=255)
                                                               cTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg          <  1:214322/69‑1 (MQ=255)
                                                               cTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGa         <  1:1752927/70‑1 (MQ=255)
                                                                  tGTTCATTGGTGATAATTGCAGAACAATATCTTCGACt                                      >  1:2236388/1‑38 (MQ=255)
                                                                  tGTTCATTGGTGATAATTGCAGAACAATATCTTCGACt                                      >  1:1755799/1‑38 (MQ=255)
                                                                  tGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTc                              <  1:830577/46‑1 (MQ=255)
                                                                  tGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGAc              >  1:975266/1‑62 (MQ=255)
                                                                     tCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAgg                    >  1:398983/1‑53 (MQ=255)
                                                                     tCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAgg                    >  1:1917265/1‑53 (MQ=255)
                                                                     tCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAgg                    >  1:279414/1‑53 (MQ=255)
                                                                     tCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACa             >  1:2431071/1‑60 (MQ=255)
                                                                     tCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCa  >  1:2258847/1‑71 (MQ=255)
                                                                      |                                                                     
CAACCGGAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCA  >  W3110S.gb/3591810‑3591949

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: