Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,869,845 G→A 14.9% intergenic (‑90/‑229) ygcW ← / → yqcE predicted deoxygluconate dehydrogenase/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,869,8450GA14.9% 124.9 / 12.5 47intergenic (‑90/‑229)ygcW/yqcEpredicted deoxygluconate dehydrogenase/predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (15/25);  new base A (3/4);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

AAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTG  >  minE/1869775‑1869912
                                                                      |                                                                   
aaGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTg                                                                     >  1:1814645/1‑71 (MQ=255)
    cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTg                                                                     >  1:2239880/1‑67 (MQ=255)
      aaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATa                                                               <  1:1804641/71‑1 (MQ=255)
         aaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAtt                                                            <  1:2455263/71‑1 (MQ=255)
          aaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAGCTTTCTCATTATCTCTGGTAACTGGAGATAAttt                                                           <  1:3046743/71‑1 (MQ=255)
                 agGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttttt                                                         >  1:2418009/1‑66 (MQ=255)
                 agGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttttt                                                         >  1:1892464/1‑66 (MQ=255)
                    tAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAtttt                                                          <  1:63067/62‑1 (MQ=18)
                     aaCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAgg                                                >  1:2417531/1‑71 (MQ=255)
                       cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTg                                                                     <  1:2274610/48‑1 (MQ=255)
                       cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGc                                              <  1:971573/71‑1 (MQ=255)
                       cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGc                                              <  1:3019708/71‑1 (MQ=255)
                             ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATa                                                               <  1:3161944/48‑1 (MQ=18)
                             ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATa                                                               <  1:358919/48‑1 (MQ=18)
                             ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATa                                                               <  1:540346/48‑1 (MQ=18)
                              cgcgAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATa                                                               <  1:2109267/47‑1 (MQ=17)
                                   gagGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCAt                                                    <  1:2046539/53‑1 (MQ=255)
                                   gagGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCAt                                                    <  1:2564393/53‑1 (MQ=255)
                                    agGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATaa                                                  <  1:302233/54‑1 (MQ=255)
                                    agGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATaa                                                  <  1:3494779/54‑1 (MQ=255)
                                     gggAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtg                                 <  1:451818/70‑1 (MQ=255)
                                        aaGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTaa                                   >  1:2668396/1‑65 (MQ=255)
                                         aGTTACCTTTCTCATTATCTCTGGTAACTGGAGATaa                                                              >  1:1540389/1‑37 (MQ=255)
                                           ttACCTTTCTCATTATCTCTGGTAACTGGAGATAAttt                                                           <  1:1835682/38‑1 (MQ=255)
                                            tACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCAt                                            <  1:3293523/52‑1 (MQ=255)
                                                tttCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTg                                       <  1:425748/53‑1 (MQ=255)
                                                  tctcATTATCTCTGGTAACTGGAGATAATTTTTAg                                                       >  1:1581150/1‑35 (MQ=255)
                                                      aTTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgtg                <  1:6439/70‑1 (MQ=255)
                                                      aTTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgt                 <  1:760187/69‑1 (MQ=255)
                                                      aTTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgt                 <  1:2383951/69‑1 (MQ=255)
                                                      aTTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAg                  <  1:987674/68‑1 (MQ=255)
                                                       ttATCTCTGGTAACTGGAGATAATTTTTAGCATaa                                                  <  1:309263/35‑1 (MQ=255)
                                                       ttATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCa                                             <  1:2783323/40‑1 (MQ=255)
                                                       ttATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCa                                             <  1:1353227/40‑1 (MQ=255)
                                                       ttATCTCCAGTTACCAGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgtgt               >  1:81907/1‑70 (MQ=37)
                                                       ttATCTCCAGTTACCAGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgtg                >  1:3003727/1‑69 (MQ=25)
                                                       ttATCTCCAGTTACCAGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgtg                >  1:2140658/1‑69 (MQ=25)
                                                         aTCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATa                                           >  1:3629166/1‑40 (MQ=255)
                                                         aTCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATa                                           >  1:1828557/1‑40 (MQ=255)
                                                         aTCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa             <  1:3043639/70‑1 (MQ=255)
                                                          tctcTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGata           >  1:394953/1‑71 (MQ=255)
                                                            tctGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATAtt         >  1:1316923/1‑71 (MQ=255)
                                                            tctGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATAtt         >  1:827119/1‑71 (MQ=255)
                                                             ctGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAt                                  >  1:1221890/1‑45 (MQ=255)
                                                             ctGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAt                                  >  1:2458626/1‑45 (MQ=255)
                                                              tGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGatat          <  1:3005014/68‑1 (MQ=255)
                                                                gTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa     <  1:2113332/71‑1 (MQ=255)
                                                                  aaCTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAtt   >  1:821376/1‑71 (MQ=255)
                                                                  aaCTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAtt   >  1:1322846/1‑71 (MQ=255)
                                                                   aCTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTg  <  1:1766303/71‑1 (MQ=255)
                                                                      |                                                                   
AAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTG  >  minE/1869775‑1869912

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: