Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,824,257 T→C 32.6% 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,2570TC32.6% 55.4 / 42.8 49F29S (TTT→TCT) nadENAD synthetase, NH3/glutamine‑dependent
Reads supporting (aligned to +/- strand):  ref base T (17/16);  new base C (5/11);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CAAATAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGA  >  W3110S.gb/1824187‑1824326
                                                                      |                                                                     
caaATAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGAttt                                                                       <  1:3069693/71‑1 (MQ=255)
       taaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGaaa                                                                <  1:1899240/71‑1 (MQ=255)
       taaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGaaa                                                                <  1:2385446/71‑1 (MQ=255)
         aaGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAg                                                              <  1:2941629/71‑1 (MQ=255)
          aGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGc                                                             <  1:2667409/71‑1 (MQ=255)
            gcgcTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTa                                                           <  1:2196547/71‑1 (MQ=255)
                       aaaaCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACtt                                                <  1:258848/71‑1 (MQ=255)
                        aaaCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACtt                                                >  1:667235/1‑70 (MQ=255)
                         aaCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTAt                                              >  1:2188528/1‑71 (MQ=255)
                          aCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTAt                                              >  1:206878/1‑70 (MQ=255)
                          aCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATc                                             >  1:1056423/1‑71 (MQ=255)
                          aCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATc                                             <  1:4610065/71‑1 (MQ=255)
                            cGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCg                                           >  1:577994/1‑71 (MQ=255)
                                   tAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCAtt                                     >  1:727318/1‑70 (MQ=255)
                                   tAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTAc                                                          >  1:3329101/1‑49 (MQ=255)
                                      tGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGaa                                                                 <  1:2341195/39‑1 (MQ=255)
                                      tGCTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTa                                    <  1:2122606/68‑1 (MQ=255)
                                        cTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCaaa                                                   >  1:2984871/1‑51 (MQ=255)
                                        cTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCa                               >  1:1361350/1‑71 (MQ=255)
                                        cTGAAGAGGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAAt                                 <  1:3320409/69‑1 (MQ=255)
                                          gAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAg                                                              <  1:561753/38‑1 (MQ=255)
                                           aaGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTAc                                                          >  1:694873/1‑41 (MQ=255)
                                              agGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCa                                                     <  1:1122691/43‑1 (MQ=255)
                                              agGAAATTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCa                                                     <  1:896442/43‑1 (MQ=255)
                                                gAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATc                                <  1:274291/62‑1 (MQ=255)
                                                gAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGc                        >  1:2611925/1‑70 (MQ=255)
                                                 aaaTTCGTCGTAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTc                      >  1:2274889/1‑71 (MQ=255)
                                                    ttCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGtt                                         >  1:611087/1‑49 (MQ=255)
                                                    ttCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGt                                          <  1:256943/48‑1 (MQ=255)
                                                       gtcgtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATc                <  1:3328475/71‑1 (MQ=255)
                                                          gtAGTGTCGATTTTCTGAAAAGCTACCTGCGAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgc             >  1:2424171/1‑71 (MQ=255)
                                                          gtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgc             >  1:117021/1‑71 (MQ=255)
                                                          gtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgc             >  1:4278517/1‑71 (MQ=255)
                                                          gtAGTGTCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgc             <  1:2752517/71‑1 (MQ=255)
                                                             gtgtCGATTTTCTGAAAAGCTACCTGCAAACTTATc                                             <  1:2743159/36‑1 (MQ=255)
                                                             gtgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTa                                    <  1:334592/45‑1 (MQ=255)
                                                              tttCGATTTTCTGAAAAGCTACCTCCAAACTTATCCGTTCa                                       <  1:869447/39‑1 (MQ=38)
                                                              tgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCAc                              <  1:2530103/50‑1 (MQ=255)
                                                              tgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGAt                 <  1:2979472/63‑1 (MQ=255)
                                                              tgtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAg              <  1:4083911/66‑1 (MQ=255)
                                                               gtCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGtt                                         <  1:695043/38‑1 (MQ=255)
                                                                tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCa                                       >  1:2301514/1‑39 (MQ=255)
                                                                tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAAt                                 >  1:2449044/1‑45 (MQ=255)
                                                                tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTgg                           >  1:3535353/1‑51 (MQ=255)
                                                                tCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGAt                 <  1:3712394/61‑1 (MQ=255)
                                                                tCGATTCTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:4137774/1‑71 (MQ=255)
                                                                  gATTTTCTGAAAATCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGa                  >  1:2174692/1‑58 (MQ=255)
                                                                  gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCggg                   <  1:3693432/57‑1 (MQ=255)
                                                                  gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGt       >  1:1341912/1‑69 (MQ=255)
                                                                  gATTCTCTGAAAAGCTACCTGCAAACATATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTc      <  1:4354453/70‑1 (MQ=255)
                                                                    ttttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTaa                                   <  1:3785552/39‑1 (MQ=255)
                                                                    ttttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTaa                                   <  1:3409563/39‑1 (MQ=255)
                                                                    ttttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgcggc          <  1:2782074/64‑1 (MQ=255)
                                                                    ttttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAgcggc          <  1:2186641/64‑1 (MQ=255)
                                                                     tttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCa                               >  1:1588235/1‑42 (MQ=255)
                                                                      ttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGa  >  1:2694644/1‑70 (MQ=255)
                                                                      |                                                                     
CAAATAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGA  >  W3110S.gb/1824187‑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: