Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 187,851 G→A 6.0% A136T (GCG→ACG)  accA → acetylCoA carboxylase, carboxytransferase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE187,8510GA 6.0% 205.1 / 2.7 66A136T (GCG→ACG) accAacetylCoA carboxylase, carboxytransferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TGATGATCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTAAGAT  >  minE/187783‑187909
                                                                    |                                                          
tgatgaTCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcg                                                          >  1:4548250/1‑71 (MQ=255)
 gatgatCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         >  1:2659915/1‑71 (MQ=255)
   tgatCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcg                                                          >  1:4411812/1‑68 (MQ=255)
     atCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      <  1:3061953/70‑1 (MQ=255)
     atCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      <  1:3738435/70‑1 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:4250914/1‑69 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:449033/1‑69 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:2545481/1‑69 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:3245823/1‑69 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:5136706/1‑69 (MQ=255)
      tCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:3874589/1‑69 (MQ=255)
       cATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:3375372/1‑68 (MQ=255)
       cATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAg                                                   >  1:3574289/1‑71 (MQ=255)
               atcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:1893757/57‑1 (MQ=255)
               atcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:792904/57‑1 (MQ=255)
               atcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:1196606/57‑1 (MQ=255)
               atcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:4548743/57‑1 (MQ=255)
               atcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:1209312/57‑1 (MQ=255)
                tcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTa                                               <  1:4091674/66‑1 (MQ=255)
                tcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCa                                          <  1:1315546/71‑1 (MQ=255)
                tcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCa                                          <  1:3319518/71‑1 (MQ=255)
                tcaAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCa                                          <  1:2180886/71‑1 (MQ=255)
                  aaaaaGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAg                                                      >  1:140825/1‑57 (MQ=255)
                  aaaaaGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGc                                           >  1:4037559/1‑68 (MQ=255)
                  aaaaaGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCaa                                         <  1:3000009/70‑1 (MQ=255)
                     aaGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcgc                                                         <  1:518170/51‑1 (MQ=255)
                      aGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcg                                                          >  1:3612459/1‑49 (MQ=255)
                      aGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAgcg                                                          >  1:1150639/1‑49 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGc                                      <  1:4127387/66‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGc                                      <  1:901508/66‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGc                                      <  1:1250590/66‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGc                                      <  1:3763868/66‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGc                                 <  1:226447/71‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGc                                 <  1:2994626/71‑1 (MQ=255)
                         tCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGc                                 <  1:3646867/71‑1 (MQ=255)
                          cGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCg                                            >  1:4463177/1‑59 (MQ=255)
                          cGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCg                                            >  1:4280256/1‑59 (MQ=255)
                           gTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGt                               <  1:4542102/71‑1 (MQ=255)
                           gTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGt                               <  1:3023350/71‑1 (MQ=255)
                           gTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGt                               <  1:62089/71‑1 (MQ=255)
                              aaaCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGa                                                     <  1:385884/46‑1 (MQ=255)
                              aaaCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGa                                                     <  1:1865887/46‑1 (MQ=255)
                              aaaCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAACGCCAg                                                      <  1:89519/45‑1 (MQ=255)
                              aaaCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAACGCCAg                                                      <  1:4685899/45‑1 (MQ=255)
                                aCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAACGCCAGaa                                                    >  1:3868678/1‑45 (MQ=255)
                                aCCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAACGCCAGaa                                                    >  1:3377324/1‑45 (MQ=255)
                                   aaaGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGt                               <  1:1292282/63‑1 (MQ=255)
                                   aaaGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCa                       <  1:3968062/71‑1 (MQ=255)
                                     aGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCAc                                    >  1:4143565/1‑56 (MQ=255)
                                     aGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCaaa                     <  1:838484/71‑1 (MQ=255)
                                      gAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCg                                <  1:2779077/59‑1 (MQ=255)
                                       aaaaaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATg                   >  1:1650341/1‑71 (MQ=255)
                                         aaaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGc                 >  1:949362/1‑71 (MQ=255)
                                         aaaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGc                 >  1:421102/1‑71 (MQ=255)
                                         aaaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGc                 >  1:1977297/1‑71 (MQ=255)
                                          aaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCt                >  1:4917518/1‑71 (MQ=255)
                                          aaaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCt                >  1:4868056/1‑71 (MQ=255)
                                           aaTTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTg               <  1:4644875/71‑1 (MQ=255)
                                                    tAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCa                                     <  1:473018/40‑1 (MQ=255)
                                                    tAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCa                                     <  1:1153744/40‑1 (MQ=255)
                                                    tAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCttt       >  1:333147/1‑70 (MQ=255)
                                                     aaCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCa                                     >  1:1858265/1‑39 (MQ=255)
                                                     aaCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCa                                     >  1:2754090/1‑39 (MQ=255)
                                                     aaCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCa                                     >  1:2241199/1‑39 (MQ=255)
                                                      aCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCttt       >  1:2011957/1‑68 (MQ=255)
                                                      aCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCttt       >  1:670208/1‑68 (MQ=255)
                                                      aCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCttt       >  1:641780/1‑68 (MQ=255)
                                                      aCTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCtt        >  1:2383940/1‑67 (MQ=255)
                                                      aCTTTCGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCttt       >  1:5098842/1‑68 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTaa     <  1:1110170/68‑1 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTaa     <  1:2906656/68‑1 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTaa     <  1:2930361/68‑1 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTaa     <  1:3457722/68‑1 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTaa     <  1:3501623/68‑1 (MQ=255)
                                                        tttGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTAAGAt  >  1:2896742/1‑71 (MQ=255)
                                                               tGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAAt                    <  1:2261116/46‑1 (MQ=255)
                                                               tGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAAt                    <  1:2298500/46‑1 (MQ=255)
                                                                    |                                                          
TGATGATCATTGGTCATCAAAAAGGTCGTGAAACCAAAGAAAAAATTCGCCGTAACTTTGGTATGCCAGCGCCAGAAGGTTACCGCAAAGCACTGCGTCTGATGCAAATGGCTGAACGCTTTAAGAT  >  minE/187783‑187909

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: