Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,993,572 G→A 35.9% V289M (GTG→ATG)  nupG → nucleoside transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,993,5720GA35.9% 121.8 / 124.2 128V289M (GTG→ATG) nupGnucleoside transporter
Reads supporting (aligned to +/- strand):  ref base G (60/22);  new base A (18/28);  total (78/50)
Fisher's exact test for biased strand distribution p-value = 2.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 128 total aligned reads displayed.

CGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCG  >  minE/1993536‑1993626
                                    |                                                      
cGCTACGGTATTAAGAACGTGATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1888849/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:398871/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1938482/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:2066652/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:844075/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:2128108/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:2136223/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:896340/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1784589/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:248366/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:169114/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1689649/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:30656/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:15726/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1540737/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1468312/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:93537/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1168997/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1203729/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1157726/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1132174/67‑1 (MQ=255)
cGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:127573/67‑1 (MQ=255)
 gCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                          <  1:1661972/66‑1 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATTGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:634848/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                  >  1:1559936/1‑63 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCgcc               >  1:263737/1‑64 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:733942/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:675212/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:785228/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:642774/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:1050756/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:1955178/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:736669/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:1689796/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:390432/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:44107/1‑66 (MQ=255)
            aaGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc               >  1:1270903/1‑66 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTg                                          >  1:451766/1‑37 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAcctg            >  1:843816/1‑65 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc             >  1:1413386/1‑66 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAcc              >  1:2063385/1‑65 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:229074/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:619797/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:367171/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:2249500/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1884471/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:651840/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:44503/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1938770/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1819412/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1625441/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1656657/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1901246/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1734876/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1756711/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1770632/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1624190/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:929112/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1522265/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:848648/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:967161/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:117827/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1132749/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1089614/1‑67 (MQ=255)
              gAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:1048980/1‑67 (MQ=255)
              gAACGTAATGATGATCACTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCg            >  1:2278383/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:277704/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:1597040/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:678834/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:69654/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:1816941/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:2199073/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:2168632/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:874155/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:2087719/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:1890518/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:797815/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:1898023/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:799378/1‑67 (MQ=255)
                       gatgatCAGTATTGTGGAGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc   >  1:1106336/1‑67 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:962853/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:740614/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:982474/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1020111/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:786912/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:633258/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:2141171/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:2006539/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1968427/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1954412/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1920084/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1867836/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1857561/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1825032/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1809228/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1586866/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1381906/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1350240/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:126078/67‑1 (MQ=255)
                        atgatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg  <  1:1122473/67‑1 (MQ=255)
                                    |                                                      
CGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCG  >  minE/1993536‑1993626

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: