Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 890,121 T→C 28.6% intergenic (‑128/+240) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE890,1210TC28.6% 57.1 / 27.3 42intergenic (‑128/+240)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base C (4/8);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGG  >  minE/890052‑890185
                                                                     |                                                                
tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGt                                                                  <  1:238150/70‑1 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTcc                                                                >  1:3929926/1‑71 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGt                                                                  <  1:1486473/67‑1 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGAc                                                             >  1:2292756/1‑71 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCg                                                           <  1:3384936/71‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCg                                                           <  1:715483/71‑1 (MQ=255)
       cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGAcctg                                                          <  1:3761446/71‑3 (MQ=255)
       cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCg                                                           >  1:560078/1‑70 (MQ=255)
       cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCgg                                                          <  1:2731966/71‑1 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAg                                                        <  1:1424921/71‑1 (MQ=255)
            cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAg                                                     >  1:1531677/1‑71 (MQ=255)
                cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTg                                                 <  1:2065692/71‑1 (MQ=255)
                    tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGc                                             <  1:1414533/71‑1 (MQ=255)
                     aaaCGGTTATTATAATCTACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGa                                                              <  1:379112/53‑1 (MQ=255)
                      aaCGGTTATTATAATCACCCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGAcc                                                            >  1:3425928/1‑54 (MQ=255)
                      aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAg                                           >  1:2716026/1‑71 (MQ=255)
                       aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgg                                          <  1:2411881/71‑1 (MQ=255)
                        cGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAggg                                         <  1:2514551/71‑1 (MQ=255)
                         ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAgggg                                        <  1:3933622/71‑1 (MQ=255)
                          gTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAggggg                                       <  1:290015/71‑1 (MQ=255)
                           ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggggg                                      <  1:2677782/71‑1 (MQ=255)
                           ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggggg                                      <  1:507280/71‑1 (MQ=255)
                           ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAgggggg                                      <  1:2974543/71‑1 (MQ=255)
                              ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCa                                            >  1:2494361/1‑62 (MQ=255)
                                 taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCtt                                >  1:2384180/1‑71 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCtttt                              <  1:2031373/71‑1 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCtttt                              <  1:2387766/71‑1 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCttt                               <  1:278093/70‑1 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCttt                               <  1:3736833/70‑1 (MQ=255)
                                    tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggg                                        >  1:3471399/1‑60 (MQ=255)
                                    tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggg                                        >  1:1880953/1‑60 (MQ=255)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTc                        <  1:3131148/71‑1 (MQ=255)
                                             cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCgg                                              <  1:494759/45‑1 (MQ=255)
                                              tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCtt                                <  1:3138295/58‑1 (MQ=255)
                                                tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCt                   <  1:2850738/69‑1 (MQ=255)
                                                   gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAAccc              <  1:2096936/71‑1 (MQ=255)
                                                    tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAggg                                         <  1:3166798/43‑1 (MQ=255)
                                                      gCTATAGTTGCCAGGCCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCctc            <  1:3553722/70‑1 (MQ=255)
                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggg                                        <  1:1583862/41‑1 (MQ=255)
                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggg                                        <  1:3863555/41‑1 (MQ=255)
                                                        tataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAggggggg                                     <  1:2637437/43‑1 (MQ=255)
                                                             ttGCCAGGCCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAAc                >  1:1449955/1‑59 (MQ=255)
                                                                ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg  <  1:3093127/70‑1 (MQ=255)
                                                                ccAGGCCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg  <  1:3353720/70‑1 (MQ=255)
                                                                     |                                                                
TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGG  >  minE/890052‑890185

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: