Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,591,880 C→T 56.9% 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,8800CT56.9% 19.1 / 65.4 51M201I (ATG→ATAyihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base T (13/16);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTAC  >  W3110S.gb/3591815‑3591947
                                                                 |                                                                   
ggAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTgg                                                                >  1:1071239/1‑71 (MQ=255)
ggAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTgg                                                                >  1:562667/1‑71 (MQ=255)
ggAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTgg                                                                >  1:425508/1‑71 (MQ=255)
ggAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTgg                                                                >  1:1362259/1‑71 (MQ=255)
 gAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGt                                                               <  1:1168856/71‑1 (MQ=255)
 gAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGt                                                               <  1:34870/71‑1 (MQ=255)
                ctAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                 <  1:1154817/70‑1 (MQ=255)
                ctAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAAc                                                 <  1:1307404/70‑1 (MQ=255)
                   aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAat                                              <  1:1155019/70‑1 (MQ=255)
                   aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAat                                              <  1:248266/70‑1 (MQ=255)
                   aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAat                                              <  1:240732/70‑1 (MQ=255)
                   aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAat                                              <  1:1904262/70‑1 (MQ=255)
                   aTCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAat                                              <  1:1467849/70‑1 (MQ=255)
                        gcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTc                                        <  1:981434/71‑1 (MQ=255)
                        gcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTc                                        <  1:56786/71‑1 (MQ=255)
                        gcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTc                                        <  1:1292494/71‑1 (MQ=255)
                        gcCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTc                                        <  1:1295019/71‑1 (MQ=255)
                           agagAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGAc                                     >  1:564976/1‑71 (MQ=255)
                           agagAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGAc                                     <  1:1213425/71‑1 (MQ=255)
                           agagAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGAc                                     <  1:1526621/71‑1 (MQ=255)
                              gaATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACttt                                  >  1:1070481/1‑71 (MQ=255)
                              gaATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACttt                                  >  1:1845677/1‑71 (MQ=255)
                              gaATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACttt                                  >  1:1750200/1‑71 (MQ=255)
                              gaATCAATGGTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACttt                                  >  1:644209/1‑71 (MQ=255)
                                  cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCt                              >  1:1608494/1‑71 (MQ=255)
                                  cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCt                              >  1:177987/1‑71 (MQ=255)
                                  cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCt                              >  1:1798420/1‑71 (MQ=255)
                                  cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCt                              >  1:1449848/1‑71 (MQ=255)
                                  cAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCt                              >  1:1465102/1‑71 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:563905/71‑1 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:167295/71‑1 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:1706518/71‑1 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:1969407/71‑1 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:192332/71‑1 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCaa                          >  1:1459469/1‑70 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCaa                          >  1:1461956/1‑70 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCaa                          >  1:109879/1‑70 (MQ=255)
                                       gTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATATTCGACTTTGCCTTCa                           >  1:598248/1‑69 (MQ=255)
                                        tAAAACTTTTTGCCGTTACTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:1082944/70‑1 (MQ=255)
                                                    ccGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGAc            <  1:614202/71‑1 (MQ=255)
                                                    ccGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGAc            <  1:1575447/71‑1 (MQ=255)
                                                         aCTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg        >  1:1437099/1‑70 (MQ=255)
                                                         aCTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg        >  1:1182236/1‑70 (MQ=255)
                                                         aCTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg        >  1:1899146/1‑70 (MQ=255)
                                                         aCTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg        >  1:1343337/1‑70 (MQ=255)
                                                         aCTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATg        >  1:1442433/1‑70 (MQ=255)
                                                          cTTTGTTTATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGa       <  1:1294247/70‑1 (MQ=255)
                                                          cTTTGTTTATTGCTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGa       <  1:426364/70‑1 (MQ=255)
                                                             tGTTTATTGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAat                         <  1:204485/48‑1 (MQ=255)
                                                               ttCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTAc  <  1:1207468/70‑1 (MQ=255)
                                                               ttCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTAc  <  1:932124/70‑1 (MQ=255)
                                                                 |                                                                   
GGAAATTTCTTTTCTTCTAATCGCGCCAGAGAATCAATGGTAAAACTTTTTGCCGTTACTTTGTTCATTGGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTAC  >  W3110S.gb/3591815‑3591947

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: