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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE349,9600CT62.5% 54.0 / 69.1 72P155L (CCC→CTC) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base C (12/15);  new base T (17/28);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAG  >  minE/349891‑350024
                                                                     |                                                                
tGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATcc                                                                  <  1:3084683/70‑1 (MQ=255)
    gatgaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCg                                                             >  1:2128141/1‑71 (MQ=255)
     atgaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGAAATTCCTGGTTTCCATCCCACCgg                                                            <  1:4434988/71‑1 (MQ=255)
           gaagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCg                                                      >  1:5323918/1‑71 (MQ=255)
            aagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCgg                                                     <  1:4510585/71‑1 (MQ=255)
            aagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCgg                                                     <  1:1837973/71‑1 (MQ=255)
             agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                    <  1:410145/71‑1 (MQ=255)
             agaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                    <  1:2776779/71‑1 (MQ=255)
              gaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCg                                                      >  1:4462012/1‑68 (MQ=255)
              gaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCg                                                      >  1:4171305/1‑68 (MQ=255)
              gaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCg                                                      >  1:5428596/1‑68 (MQ=255)
              gaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCa                                                   <  1:3250166/71‑1 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCt                                             >  1:2108580/1‑71 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                               >  1:4541229/1‑53 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                               >  1:3797272/1‑53 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                               >  1:106318/1‑53 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                               >  1:647987/1‑53 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                               >  1:1276846/1‑53 (MQ=255)
                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGgcgc                                              <  1:2635219/70‑1 (MQ=255)
                     tattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCt                                             <  1:3965963/70‑1 (MQ=255)
                       ttaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGa                                           <  1:153782/70‑1 (MQ=255)
                              ttGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCgg                                                     <  1:238731/53‑1 (MQ=255)
                                gAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattgattg                                  <  1:5645733/70‑1 (MQ=255)
                                gAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattgattg                                  <  1:1561460/70‑1 (MQ=255)
                                 aaCCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                <  1:1674394/71‑1 (MQ=255)
                                  aCCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGa                                                           <  1:2145037/43‑1 (MQ=255)
                                  aCCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGa                                                           <  1:3210566/43‑1 (MQ=255)
                                  aCCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                               <  1:1315247/71‑1 (MQ=255)
                                    cTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGaaa                                         <  1:3873809/59‑1 (MQ=255)
                                     tATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAat                                        <  1:2055325/59‑1 (MQ=255)
                                     tATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAat                                        <  1:4433608/59‑1 (MQ=255)
                                      aTTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgt                           >  1:1356333/1‑71 (MQ=255)
                                       ttAACCATTTTGATATTCCTGGTTTCCATCTCACCgg                                                            <  1:5573606/37‑1 (MQ=255)
                                       ttAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCgg                                                     <  1:444620/44‑1 (MQ=255)
                                       ttAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgtt                          >  1:3187932/1‑71 (MQ=255)
                                         aaCCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                    <  1:424486/43‑1 (MQ=255)
                                         aaCCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                    <  1:1291537/43‑1 (MQ=255)
                                         aaCCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGc                                                    <  1:1310847/43‑1 (MQ=255)
                                         aaCCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCa                                                   >  1:4979480/1‑44 (MQ=255)
                                         aaCCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATtgatga                              >  1:356831/1‑65 (MQ=255)
                                           ccATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                               <  1:2255553/62‑1 (MQ=255)
                                            cATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCt                                             <  1:2655740/47‑1 (MQ=255)
                                             aTTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGaaa                                         >  1:2365573/1‑50 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGgcg                                               >  1:1461813/1‑43 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGa                                           <  1:4691159/47‑1 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGa                                           <  1:417159/47‑1 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                <  1:1070845/58‑1 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggt                       <  1:5512542/67‑1 (MQ=255)
                                              ttttGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATtgatga                              <  1:3733347/60‑1 (MQ=255)
                                               tttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCa                                                   <  1:819388/38‑1 (MQ=255)
                                               tttGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctg                  >  1:2654870/1‑71 (MQ=255)
                                               tttGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAat                                        >  1:1777536/1‑49 (MQ=255)
                                                  gATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGatgttggtt                         <  1:1022585/61‑4 (MQ=255)
                                                  gATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATAGatgttggtt                         <  1:4770937/61‑4 (MQ=255)
                                                   atatTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt              <  1:5195275/71‑1 (MQ=255)
                                                   atatTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt              <  1:97888/71‑1 (MQ=255)
                                                   atatTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt              <  1:4522861/71‑1 (MQ=255)
                                                    tatTCCTGGTTTCCATCTCACCGGAACTCGGCAGgcgc                                              >  1:3183470/1‑38 (MQ=255)
                                                    tatTCCTGGTTTCCATCTCACCGGAACTCGGCAGgcgc                                              >  1:1125974/1‑38 (MQ=255)
                                                      ttCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAattg                                      <  1:1283683/44‑1 (MQ=255)
                                                        ccTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGa          <  1:2398437/70‑1 (MQ=255)
                                                        ccTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                >  1:2635514/1‑48 (MQ=255)
                                                        ccTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                >  1:1827473/1‑48 (MQ=255)
                                                        ccTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGa          <  1:3029184/70‑1 (MQ=255)
                                                           ggTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaaaa       >  1:2277117/1‑70 (MQ=255)
                                                           ggTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAt      >  1:500080/1‑71 (MQ=255)
                                                            gTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAt      <  1:5138495/70‑1 (MQ=255)
                                                               tCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCa   <  1:2532224/70‑1 (MQ=255)
                                                                ccATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAg  >  1:1544460/1‑70 (MQ=255)
                                                                ccATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttg                         >  1:4887503/1‑47 (MQ=255)
                                                                ccATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttg                         >  1:772974/1‑47 (MQ=255)
                                                                  aTCCCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                <  1:5487573/38‑1 (MQ=255)
                                                                   tCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgt                           >  1:1965698/1‑42 (MQ=255)
                                                                     ccACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                               >  1:5064456/1‑36 (MQ=255)
                                                                     ccACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaa         <  1:2259166/58‑1 (MQ=255)
                                                                     |                                                                
TGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAG  >  minE/349891‑350024

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: