Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,769 A→G 10.9% S361S (AGT→AGC)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,7690AG10.9% 139.7 / 8.9 55S361S (AGT→AGC) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (31/18);  new base G (2/4);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

AACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093700‑2093838
                                                                     |                                                                     
aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                                      <  1:3014455/71‑1 (MQ=255)
aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                                      >  1:513882/1‑71 (MQ=255)
aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                                      <  1:1732695/71‑1 (MQ=255)
 aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                                     >  1:2419385/1‑71 (MQ=255)
 aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                                     >  1:1686890/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)
                    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)
                      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)
                             gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:2008961/1‑60 (MQ=255)
                               tCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgc                                       >  1:2338377/1‑71 (MQ=255)
                                ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                                      >  1:3727431/1‑71 (MQ=255)
                                ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                                      >  1:902015/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:790937/1‑71 (MQ=255)
                                     gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:1352718/1‑71 (MQ=255)
                                     gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:2618907/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:3228860/1‑69 (MQ=255)
                                           tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                             >  1:2260235/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)
                                                    ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                   >  1:3167169/1‑70 (MQ=255)
                                                    atGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacgccacgcc                  <  1:788451/70‑1 (MQ=255)
                                                      gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTg                                           >  1:1979594/1‑44 (MQ=255)
                                                      gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCaa                >  1:3864670/1‑71 (MQ=255)
                                                               cgcCGGACTGCCGCCCAGTATTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt       <  1:1525452/71‑1 (MQ=255)
                                                                    gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAg  <  1:1332275/71‑1 (MQ=255)
                                                                     |                                                                     
AACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093700‑2093838

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