Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,390,267 A→G 12.1% N65D (AAC→GAC)  yjeN → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,390,2670AG12.1% 81.2 / 4.9 33N65D (AAC→GAC) yjeNhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/9);  new base G (2/2);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TAAAGTGGCGATGGTGGTGGGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAATTAGTGG  >  W3110S.gb/4390200‑4390330
                                                                   |                                                               
tAAAGTGGCGATGGTGGTGGGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGaaa                                                                <  1:1382213/69‑1 (MQ=255)
      ggCGATGGTGGTGGGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGc                                                           >  1:3462539/1‑68 (MQ=255)
            ggtggtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGata                                                     >  1:326405/1‑68 (MQ=255)
            ggtggtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGata                                                     >  1:2962071/1‑68 (MQ=255)
             gtggtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtatt                                                   >  1:3349014/1‑69 (MQ=255)
             gtggtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtatt                                                   <  1:2973012/69‑1 (MQ=255)
             gtggtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtatt                                                   >  1:2471717/1‑69 (MQ=255)
                gtggGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtattt                                                  <  1:253486/67‑1 (MQ=255)
                      aGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAGACCTGCCGGATATTTATTGCGct                                          <  1:103969/69‑1 (MQ=255)
                      aGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtattt                                                  >  1:1115706/1‑61 (MQ=255)
                      aGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtattt                                                  >  1:3164606/1‑61 (MQ=255)
                      aGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtattt                                                  >  1:2025338/1‑61 (MQ=255)
                      aGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtattt                                                  >  1:1844247/1‑61 (MQ=255)
                       gCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGctt                                         <  1:3429854/69‑1 (MQ=255)
                       gCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGctt                                         <  1:2728298/69‑1 (MQ=255)
                        cACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGAtatt                                                   >  1:965699/1‑58 (MQ=255)
                         actaTTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGct                                          <  1:1135845/63‑1 (MQ=255)
                         aCCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttat                                       >  1:2131876/1‑69 (MQ=255)
                            aTTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttatctt                                    >  1:68854/1‑69 (MQ=255)
                               aaCACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttatcttat                                  >  1:3104327/1‑68 (MQ=255)
                               aaCACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttatcttat                                  >  1:1642523/1‑68 (MQ=255)
                               aaCACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttatcttat                                  >  1:581941/1‑68 (MQ=255)
                               aaCACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGctt                                         <  1:2568215/61‑1 (MQ=255)
                                          aGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGt                      >  1:1176580/1‑69 (MQ=255)
                                          aGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGt                      >  1:3519633/1‑69 (MQ=255)
                                             tCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttat                                       <  1:3128151/49‑1 (MQ=255)
                                             tCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGcttat                                       >  1:2203224/1‑49 (MQ=255)
                                             tCGCGATGGTAAACGGTGGGGAAACCTGACGGATATTTATTGCGcttat                                       >  1:693134/1‑49 (MQ=255)
                                              cgcgATGGTAAACGGTGGGGAGACCTGCCGGAtatt                                                   >  1:2624381/1‑36 (MQ=255)
                                              cgcgATGGTAAACGGTGGGGAGACCTGCCGGAtatt                                                   >  1:916877/1‑36 (MQ=255)
                                                          cGGTGGGGAGACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAATTa      <  1:1425929/69‑1 (MQ=255)
                                                          cGGTGGGGAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAAtt       <  1:1880529/68‑1 (MQ=255)
                                                              ggggAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAATTAGTgg  <  1:3216072/69‑1 (MQ=255)
                                                              ggggAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAATTAGTgg  >  1:741730/1‑69 (MQ=255)
                                                                   |                                                               
TAAAGTGGCGATGGTGGTGGGTAGCACCATTAACACCCTCAAAGATCGCGATGGTAAACGGTGGGGAAACCTGCCGGATATTTATTGCGCTTATCTTATCCGATGTCTGGTTTTTCGCGGGGAATTAGTGG  >  W3110S.gb/4390200‑4390330

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: