Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,002,862 (T)7→8 intergenic (‑87/‑145) pitA ← / → yhiN phosphate transporter, low‑affinity/predicted oxidoreductase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,002,8551.T98.4% 211.3 / ‑3.2 64intergenic (‑80/‑152)pitA/yhiNphosphate transporter, low‑affinity/predicted oxidoreductase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (33/30);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGG  >  W3110S.gb/4002806‑4002917
                                                  |                                                              
gcggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGc                                                  <  1:36732/65‑1 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:957883/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:2093991/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:2303370/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:2601735/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:2629867/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:1451344/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:2853453/1‑65 (MQ=255)
 cggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGcc                                                 >  1:327112/1‑65 (MQ=255)
   gcgcATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGccc                                                <  1:1398835/64‑1 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGa                                      >  1:1604940/1‑63 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2878960/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2801272/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:3069287/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2710193/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2603383/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:371463/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2470766/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:78773/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2286700/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:216187/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2272811/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2004833/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:1097067/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:149605/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:1422251/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:2004831/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGAt                                     >  1:167616/1‑64 (MQ=255)
              gTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCACTGTATGGAt                                     >  1:2217872/1‑64 (MQ=38)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:190826/1‑64 (MQ=255)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:868540/1‑64 (MQ=255)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:115238/1‑64 (MQ=255)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:772732/1‑64 (MQ=255)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:612407/1‑64 (MQ=255)
                     tCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTc                              >  1:242400/1‑64 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1267490/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1284359/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:3012238/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:3010292/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1461377/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:316375/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1263290/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1262529/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1227072/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:638616/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:699465/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:728433/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1293656/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2033339/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2106233/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:215749/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1802592/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2232660/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1778210/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2398090/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1637006/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1566872/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2625504/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2664953/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:283997/65‑1 (MQ=255)
                          gcgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2866091/65‑1 (MQ=255)
                           cgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:1735876/64‑1 (MQ=255)
                           cgcggGCAAAGTGAAATATAGACTTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCaaaaaa                        <  1:2138511/64‑1 (MQ=255)
                                               gACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAgg  <  1:1020121/65‑1 (MQ=255)
                                                  |                                                              
GCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGG  >  W3110S.gb/4002806‑4002917

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: