Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,041,438 G→A 67.5% D32N (GAT→AAT)  yqiI → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,041,4380GA67.5% 89.9 / 67.4 83D32N (GAT→AAT) yqiIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/17);  new base A (28/28);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCT  >  minE/2041373‑2041507
                                                                 |                                                                     
tttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATAcc                                                                  >  1:2141042/1‑71 (MQ=255)
tttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAAt                                                                     <  1:4047506/68‑1 (MQ=255)
tttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATAcc                                                                  >  1:1793979/1‑71 (MQ=255)
tttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATAc                                                                   >  1:2315947/1‑70 (MQ=255)
  tttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCg                                                                 <  1:3563194/70‑1 (MQ=255)
    tATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGtt                                                               <  1:3524237/70‑1 (MQ=255)
    tATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGtt                                                               <  1:1775562/70‑1 (MQ=255)
         gTTTTAGTTGATTAGCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCAtt                                                         <  1:1062939/71‑1 (MQ=255)
         gTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCAtt                                                         <  1:3521847/71‑1 (MQ=255)
          ttttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTg                                                        >  1:4844340/1‑71 (MQ=255)
            ttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      >  1:4657081/1‑71 (MQ=255)
            ttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      >  1:5697283/1‑71 (MQ=255)
             tAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGt                                                      <  1:2891877/70‑1 (MQ=255)
             tAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGt                                                      <  1:2857530/70‑1 (MQ=255)
             tAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGt                                                      <  1:5423064/70‑1 (MQ=255)
             tAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      <  1:1373246/70‑1 (MQ=255)
                ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCg                                                                 <  1:2985009/56‑1 (MQ=255)
                ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCg                                                                 <  1:3899831/56‑1 (MQ=255)
                ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      <  1:1584279/67‑1 (MQ=255)
                ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  <  1:277540/71‑1 (MQ=255)
                 tCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      <  1:5251869/66‑1 (MQ=255)
                 tCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  >  1:2151480/1‑70 (MQ=255)
                 tCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  >  1:1926305/1‑70 (MQ=255)
                  cATTACCTGCATGGGCAATGGACTGCTATTCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGt                                                 <  1:881452/70‑1 (MQ=255)
                  cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      >  1:3977535/1‑65 (MQ=255)
                  cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      >  1:5031913/1‑65 (MQ=255)
                  cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGt                                                      >  1:5289081/1‑65 (MQ=255)
                  cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTa                                                >  1:5625310/1‑71 (MQ=255)
                   aTTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGt                                                 <  1:1294255/69‑1 (MQ=255)
                      aCCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTAcc                                              >  1:2132674/1‑69 (MQ=255)
                      aCCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTAc                                               >  1:3332393/1‑68 (MQ=255)
                       ccTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAg                                           <  1:1276411/71‑1 (MQ=255)
                       ccTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAg                                           <  1:1693992/71‑1 (MQ=255)
                       ccTGCATGGGCAATGGACTGCTATGCTGAACATCAAGGTGGAAATACCGTTGTCATTGGTTACGTAcc                                              >  1:4927517/1‑68 (MQ=255)
                        cTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGa                                          <  1:4356270/71‑1 (MQ=255)
                         tGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGaa                                         >  1:2532623/1‑71 (MQ=255)
                          gCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCa                                                           >  1:4558904/1‑52 (MQ=255)
                           cATGGGCAATGGACTGCTATGCTGCACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAAt                                        <  1:4288348/70‑1 (MQ=255)
                           cATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGt                                                      >  1:2008738/1‑56 (MQ=255)
                           cATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTa                                                >  1:5680257/1‑62 (MQ=255)
                           cATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAAtc                                       >  1:1402379/1‑71 (MQ=255)
                           cATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAAt                                        <  1:3027194/70‑1 (MQ=255)
                            aTGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAAtct                                      <  1:4328076/71‑1 (MQ=255)
                             tGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCg                                                                 >  1:2332459/1‑43 (MQ=255)
                              gggCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTa                                                >  1:3046636/1‑59 (MQ=255)
                               ggCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTa                                   >  1:386393/1‑71 (MQ=255)
                                   aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  >  1:2094443/1‑52 (MQ=255)
                                   aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  >  1:3424723/1‑52 (MQ=255)
                                   aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  >  1:3368375/1‑52 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTAGGTACCAAGAATctct                                    >  1:4346173/1‑65 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  <  1:2600859/51‑1 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  <  1:3910435/51‑1 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  <  1:4552632/51‑1 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACg                                                  <  1:3335632/51‑1 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATcccc                              <  1:2980024/71‑1 (MQ=255)
                                    tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATccc                               >  1:4235889/1‑70 (MQ=255)
                                     ggACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAAtctc                                     <  1:1809469/63‑1 (MQ=255)
                                     ggACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATccc                               >  1:3896473/1‑69 (MQ=255)
                                       aCTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                           <  1:5557865/71‑1 (MQ=255)
                                       aCTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                           <  1:3692540/71‑1 (MQ=255)
                                        cTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                           <  1:3920871/70‑1 (MQ=255)
                                          gCTATGCTGAACATGAAGGTGGAGATACCGTTGTCAt                                                          >  1:2204896/1‑37 (MQ=255)
                                           cTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATg                       >  1:3172588/1‑71 (MQ=255)
                                            tATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATg                       <  1:3706266/70‑1 (MQ=255)
                                              tGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTAc                                               <  1:281847/44‑1 (MQ=255)
                                              tGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTa                                                >  1:3923171/1‑43 (MQ=255)
                                              tGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCa                             >  1:2006852/1‑62 (MQ=255)
                                                cTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaa                   <  1:2136776/70‑1 (MQ=255)
                                                 tGAACATGAAGGTGGAGATACCGTTGTCAGTGGTTACGTACCAAGAAtct                                      >  1:5163355/1‑50 (MQ=255)
                                                   aaCATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaaa                 <  1:2969979/69‑1 (MQ=255)
                                                   aaCATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaaa                 <  1:1314020/69‑1 (MQ=255)
                                                      aTGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaaaa                <  1:1281944/67‑1 (MQ=255)
                                                       tGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGaa                                         >  1:1531711/1‑41 (MQ=255)
                                                       tGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                           <  1:2922496/55‑1 (MQ=255)
                                                       tGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCACCAGc                           <  1:4027550/55‑1 (MQ=255)
                                                       tGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATctct                                    <  1:4099404/46‑1 (MQ=255)
                                                        gAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                           <  1:4513307/54‑1 (MQ=255)
                                                        gAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAg              >  1:4064851/1‑67 (MQ=255)
                                                         aaGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGa                         >  1:3931014/1‑55 (MQ=255)
                                                         aaGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGa                         >  1:951085/1‑55 (MQ=255)
                                                          aGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATa        <  1:3735510/71‑1 (MQ=255)
                                                           ggTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATaa       <  1:3882761/71‑1 (MQ=255)
                                                            gTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATaaa      <  1:709657/71‑1 (MQ=255)
                                                               gagaTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAAt    <  1:4560122/70‑1 (MQ=255)
                                                                agaTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCt  >  1:852025/1‑71 (MQ=255)
                                                                 |                                                                     
TTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAGATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCT  >  minE/2041373‑2041507

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: