Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,751 A→G 15.7% N367N (AAT→AAC) ‡ tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,7510AG15.7% 164.9 / 22.4 70N367N (AAT→AAC) ‡tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (2/9);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

ATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCA  >  minE/2093683‑2093818
                                                                    |                                                                   
aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                                   <  1:755610/71‑1 (MQ=255)
aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                                   >  1:1312345/1‑71 (MQ=255)
aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTACGCCGCAtt                                                                   <  1:3689283/71‑1 (MQ=255)
  ttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                                  <  1:2937035/70‑1 (MQ=255)
   tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                                  <  1:2619489/69‑1 (MQ=255)
     tACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCac                                                               <  1:2748418/70‑2 (MQ=255)
     tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                               <  1:660514/70‑1 (MQ=255)
     tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                               <  1:665244/70‑1 (MQ=255)
     tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                               <  1:1846827/70‑1 (MQ=255)
       ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                                            <  1:3773520/71‑1 (MQ=255)
       ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                                            <  1:621551/71‑1 (MQ=255)
        cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                                           <  1:3419973/71‑1 (MQ=255)
        cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                                           <  1:3625457/71‑1 (MQ=255)
           tCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTg                                                          <  1:3655475/69‑1 (MQ=9)
             cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                       <  1:2887763/70‑1 (MQ=255)
               gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCAATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                                     >  1:3141530/1‑70 (MQ=255)
                cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                                   <  1:1670782/69‑2 (MQ=255)
                cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                                   <  1:2687725/69‑2 (MQ=255)
                 aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                  <  1:1732695/71‑1 (MQ=255)
                 aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                  >  1:513882/1‑71 (MQ=255)
                 aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                  <  1:3014455/71‑1 (MQ=255)
                  aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                 >  1:1686890/1‑71 (MQ=255)
                  aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                 >  1:2419385/1‑71 (MQ=255)
                   ccGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCtt                                                >  1:2433961/1‑70 (MQ=255)
                    cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgc                                               >  1:613875/1‑71 (MQ=255)
                          gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGccc                                          <  1:838594/70‑1 (MQ=255)
                           cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                                        >  1:2013484/1‑71 (MQ=255)
                               tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt                                    <  1:103575/71‑1 (MQ=255)
                                 aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTTCATACCTGCGCCGGACTGCCGCCCAGTAGTTc                                  <  1:1713055/71‑1 (MQ=255)
                                 aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                   >  1:328614/1‑70 (MQ=255)
                                  gCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                                 <  1:1277156/71‑1 (MQ=255)
                                   cccAAGGTTATGGTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                       >  1:2891604/1‑48 (MQ=255)
                                     cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                                 <  1:3106494/68‑1 (MQ=255)
                                     cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                              <  1:2743762/71‑1 (MQ=255)
                                     cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGt                               >  1:1230071/1‑70 (MQ=255)
                                      aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagt                                       <  1:3016703/61‑1 (MQ=255)
                                       aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                       >  1:833009/1‑44 (MQ=255)
                                       aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                            >  1:143060/1‑71 (MQ=255)
                                       aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                             >  1:76089/1‑70 (MQ=255)
                                        ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAc                           <  1:265030/71‑1 (MQ=255)
                                         gTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGtt                                   <  1:3850049/62‑1 (MQ=255)
                                         gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTg                                                >  1:1920682/1‑49 (MQ=255)
                                          ttATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAccc                         <  1:306517/71‑1 (MQ=255)
                                            aTGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCt                        <  1:2570331/70‑1 (MQ=255)
                                             tGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTg                       <  1:995383/70‑1 (MQ=255)
                                              gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                                                      <  1:2275269/38‑1 (MQ=255)
                                              gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                >  1:2008961/1‑60 (MQ=255)
                                                tCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgc                   >  1:2338377/1‑71 (MQ=255)
                                                 ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                  >  1:902015/1‑71 (MQ=255)
                                                 ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                  >  1:3727431/1‑71 (MQ=255)
                                                  cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                 >  1:115773/1‑71 (MQ=255)
                                                   aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgcc                     >  1:1887974/1‑66 (MQ=255)
                                                    tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgcc               <  1:2749942/71‑1 (MQ=255)
                                                      gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCa              <  1:3865026/70‑1 (MQ=255)
                                                      gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt             >  1:1352718/1‑71 (MQ=255)
                                                      gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt             >  1:2618907/1‑71 (MQ=255)
                                                      gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt             >  1:790937/1‑71 (MQ=255)
                                                       cgATTTCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCAtt            <  1:31148/71‑1 (MQ=255)
                                                        gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                   >  1:869469/1‑47 (MQ=255)
                                                         aTTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTg           >  1:1832840/1‑70 (MQ=255)
                                                          tttCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCAt             <  1:546346/67‑1 (MQ=255)
                                                           ttCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGag         >  1:3245122/1‑69 (MQ=255)
                                                           ttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                              <  1:3523834/49‑1 (MQ=255)
                                                           ttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                            >  1:3488973/1‑51 (MQ=255)
                                                            tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                            >  1:319873/1‑50 (MQ=255)
                                                            tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa         >  1:2260235/1‑69 (MQ=255)
                                                            tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa         >  1:3228860/1‑69 (MQ=255)
                                                                ccgcATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                   >  1:2083737/1‑39 (MQ=255)
                                                                ccgcATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                   >  1:2461100/1‑39 (MQ=255)
                                                                 cgcATTGCATACCTGCGCCGGACTTCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCa  >  1:1650097/1‑71 (MQ=255)
                                                                    |                                                                   
ATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCA  >  minE/2093683‑2093818

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: ‑