Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 907,713 G→A 33.3% T331I (ACC→ATC)  ybjT ← conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb907,7130GA33.3% 34.9 / 24.7 30T331I (ACC→ATC) ybjTconserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base A (5/5);  total (16/14)
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.67e-01

ACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTA  >  W3110S.gb/907646‑907780
                                                                   |                                                                   
aCGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTa                                                                   <  1:4267886/70‑1 (MQ=255)
     ctTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAAcccc                                                             >  1:4762872/1‑71 (MQ=255)
     ctTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAAcccc                                                             >  1:3587864/1‑71 (MQ=255)
        tACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCg                                                            <  1:2676225/69‑1 (MQ=255)
        tACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGcc                                                          >  1:2765085/1‑71 (MQ=255)
        tACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATaaa                                                                 >  1:2188585/1‑64 (MQ=255)
                     gTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAAt                                             <  1:4500957/71‑1 (MQ=255)
                        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAAcccc                                                             <  1:2500913/52‑1 (MQ=255)
                        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAAc                                          >  1:4702469/1‑71 (MQ=255)
                        cACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATaa                                           <  1:1384081/70‑1 (MQ=255)
                          cTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGc                                                           >  1:984699/1‑52 (MQ=255)
                          cTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCg                                        <  1:681149/71‑1 (MQ=255)
                            aCCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGt                                       >  1:4195070/1‑70 (MQ=255)
                              cTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACt                                    <  1:1215352/71‑1 (MQ=255)
                              cTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACt                                    <  1:4248444/71‑1 (MQ=255)
                                gCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTcc                                  >  1:1767191/1‑71 (MQ=255)
                                 ccATAAAGCAGCAAGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTcc                                  <  1:850419/70‑1 (MQ=255)
                                 ccATAAAGCAGCAAGGCTAGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTcc                                  <  1:3585995/70‑1 (MQ=255)
                                  cATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCgg                                >  1:1886989/1‑71 (MQ=255)
                                  cATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCgg                                <  1:115444/71‑1 (MQ=255)
                                   aTAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGt                               >  1:3087919/1‑71 (MQ=255)
                                      aaGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGc                            >  1:1708845/1‑71 (MQ=255)
                                            caAGGCTGGCGGACGTTTTAACGGTAAACCCCGGCTGTTTGGCAAAATAAc                                          >  1:2572915/1‑51 (MQ=255)
                                             aaGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTAc                                     >  1:745855/1‑55 (MQ=255)
                                             aaGGCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCa                          >  1:1438320/1‑66 (MQ=255)
                                               ggCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCa                   <  1:3088618/71‑1 (MQ=255)
                                                gCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                  <  1:1920712/71‑1 (MQ=255)
                                                gCTGGCGGACGTTTTAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCa                   >  1:2641175/1‑70 (MQ=255)
                                                             ttAACGATAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCa                          >  1:1595526/1‑50 (MQ=255)
                                                                aCGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTa  <  1:3900620/71‑1 (MQ=255)
                                                                   |                                                                   
ACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTA  >  W3110S.gb/907646‑907780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: