Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,910,338 C→T 8.3% *48* (TGA→TAA)  yobG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,910,3380CT 8.3% 145.8 / 2.9 48*48* (TGA→TAA) yobGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base T (2/2);  total (27/21)
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.85e-01

CACCTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCT  >  W3110S.gb/1910268‑1910392
                                                                      |                                                      
caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:2646661/71‑1 (MQ=255)
caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:1516310/71‑1 (MQ=255)
caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:1944764/71‑1 (MQ=255)
     ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:69839/1‑71 (MQ=255)
     ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:684058/1‑71 (MQ=255)
     ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:578831/1‑71 (MQ=255)
                cgcTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTg                                        >  1:733668/1‑71 (MQ=255)
                 gcTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTgg                                       >  1:644541/1‑71 (MQ=255)
                    aTTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGtt                                     >  1:728226/1‑70 (MQ=255)
                     ttCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcac                                                      <  1:2105281/52‑1 (MQ=255)
                     ttCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:1570573/1‑71 (MQ=255)
                       cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:2997901/1‑71 (MQ=255)
                       cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:2564313/1‑71 (MQ=255)
                       cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:1757896/1‑71 (MQ=255)
                            aCTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGcaca                            <  1:737596/71‑1 (MQ=255)
                              tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:46304/1‑71 (MQ=255)
                              tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:4059665/1‑71 (MQ=255)
                              tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:2579260/1‑71 (MQ=255)
                              agctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATCGCACaaa                          >  1:2583904/2‑71 (MQ=255)
                               gctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAAt                         >  1:5239/1‑71 (MQ=255)
                                ctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATaaa                                           >  1:347112/1‑52 (MQ=255)
                                ctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATaaa                                           >  1:948052/1‑52 (MQ=255)
                                          agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:925949/1‑49 (MQ=255)
                                          agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:115167/1‑49 (MQ=255)
                                          agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:2188871/1‑49 (MQ=255)
                                          agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa              <  1:859078/71‑1 (MQ=255)
                                           gaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATa                                             <  1:3030000/39‑1 (MQ=255)
                                            aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa              <  1:3758526/69‑1 (MQ=255)
                                             aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATa                                             <  1:1616271/37‑1 (MQ=255)
                                             aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:3953363/1‑47 (MQ=255)
                                             aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:3108795/1‑47 (MQ=255)
                                              aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaa               <  1:2713563/66‑1 (MQ=255)
                                              aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:3706595/71‑1 (MQ=255)
                                              aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:2976381/71‑1 (MQ=255)
                                              aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          >  1:779246/1‑71 (MQ=255)
                                               aTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:1564016/70‑1 (MQ=255)
                                                  tAGTGCTGAAAAAATGATATCACCACGGGATAAACTgg                                       >  1:3388313/1‑38 (MQ=255)
                                                    gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:1491659/1‑40 (MQ=255)
                                                    gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       <  1:2083158/68‑1 (MQ=255)
                                                     tGCTGAAAAAATGATATTACCACGGGATAAACTGGtt                                     <  1:97111/37‑1 (MQ=255)
                                                     tGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTa        <  1:2520392/66‑1 (MQ=255)
                                                       cTGAAAAAATGATATTACCACGGGATAAACTGGTTaa                                   >  1:3187843/1‑37 (MQ=255)
                                                       cTGAAAAAATGATATTACCACGGGATAAACTGGTTaa                                   >  1:125420/1‑37 (MQ=255)
                                                       cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  <  1:1896717/70‑1 (MQ=255)
                                                       cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  >  1:577635/1‑70 (MQ=255)
                                                       cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  <  1:746556/70‑1 (MQ=255)
                                                               aTGATATTACCACGGGATAAACTGGTTAATGGCACaaa                          <  1:2404575/38‑1 (MQ=255)
                                                               aTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       >  1:1164366/1‑57 (MQ=255)
                                                               aTGATAGCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       >  1:522408/1‑57 (MQ=255)
                                                                 gATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa             <  1:3152796/49‑1 (MQ=255)
                                                                 gATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa             <  1:2241729/49‑1 (MQ=255)
                                                                     tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:607262/48‑1 (MQ=255)
                                                                     tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:882785/48‑1 (MQ=255)
                                                                     tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:887647/48‑1 (MQ=255)
                                                                     tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:2681780/48‑1 (MQ=255)
                                                                      |                                                      
CACCTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCT  >  W3110S.gb/1910268‑1910392

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: