Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,102,502 T→C 7.1% F42F (TTT→TTC aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,102,5020TC 7.1% 153.5 / 3.1 56F42F (TTT→TTCaroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (17/35);  new base C (2/2);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATGGCGG  >  minE/1102439‑1102566
                                                               |                                                                
gCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATAtt                                                            <  1:1724725/70‑1 (MQ=255)
gCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATAtt                                                            <  1:1649164/70‑1 (MQ=255)
    aTGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGa                                                       <  1:1332209/71‑1 (MQ=255)
     tGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGaa                                                      <  1:2990502/71‑1 (MQ=255)
        ctaaaGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATgg                                                   >  1:1688741/3‑71 (MQ=255)
         gAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTCGATATTCTGGAATGGc                                                  <  1:2873579/71‑1 (MQ=255)
         gAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTCGATATTCTGGAATGGc                                                  <  1:3183904/71‑1 (MQ=255)
           aaGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATg                                                    >  1:1986589/1‑67 (MQ=255)
            aGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTg                                                         >  1:116873/1‑61 (MQ=255)
            aGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATg                                                    <  1:3248834/66‑1 (MQ=255)
              atatCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCtcg                                                        >  1:3949393/1‑58 (MQ=255)
               tatCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGTCTTTGATATTCTGGAATg                                                    >  1:541641/1‑63 (MQ=255)
               tatCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGtgtg                                             <  1:3548342/70‑1 (MQ=255)
                   gCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCg                                                 <  1:2322615/62‑1 (MQ=255)
                   gCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCa                                        <  1:2632548/71‑1 (MQ=255)
                   gCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCa                                        <  1:2580245/71‑1 (MQ=255)
                    ccAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTg                                                         <  1:474588/53‑1 (MQ=255)
                     cAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATAtt                                                            <  1:45094/49‑1 (MQ=255)
                     cAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGtgtg                                             <  1:316982/64‑1 (MQ=255)
                      aGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATAtt                                                            <  1:2676444/48‑1 (MQ=255)
                      aGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTa                                     <  1:3230437/71‑1 (MQ=255)
                      aGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTa                                     <  1:494037/71‑1 (MQ=255)
                      aGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTa                                     <  1:1654465/71‑1 (MQ=255)
                       gCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGat                                                               >  1:3947823/1‑44 (MQ=255)
                        cGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATg                                   >  1:2806927/1‑71 (MQ=255)
                        cGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATg                                   <  1:35060/71‑1 (MQ=255)
                         gTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGtgtg                                             <  1:3323431/60‑1 (MQ=255)
                          tGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGcc                                 >  1:3259637/1‑71 (MQ=255)
                          tGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGcc                                 <  1:994467/71‑1 (MQ=255)
                            aaaTCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCg                                <  1:3167556/70‑1 (MQ=255)
                             aaTCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGAc                              >  1:3673310/1‑71 (MQ=255)
                              aTCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTgg                                                        <  1:1851695/44‑1 (MQ=255)
                                ccGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGaa                                                      >  1:3814634/1‑44 (MQ=255)
                                   aaGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCtctc                         >  1:2843100/1‑70 (MQ=255)
                                   aaGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTcc                        <  1:1094597/71‑1 (MQ=255)
                                    aGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCa                       <  1:961549/71‑1 (MQ=255)
                                      ctctCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAt                     <  1:2237551/71‑1 (MQ=255)
                                         tcGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGa                               >  1:1544073/1‑58 (MQ=255)
                                         tcGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAtgtg                  >  1:426251/1‑71 (MQ=255)
                                         tcGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAtg                    <  1:80234/69‑1 (MQ=255)
                                           gCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAtgtg                  <  1:932337/69‑1 (MQ=255)
                                           gCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGa                <  1:1251990/71‑1 (MQ=255)
                                           gCCTATCGTGAAGCGGACTTCGATATTCTGGAATGGCGTGTGGa                                           >  1:2585012/1‑44 (MQ=255)
                                           gCCTATCGTGAAGCGGACTTCGATATTCTGGAATGGCGTGTGGAc                                          >  1:1048390/1‑45 (MQ=255)
                                             cTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAt                     >  1:3885462/1‑64 (MQ=255)
                                              tATCGTGAAGCGGACTTTGATATTCTGGAATGGCgtgt                                              <  1:767894/38‑1 (MQ=255)
                                                    gAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCAt       >  1:2549757/1‑71 (MQ=255)
                                                     aaGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCAt       <  1:619583/70‑1 (MQ=255)
                                                     aaGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATg      <  1:2470843/71‑1 (MQ=255)
                                                       gCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAAt                     <  1:244540/54‑1 (MQ=255)
                                                        cGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGAcc                             <  1:763874/45‑1 (MQ=255)
                                                         ggACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGa                <  1:1583703/57‑1 (MQ=255)
                                                         ggACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGa                <  1:2378391/57‑1 (MQ=255)
                                                         ggACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATggcgg  <  1:1651369/71‑1 (MQ=255)
                                                           aCTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCAt       >  1:1354081/1‑64 (MQ=255)
                                                            cTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCt            >  1:1117190/1‑58 (MQ=255)
                                                             tttGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGt          <  1:1986390/59‑1 (MQ=255)
                                                             tttGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCAt       <  1:2274343/62‑1 (MQ=255)
                                                             tttGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCAt       <  1:1871847/62‑1 (MQ=255)
                                                             tttGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATggcgg  <  1:524841/67‑1 (MQ=255)
                                                             tttGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATggcgg  <  1:3575078/67‑1 (MQ=255)
                                                               tGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCa        >  1:1356865/1‑59 (MQ=255)
                                                               tGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTATGTCa        >  1:1490980/1‑59 (MQ=255)
                                                               |                                                                
GCTGATGGCGAAAGATATCGCCAGCGTGAAATCCGAAGCTCTCGCCTATCGTGAAGCGGACTTTGATATTCTGGAATGGCGTGTGGACCACTATGCCGACCTCTCCAATGTGGAGTCTGTCATGGCGG  >  minE/1102439‑1102566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: