Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,819 C→G 8.1% G11G (GGG→GGC tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,8190CG 8.1% 216.5 / 6.8 74G11G (GGG→GGCtdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (44/24);  new base G (2/4);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTG  >  minE/2094751‑2094884
                                                                    |                                                                 
tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                                                  >  1:844217/1‑70 (MQ=255)
   aGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                                              <  1:2269486/71‑1 (MQ=255)
     ccGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAttt                                                            <  1:2829431/71‑1 (MQ=255)
          cAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                                                >  1:2769049/1‑62 (MQ=255)
               aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                                                  >  1:1161045/1‑71 (MQ=255)
                 ttgtttaCCCGCTTTCATTGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCg                                                 >  1:716263/4‑70 (MQ=255)
                  cTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:1424773/1‑71 (MQ=255)
                  cTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:1677234/1‑71 (MQ=255)
                   tttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                                                <  1:666668/53‑1 (MQ=255)
                    ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                                                >  1:1399229/1‑52 (MQ=255)
                    ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                             >  1:2781429/1‑71 (MQ=255)
                    ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCACAATTTTGAAAATATCGAATg                                             >  1:2286589/1‑71 (MQ=255)
                     tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCaa                                                               >  1:3692426/1‑52 (MQ=255)
                     tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                             >  1:691851/1‑70 (MQ=255)
                     tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                             >  1:2389924/1‑70 (MQ=255)
                           gctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAAt                                      >  1:2976114/4‑71 (MQ=255)
                            cattcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                             >  1:2490450/1‑63 (MQ=255)
                             tttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGaa                                               <  1:3824543/59‑1 (MQ=255)
                             attcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGa                                                         >  1:1959186/1‑50 (MQ=255)
                              ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                      >  1:324940/1‑52 (MQ=255)
                               tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata                                  >  1:3131267/1‑71 (MQ=255)
                                cattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat                                 >  1:919355/1‑71 (MQ=255)
                                 attGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAtac                                            <  1:2428345/59‑3 (MQ=255)
                                 attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                                 <  1:430453/54‑1 (MQ=255)
                                  ttGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                     <  1:1707496/49‑1 (MQ=255)
                                   tGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAata                                                    >  1:2180987/1‑49 (MQ=255)
                                   tGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCca                              >  1:1812294/1‑71 (MQ=255)
                                      ccccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                      <  1:1680818/44‑1 (MQ=255)
                                       cccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaa                                                       <  1:514929/42‑1 (MQ=255)
                                        ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                        <  1:2947819/40‑1 (MQ=255)
                                           cGGTATGCGAAGTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCcaca                            >  1:203314/1‑65 (MQ=255)
                                           cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                      >  1:2068757/1‑39 (MQ=255)
                                           cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:1134163/1‑46 (MQ=255)
                                           cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCtag                      >  1:807733/1‑69 (MQ=255)
                                           cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCg                      <  1:2638756/71‑1 (MQ=255)
                                           cGGTATGCGAACTGGAGGGACCAATCCAAATTTTGAAAATATAGa                                                >  1:2155893/1‑45 (MQ=38)
                                             gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                     >  1:1159900/1‑38 (MQ=255)
                                             gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                     >  1:2169022/1‑38 (MQ=255)
                                             gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                     >  1:1799083/1‑38 (MQ=255)
                                             gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTa                                        <  1:3131447/51‑1 (MQ=255)
                                             gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGa                    <  1:2363500/71‑1 (MQ=255)
                                              tATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACAccc                         <  1:1408336/65‑1 (MQ=255)
                                              tATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGAt                   <  1:908722/71‑1 (MQ=255)
                                               aTGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGAtt                  >  1:904506/1‑71 (MQ=255)
                                               aTGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGAtt                  >  1:2188270/1‑71 (MQ=255)
                                                 gCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                                >  1:1757316/1‑55 (MQ=255)
                                                  cGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGcc               >  1:2267187/1‑71 (MQ=255)
                                                   gAACTTGAGGGACCAATCCCAATTTTGAAAATATCg                                                 <  1:2000894/36‑1 (MQ=255)
                                                   gAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:2172706/1‑38 (MQ=255)
                                                   gAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCa                                    <  1:939699/49‑1 (MQ=19)
                                                   gAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCa                                    <  1:1218923/49‑1 (MQ=19)
                                                    aaCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata                                  >  1:1754043/1‑50 (MQ=255)
                                                     aCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt            >  1:2097949/1‑71 (MQ=255)
                                                     aCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt            >  1:2632604/1‑71 (MQ=255)
                                                      cTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:2582/1‑35 (MQ=255)
                                                      cTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                               >  1:3097608/1‑35 (MQ=255)
                                                      cTGGAGGGACCAATCCCAATATTGAAAATATCGaa                                               >  1:1620267/1‑35 (MQ=255)
                                                         gAGGGACCAATCCCAATTTTGAAAATATCGAATGCAc                                          <  1:2924990/37‑1 (MQ=255)
                                                         gAGGGACCAATCCCAATTTTGAAAATATCGAATGCAc                                          <  1:3468271/37‑1 (MQ=255)
                                                         gAGGGACCAATCCCAATTTTGAAAATATCGAATGCAc                                          <  1:3471092/37‑1 (MQ=255)
                                                         gAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGa                    >  1:1126004/1‑59 (MQ=255)
                                                          aGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTg       <  1:1981665/71‑1 (MQ=255)
                                                           gggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTg       >  1:284978/1‑70 (MQ=255)
                                                            ggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata                                  <  1:2819680/42‑1 (MQ=255)
                                                            ggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTg                 >  1:1828808/1‑59 (MQ=255)
                                                             gACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAg    <  1:3177156/71‑1 (MQ=255)
                                                              aCCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAg    <  1:3632857/70‑1 (MQ=255)
                                                              aCCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAg    <  1:2163052/70‑1 (MQ=255)
                                                              aCCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAg    <  1:3016498/70‑1 (MQ=255)
                                                              aCCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGt   <  1:3189754/71‑1 (MQ=255)
                                                                cAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                                >  1:3277506/1‑40 (MQ=255)
                                                                cAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTc           >  1:2891123/1‑61 (MQ=255)
                                                                cAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTg  >  1:3613049/1‑70 (MQ=255)
                                                                cAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTg  >  1:1177845/1‑70 (MQ=255)
                                                                   tCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTc                       >  1:3364384/1‑46 (MQ=255)
                                                                    cccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGt             >  1:735780/1‑55 (MQ=255)
                                                                    |                                                                 
TCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTG  >  minE/2094751‑2094884

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: