Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,111,845 A→C 24.1% L702R (CTG→CGG)  yrfF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,111,8450AC24.1% 63.4 / 44.5 54L702R (CTG→CGG) yrfFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/13);  new base C (7/6);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

CAGGGTAGCATAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGT  >  W3110S.gb/4111784‑4111911
                                                             |                                                                  
cAGGGTAGCATAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGt                                                           <  1:731276/71‑1 (MQ=255)
         aTAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGc                                                  >  1:666629/1‑71 (MQ=255)
         aTAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGc                                                  >  1:822578/1‑71 (MQ=255)
         aTAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGc                                                  >  1:1511155/1‑71 (MQ=255)
         aTAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGc                                                  >  1:1428020/1‑71 (MQ=255)
          tAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCt                                                 >  1:1074044/1‑71 (MQ=255)
          tAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCt                                                 >  1:409005/1‑71 (MQ=255)
          tAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCt                                                 >  1:1661219/1‑71 (MQ=255)
          tAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCt                                                 >  1:1192130/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCttt                                               >  1:271740/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCttt                                               >  1:259561/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCttt                                               >  1:750842/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCt                                                 >  1:1317608/1‑69 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1377353/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:391287/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:315084/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:269860/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:25307/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1054173/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:205106/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1779291/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:765350/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:905684/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:874582/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:847180/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1603925/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1460382/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1441007/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1286309/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCttt                                               >  1:1412257/1‑71 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCtt                                                >  1:790063/1‑70 (MQ=255)
            aCCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCtt                                                >  1:457540/1‑70 (MQ=255)
                 ccGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCTTTCat                                            <  1:79944/69‑1 (MQ=255)
                 ccGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCTTTCat                                            <  1:278972/69‑1 (MQ=255)
                 ccGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCata                                           <  1:1730377/70‑1 (MQ=255)
                 ccGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCata                                           <  1:1237592/70‑1 (MQ=255)
                  cgcgCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCTTTCat                                            <  1:951102/68‑1 (MQ=255)
                        aaCGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCt                                   >  1:1118795/1‑71 (MQ=255)
                          cGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGa                                 <  1:1724813/71‑1 (MQ=255)
                          cGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGa                                 <  1:523939/71‑1 (MQ=255)
                             gtTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGACAGCTTTCata                                           <  1:1704333/58‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:1757321/71‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:1738205/71‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:493028/71‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:859408/71‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:1215779/71‑1 (MQ=255)
                                     aTAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCcatcat                      <  1:1213262/71‑1 (MQ=255)
                                        aGGCTTTCTGAAGGGGTGATCCGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGa                                 <  1:890777/57‑1 (MQ=255)
                                                  aaGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGAt           <  1:1743192/69‑1 (MQ=255)
                                                  aaGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGAt           <  1:1717109/69‑1 (MQ=255)
                                                          aTCCGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGt  >  1:574318/1‑70 (MQ=255)
                                                          aTCCGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGt  >  1:706849/1‑70 (MQ=255)
                                                          aTCCGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGt  >  1:1648973/1‑70 (MQ=255)
                                                          aTCCGTTGCGGGTTCAGACAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGt  <  1:236814/70‑1 (MQ=255)
                                                             |                                                                  
CAGGGTAGCATAACCTGCCGCGCAAACGTGTTATTCGATAAGGCTTTCTGAAGGGGTGATCAGTTGCGGGTTCAGGCAGCTTTCATAATAGGCCTGAATCTCCATCATGCGAGTGCGATGACGCTGGT  >  W3110S.gb/4111784‑4111911

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: