Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 349,960 C→T 56.4% P155L (CCC→CTC)  hyi → hydroxypyruvate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE349,9600CT56.4% 17.6 / 65.9 55P155L (CCC→CTC) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base T (13/18);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGA  >  minE/349893‑350029
                                                                   |                                                                     
cTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCa                                                                     <  1:1748389/70‑1 (MQ=255)
     gaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGaa                                                               >  1:2860933/1‑71 (MQ=255)
       aagaagATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACt                                                             <  1:1309045/71‑1 (MQ=255)
       aagaagATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACt                                                             >  1:24312/1‑71 (MQ=255)
       aagaagATATTTTATTACTGATTGAACCTATTAACCATTTTGAAATTCCTGGTTTCCATCCCACCGGAACt                                                             <  1:2973174/71‑1 (MQ=255)
           agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                         <  1:568899/71‑1 (MQ=255)
           agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                         <  1:3818376/71‑1 (MQ=255)
           agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCgg                                                          >  1:2274141/1‑70 (MQ=255)
           agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGc                                                         <  1:102913/71‑1 (MQ=255)
                ttttATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACt                                                             <  1:2322704/62‑1 (MQ=255)
                ttttATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACt                                                             <  1:236470/62‑1 (MQ=255)
                  ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGgcgc                                                   <  1:3474991/70‑1 (MQ=255)
                   tattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCt                                                  >  1:2047750/1‑70 (MQ=255)
                     ttaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGa                                                <  1:1840933/70‑1 (MQ=255)
                     ttaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGa                                                <  1:1174318/70‑1 (MQ=255)
                         tGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCa                                                        >  1:2404434/1‑58 (MQ=255)
                         tGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCa                                                        >  1:657771/1‑58 (MQ=255)
                           aTTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattgat                                         >  1:2519822/1‑71 (MQ=255)
                           aTTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattg                                           <  1:2928148/69‑1 (MQ=255)
                           aTTGAACCTATTAACCATTTTGATATTCCTGGGTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattg                                           <  1:3450762/69‑1 (MQ=255)
                             tGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGa                                                >  1:1921451/1‑62 (MQ=255)
                                 ccTATTAACCATTTTGATATTCCTGGTTTCCATCTCAc                                                                    <  1:3730607/38‑1 (MQ=255)
                                 ccTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATtgatga                                   <  1:1230040/71‑1 (MQ=255)
                                    aTTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCt                                                  <  1:1960218/53‑1 (MQ=255)
                                     ttAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                                    <  1:1647709/66‑1 (MQ=255)
                                       aaCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttgg                             >  1:2126811/1‑71 (MQ=255)
                                       aaCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttg                              <  1:1356159/70‑1 (MQ=255)
                                        aCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttgg                             <  1:225850/70‑1 (MQ=255)
                                          cATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGaaa                                              >  1:2938307/1‑51 (MQ=255)
                                          cATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggttg                          <  1:1703518/71‑1 (MQ=255)
                                          cATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggttg                          <  1:2509976/71‑1 (MQ=255)
                                             tttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggt                            >  1:1881712/1‑66 (MQ=255)
                                              ttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGgcg                                                    >  1:869941/1‑41 (MQ=255)
                                              ttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctg                       >  1:34166/1‑70 (MQ=255)
                                               tGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctg                       >  1:2919281/1‑69 (MQ=255)
                                                gATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCa                     <  1:2875392/70‑1 (MQ=255)
                                                gATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCa                     <  1:2687120/70‑1 (MQ=255)
                                                 atatTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt                   <  1:2867589/71‑1 (MQ=255)
                                                  tatTCCTGGTTTCCATCTCACCGGAACTCGGCAgg                                                      >  1:3636027/1‑35 (MQ=255)
                                                    ttCCTGGTTACCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAttt                 <  1:1438856/70‑1 (MQ=255)
                                                    ttCATGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAttt                 >  1:1158747/1‑70 (MQ=255)
                                                      ccTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                                    >  1:1393544/1‑49 (MQ=255)
                                                       cTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaa              <  1:3087138/70‑1 (MQ=255)
                                                         ggTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAt           >  1:1331072/1‑71 (MQ=255)
                                                          gTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttgg                             >  1:3616178/1‑52 (MQ=255)
                                                             tCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:1603563/1‑71 (MQ=255)
                                                             tCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:1514878/1‑71 (MQ=255)
                                                             tCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:581913/1‑71 (MQ=255)
                                                             tCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:281099/1‑71 (MQ=255)
                                                             tCCATACCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:824338/1‑71 (MQ=255)
                                                             tCCATACCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:828791/1‑71 (MQ=255)
                                                              ccATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       >  1:3339411/1‑70 (MQ=255)
                                                              ccATCCCACAGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg       <  1:3559172/70‑1 (MQ=255)
                                                               cATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGt      <  1:639935/70‑1 (MQ=255)
                                                                aTCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt                   <  1:3751938/56‑1 (MQ=255)
                                                                   ccACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGa  >  1:2126883/1‑70 (MQ=255)
                                                                   |                                                                     
CTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGA  >  minE/349893‑350029

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: