Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,000 T→C 21.0% A284A (GCA→GCG tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,0000TC21.0% 149.0 / 36.1 76A284A (GCA→GCGtdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base C (8/8);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

TACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCA  >  minE/2093931‑2094065
                                                                     |                                                                 
tACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGt                                                                   <  1:746612/70‑1 (MQ=255)
  ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                                                <  1:1854728/71‑1 (MQ=255)
  ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                                                <  1:3129887/71‑1 (MQ=255)
   cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                                               <  1:1144772/71‑1 (MQ=255)
   cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                                               >  1:2878945/1‑71 (MQ=255)
   cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                                               >  1:893794/1‑71 (MQ=255)
   cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                                               >  1:3357418/1‑71 (MQ=255)
    ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                              <  1:989600/71‑1 (MQ=255)
     gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                              >  1:3551870/1‑70 (MQ=255)
     gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                              >  1:409085/1‑70 (MQ=255)
       cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                                           <  1:3261656/71‑1 (MQ=255)
       cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                                           <  1:1173665/71‑1 (MQ=255)
        gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                            <  1:3471859/69‑1 (MQ=255)
        gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                          <  1:3532993/71‑1 (MQ=255)
        gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                          <  1:84647/71‑1 (MQ=255)
        gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                          <  1:1039982/71‑1 (MQ=255)
           taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                        <  1:711265/68‑1 (MQ=255)
            gAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                      <  1:186443/71‑1 (MQ=255)
             aaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCAcgc                                                      <  1:163734/70‑1 (MQ=255)
             aaCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgcg                                                     >  1:99484/1‑71 (MQ=255)
             aaCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgcg                                                     >  1:279286/1‑71 (MQ=255)
               cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                    >  1:389199/1‑70 (MQ=255)
               cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                   >  1:2515975/1‑71 (MQ=255)
               cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                   >  1:2520396/1‑71 (MQ=255)
                ttATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg                                                       <  1:3082986/66‑1 (MQ=255)
                 tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgaccg                                                 >  1:3615862/3‑71 (MQ=255)
                 tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                  >  1:3278248/3‑70 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:1591342/1‑71 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:73143/1‑71 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:1539540/1‑71 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:1624712/1‑71 (MQ=255)
                   tCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgcc                                               <  1:2693478/71‑1 (MQ=255)
                   tCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                <  1:2267458/70‑1 (MQ=255)
                    cGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:1069282/1‑69 (MQ=255)
                    cGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                >  1:3070272/1‑69 (MQ=255)
                     gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                            >  1:1458611/1‑56 (MQ=255)
                     gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAg                                             >  1:2664009/1‑71 (MQ=255)
                     gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAg                                             <  1:885111/71‑1 (MQ=255)
                      taataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGc                                            >  1:2931244/1‑71 (MQ=255)
                       aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgccc                                                  >  1:103838/1‑64 (MQ=255)
                           cgCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGt                                       <  1:878775/70‑1 (MQ=255)
                            gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                                      <  1:3639644/71‑1 (MQ=255)
                            gCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCTCAGTCACCACGCGACCGCCAGCAGCGtt                                      <  1:1300937/71‑1 (MQ=255)
                             ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                     <  1:344529/71‑2 (MQ=255)
                              cAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCAcgcgac                                                   >  1:1228213/1‑54 (MQ=255)
                              cAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttct                                    >  1:3765492/1‑69 (MQ=255)
                                gaaCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAg                                             >  1:636282/3‑60 (MQ=255)
                                gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                  >  1:3707312/1‑71 (MQ=255)
                                gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                  >  1:2667918/1‑71 (MQ=255)
                                 tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                 <  1:3127909/71‑1 (MQ=255)
                                 aaCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctctt                                 >  1:3484239/2‑71 (MQ=255)
                                  aCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                  >  1:3847039/1‑53 (MQ=255)
                                   cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                               <  1:1568880/71‑2 (MQ=255)
                                   cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttcttc                                >  1:1453895/1‑70 (MQ=255)
                                     gCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTtctc                                   <  1:346085/65‑1 (MQ=255)
                                     gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                 >  1:1656806/1‑67 (MQ=255)
                                      ccGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGtt                                      >  1:986529/1‑61 (MQ=255)
                                       cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                           <  1:421292/71‑1 (MQ=255)
                                       cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                           <  1:2554477/71‑1 (MQ=255)
                                        ggAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAc                          >  1:3107302/1‑71 (MQ=255)
                                           ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgcc                                               >  1:3078741/1‑47 (MQ=255)
                                                gCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                       <  1:1832254/66‑1 (MQ=255)
                                                   gcacgcacCGTTAGTCGGTGCGGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                >  1:777986/1‑70 (MQ=255)
                                                   gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAgcgc                 >  1:1339390/1‑69 (MQ=255)
                                                    cacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGt               >  1:2436712/1‑70 (MQ=255)
                                                       gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                       <  1:1597629/59‑1 (MQ=255)
                                                       gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAgcg                  <  1:1285285/64‑1 (MQ=255)
                                                       gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACAt           >  1:1447402/1‑71 (MQ=255)
                                                        cacCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAGTGCAAACATg          <  1:1704461/71‑1 (MQ=255)
                                                           cGTTAGTCGGTGCCGTGACCACGCGCCCGCCAGCTGCGTTttct                                  >  1:56041/1‑44 (MQ=255)
                                                            gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACt                             <  1:1731682/48‑1 (MQ=255)
                                                            gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCTCCAGCGCGTACATGTTGa      <  1:3195049/71‑1 (MQ=255)
                                                              tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc    >  1:3610100/1‑71 (MQ=255)
                                                               aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc    <  1:745603/70‑1 (MQ=255)
                                                               aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc    <  1:2531017/70‑1 (MQ=255)
                                                                gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa  <  1:1597341/71‑1 (MQ=255)
                                                                 tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGtt        <  1:844842/64‑1 (MQ=255)
                                                                     |                                                                 
TACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCA  >  minE/2093931‑2094065

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: