Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,824,257 T→C 13.3% F29S (TTT→TCT)  nadE → NAD synthetase, NH3/glutamine‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,824,2570TC13.3% 70.4 / 5.1 30F29S (TTT→TCT) nadENAD synthetase, NH3/glutamine‑dependent
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base C (2/2);  total (15/15)
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.90e-01

TAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGA  >  W3110S.gb/1824191‑1824326
                                                                  |                                                                     
taataaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCt                                                                    >  1:1024939/1‑70 (MQ=255)
       ggCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCt                                                            <  1:348268/71‑1 (MQ=255)
       ggCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCt                                                            <  1:1532657/71‑1 (MQ=255)
        gcgcTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCt                                                            <  1:1222441/70‑1 (MQ=255)
               gcgcAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCaa                                                    <  1:285647/71‑1 (MQ=255)
                             aTTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCt                                                            <  1:2176043/49‑1 (MQ=255)
                                aaTGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCt                                                            <  1:2526648/46‑1 (MQ=255)
                                aaTGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTa                                    <  1:700487/70‑1 (MQ=255)
                                    cTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTAc                                                          >  1:2278963/1‑44 (MQ=255)
                                    cTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTAc                                                          >  1:2282194/1‑44 (MQ=255)
                                       aaGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCAc                              <  1:1638445/69‑1 (MQ=255)
                                               aTTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCgg                    >  1:2312667/1‑71 (MQ=255)
                                                ttCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGtt                                         <  1:2765203/49‑1 (MQ=255)
                                                ttCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCgg                    >  1:1520647/1‑70 (MQ=255)
                                                ttCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCggg                   <  1:8698/71‑1 (MQ=255)
                                                      gtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgc             <  1:1925938/71‑1 (MQ=255)
                                                         gtgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTaa                                   <  1:1785783/46‑1 (MQ=255)
                                                         gtgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCggg                   <  1:2268302/62‑1 (MQ=255)
                                                         gtgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCggg                   <  1:2271363/62‑1 (MQ=255)
                                                          tgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTc                      <  1:2313951/58‑1 (MQ=255)
                                                            tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCAt                                      >  1:1450425/1‑40 (MQ=255)
                                                            tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCAt                                      >  1:1218468/1‑40 (MQ=255)
                                                            tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCAGTgg                           >  1:2011973/1‑51 (MQ=255)
                                                            tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:2718509/1‑71 (MQ=255)
                                                            tCGATTCTCTGAAAAGCTACCTGCAATCTTATCCGTTCATTAAATCAc                              >  1:1839442/1‑48 (MQ=255)
                                                              gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGcggcgg        >  1:2620304/1‑68 (MQ=255)
                                                              gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:149334/1‑69 (MQ=255)
                                                              gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:1340230/1‑69 (MQ=255)
                                                              gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:660054/1‑69 (MQ=255)
                                                              gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:1065531/1‑69 (MQ=255)
                                                                 tttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCATGa  >  1:85649/1‑71 (MQ=255)
                                                                  ttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTg                            >  1:1253014/1‑44 (MQ=255)
                                                                  |                                                                     
TAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGA  >  W3110S.gb/1824191‑1824326

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: