Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,065,027 G→T 55.8% A145D (GCT→GAT)  erfK ← 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.gb2,065,0270GT55.8% 4.7 / 55.1 34A145D (GCT→GAT) erfKconserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (11/4);  new base T (11/8);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGGCAAACTCTCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCGCCCAGCCTGGCCGA  >  W3110S.gb/2064959‑2065090
                                                                    |                                                               
aTGCGGGCAAACTCTCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATct                                                               >  1:1688612/1‑71 (MQ=255)
   cGGGCAAACTCTCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGcttct                                                            >  1:1435113/1‑71 (MQ=255)
       cAAACTCTCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGtt                                                        <  1:2056364/71‑1 (MQ=255)
           ctctCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACg                                                    >  1:1260488/1‑71 (MQ=255)
                     gTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTa                                          <  1:2260497/71‑1 (MQ=255)
                           cATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCAcc                                    >  1:89440/1‑71 (MQ=255)
                            atatTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGtt                                                        >  1:1205001/1‑50 (MQ=255)
                             tatTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCAcc                                    >  1:2131135/1‑69 (MQ=255)
                             tatTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCa                                  <  1:1116244/71‑1 (MQ=255)
                             tatTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCa                                  <  1:1523685/71‑1 (MQ=255)
                               ttCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTTGAGCTTCTTTTTTACGTTCAACGGTAGTCACCCAGt                                <  1:1399530/71‑1 (MQ=255)
                               ttCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGt                                >  1:1921515/1‑71 (MQ=255)
                               ttCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGt                                >  1:575084/1‑71 (MQ=255)
                               ttCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCatt                                <  1:2609776/71‑3 (MQ=255)
                               ttCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGt                                >  1:1675479/1‑71 (MQ=255)
                                tCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGt                                                         >  1:2553646/1‑45 (MQ=255)
                                tCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGt                                >  1:2595201/1‑70 (MQ=255)
                                tCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGt                                >  1:2057050/1‑70 (MQ=255)
                                 cgcgCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGtt                                                        >  1:2724063/1‑45 (MQ=255)
                                 cgcgCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTa                                          >  1:378845/1‑59 (MQ=255)
                                            ttCGGCTTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCa                                                <  1:2777780/42‑1 (MQ=38)
                                            ttCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGttt                    >  1:1286051/1‑70 (MQ=255)
                                            ttCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTtc                   >  1:704697/1‑71 (MQ=255)
                                              cGGCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAAc                                              >  1:1811649/1‑42 (MQ=255)
                                              cGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTtctc                 >  1:998018/1‑71 (MQ=255)
                                               ggCGTTGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCg                <  1:992564/71‑1 (MQ=255)
                                                   ttGGCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCGCCCa            <  1:2241088/71‑1 (MQ=255)
                                                     ggCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCa                                      >  1:1941184/1‑43 (MQ=255)
                                                      gCGTCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCGCCCAGc          <  1:1193927/70‑1 (MQ=255)
                                                         tCCAGGTTGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTtctc                 >  1:1787689/1‑60 (MQ=255)
                                                               ttGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGttt                    <  1:1738910/51‑1 (MQ=255)
                                                               ttGGATCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCGCCCAGCCTGGCCGa  >  1:656277/1‑69 (MQ=255)
                                                                 ggAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGt                      >  1:1021241/1‑47 (MQ=255)
                                                                  gAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTAcgcg                          <  1:7028/42‑1 (MQ=255)
                                                                    |                                                               
ATGCGGGCAAACTCTCCCCTCGTTTCGCATATTCGCGCCGAGTGTTCGGCGTTGGCGTCCAGGTTGGAGCTTCTTGTTTACGTTCAACGGTAGTCACCCAGTTACGCGGGGTTTCTCGCCCAGCCTGGCCGA  >  W3110S.gb/2064959‑2065090

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: