Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,783,965 C→T 35.9% G13D (GGC→GAC)  yidI ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,783,9650CT35.9% 36.4 / 35.1 39G13D (GGC→GAC) yidIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base T (8/6);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CGGAGAAAGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAA  >  W3110S.gb/3783900‑3784027
                                                                 |                                                              
cGGAGAAAGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAAt                                                           >  1:276138/1‑71 (MQ=255)
cGGAGAAAGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAAt                                                           >  1:606179/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    >  1:687158/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:217014/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:433806/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:1637931/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    >  1:771393/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:1401980/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    >  1:920030/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    >  1:1137250/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    >  1:962341/1‑71 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCAAATGAAc                                                       >  1:1835473/1‑52 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCTATGAACTGATTTTATTTTTAGCTAt                                    <  1:1412836/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:1023449/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:581528/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:921261/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:277854/71‑1 (MQ=255)
                       tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:1701421/71‑1 (MQ=255)
                        ggATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:1002790/70‑1 (MQ=255)
                         gATAATCCCCATCATTAAAGCGATGGCACCGAAGAGCATGCCCAATGAA‑‑GATTTTATTTTGAGCTATGATa                                >  1:635407/1‑71 (MQ=255)
                         gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGa                                  >  1:962324/1‑71 (MQ=255)
                         gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGa                                  >  1:959458/1‑71 (MQ=255)
                         gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGa                                  >  1:471741/1‑71 (MQ=255)
                         gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGa                                  >  1:304878/1‑71 (MQ=255)
                         gATAATCCACATCATTAAAGCGATGTCACCGAACATCATGCCCAATGAACTGATTTTATTTTGAGatat                                    <  1:878813/69‑1 (MQ=255)
                           tAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:607724/67‑1 (MQ=255)
                              tCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGattt                                                 <  1:382127/51‑1 (MQ=255)
                              tCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGattt                                                 <  1:649816/51‑1 (MQ=255)
                              tCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGattt                                                 <  1:1840696/51‑1 (MQ=255)
                                cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGatttta                                               >  1:1873597/1‑51 (MQ=255)
                                  catcatTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGattt                                                 <  1:298165/47‑1 (MQ=255)
                                   atcatTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTg                        <  1:1253100/71‑1 (MQ=255)
                                      atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAcc                     >  1:1619816/1‑71 (MQ=255)
                                      atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAcc                     >  1:1464367/1‑71 (MQ=255)
                                      atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAcc                     >  1:336779/1‑71 (MQ=255)
                                      atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAc                      >  1:133555/1‑70 (MQ=255)
                                                          cAGCATGTCCAATGAACTGATTTTATTTTGAGCTAt                                    <  1:478095/36‑1 (MQ=255)
                                                          cAGCATGTCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCaa  >  1:1452568/1‑70 (MQ=255)
                                                          cAGCATGTCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCaa  >  1:81470/1‑70 (MQ=255)
                                                                 |                                                              
CGGAGAAAGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAA  >  W3110S.gb/3783900‑3784027

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: