Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,822 A→G 7.3% I10I (ATT→ATC tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,8220AG 7.3% 228.7 / 7.6 82I10I (ATT→ATCtdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (49/27);  new base G (2/4);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCGGAA  >  minE/2094754‑2094892
                                                                    |                                                                      
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:1677234/1‑71 (MQ=255)
               cTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                                       >  1:1424773/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:2389924/1‑70 (MQ=255)
                  tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                                     >  1:691851/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:2169022/1‑38 (MQ=255)
                                          gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                             >  1:1159900/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:2188270/1‑71 (MQ=255)
                                            aTGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGAtt                          >  1:904506/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:2632604/1‑71 (MQ=255)
                                                  aCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt                    >  1:2097949/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:3468271/37‑1 (MQ=255)
                                                      gAGGGACCAATCCCAATTTTGAAAATATCGAATGCAc                                                  <  1:3471092/37‑1 (MQ=255)
                                                      gAGGGACCAATCCCAATTTTGAAAATATCGAATGCAc                                                  <  1:2924990/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)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCcaca                                    >  1:1604397/1‑39 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt                    >  1:1896096/1‑55 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCa                  >  1:2932591/1‑57 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCgg    <  1:233942/71‑1 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCg     <  1:3521026/70‑1 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCg     >  1:1944920/1‑70 (MQ=255)
                                                                  ccAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCg     >  1:1806156/1‑70 (MQ=255)
                                                                  acAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCa                  >  1:376260/2‑57 (MQ=255)
                                                                   cAATTTTGAAAATATCGAATGCACTAATCATATcc                                       <  1:2460406/35‑1 (MQ=255)
                                                                   cAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCgg    >  1:1523832/1‑70 (MQ=255)
                                                                   cAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCGGa   <  1:2499155/71‑1 (MQ=255)
                                                                    aaTTTTGAAAATATCGAATGCACTAATCATTTCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAg       >  1:1383961/1‑66 (MQ=255)
                                                                    aaTTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCGGaa  >  1:653196/1‑71 (MQ=255)
                                                                    aaTTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCGGaa  >  1:1920585/1‑71 (MQ=255)
                                                                    |                                                                      
AGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGAAGTGGAGCGGAA  >  minE/2094754‑2094892

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: