Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,531,609 A→G 24.0% I14I (ATT→ATC coaD ← pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,531,6090AG24.0% 138.3 / 47.3 75I14I (ATT→ATCcoaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/33);  new base G (2/16);  total (26/49)
Fisher's exact test for biased strand distribution p-value = 2.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCA  >  minE/2531547‑2531677
                                                              |                                                                    
gacAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATc                                                              >  1:1902855/3‑71 (MQ=255)
gCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGTGATc                                                              >  1:1556255/1‑71 (MQ=255)
gCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGAt                                                               >  1:1672744/1‑70 (MQ=255)
gCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATc                                                              >  1:3415800/1‑71 (MQ=255)
gCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATc                                                              >  1:1561243/1‑71 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:2411482/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:2371389/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:3392291/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:3809034/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:3423360/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:3438025/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:3757157/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:278756/71‑1 (MQ=255)
  cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:3121255/71‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGaa                                                           <  1:3705385/70‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGaa                                                           <  1:2254049/70‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGaa                                                           <  1:3078958/70‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGaa                                                           <  1:536904/70‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:1640663/69‑1 (MQ=255)
    gAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:3196299/69‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAg                                                         <  1:951877/71‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAg                                                         <  1:3808997/71‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAg                                                         <  1:383703/71‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAg                                                         <  1:2138914/71‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAg                                                         <  1:332216/71‑1 (MQ=255)
     aaTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAg                                                         <  1:448639/71‑1 (MQ=255)
      aTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGaaa                                                          >  1:3798511/1‑69 (MQ=255)
      aTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGaaa                                                          >  1:3106751/1‑69 (MQ=255)
      aTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGaaa                                                          >  1:1122306/1‑69 (MQ=255)
      aTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGaaa                                                          >  1:1941679/1‑69 (MQ=255)
      aTAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGt                                                        <  1:550391/71‑1 (MQ=255)
       tAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                            <  1:1593202/66‑1 (MQ=255)
         aCGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:4005919/64‑1 (MQ=255)
                cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTaa                                                                     >  1:503985/1‑48 (MQ=255)
                cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGAATCGAAAGTACCCGGATa                                               >  1:4053039/1‑70 (MQ=255)
                 gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATaa                                              >  1:439961/1‑70 (MQ=255)
                    cATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                            >  1:253011/1‑69 (MQ=255)
                    cATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                            >  1:1403806/1‑69 (MQ=255)
                    cATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                            >  1:3151919/1‑69 (MQ=255)
                      tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:1516068/51‑1 (MQ=255)
                      tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                            <  1:7840/51‑1 (MQ=255)
                        tGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCg                                      >  1:2258684/1‑71 (MQ=255)
                        tGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCg                                      >  1:1002033/1‑71 (MQ=255)
                        tGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCg                                      >  1:3993322/1‑71 (MQ=255)
                          cGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                    <  1:3448548/71‑1 (MQ=255)
                            tGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtttt                                  <  1:3427920/71‑1 (MQ=255)
                             ggCGCGCGTCACGATATCGATATGACCATTGGTaa                                                                     >  1:3701942/1‑35 (MQ=255)
                              gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTg                                <  1:2100634/71‑1 (MQ=255)
                              gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTg                                <  1:874319/71‑1 (MQ=255)
                              gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTg                                <  1:3537027/71‑1 (MQ=255)
                                   cgTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATa                            >  1:236997/1‑70 (MQ=255)
                                     tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                              <  1:2729264/50‑1 (MQ=255)
                                     tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                              <  1:2189321/50‑1 (MQ=255)
                                     tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                              <  1:264368/50‑1 (MQ=255)
                                     tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                              <  1:3998121/50‑1 (MQ=255)
                                      cACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                              <  1:2349042/49‑1 (MQ=255)
                                        cgatatcgatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAAcc                      <  1:1432741/71‑1 (MQ=255)
                                        cgatatcgatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAAcc                      <  1:2218775/71‑1 (MQ=255)
                                         gatatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCt                     <  1:2546353/71‑1 (MQ=255)
                                         gatatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCt                     <  1:2119231/71‑1 (MQ=255)
                                         gatatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCt                     <  1:3144935/71‑1 (MQ=255)
                                                 tatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGt             <  1:2734720/71‑1 (MQ=255)
                                                     aCCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                            >  1:1286116/1‑36 (MQ=255)
                                                     aCCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                            >  1:2555041/1‑36 (MQ=255)
                                                     aCCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGttt                                   >  1:1652642/1‑45 (MQ=255)
                                                     aCCATTGGTAATGGGATCGAAAATACCCGGATAAAt                                            >  1:3289777/1‑36 (MQ=255)
                                                       cATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCgg       <  1:3842631/71‑1 (MQ=255)
                                                         ttGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTg     <  1:2802686/71‑1 (MQ=255)
                                                         ttGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTg     <  1:2468972/71‑1 (MQ=255)
                                                           ggTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTgg    >  1:2020792/1‑70 (MQ=255)
                                                           ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATa                            >  1:2979884/1‑46 (MQ=255)
                                                           ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATa                            >  1:3554480/1‑46 (MQ=255)
                                                           ggTAATGGGATCGAAAGTACCCGGATAAATCGCCAGTTTTTGCATa                            >  1:984185/1‑46 (MQ=255)
                                                            gTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCa  <  1:2721507/71‑1 (MQ=255)
                                                            gTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCa  <  1:1204361/71‑1 (MQ=255)
                                                             tAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCa                              <  1:555529/42‑1 (MQ=255)
                                                             tAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCa                              <  1:1776618/42‑1 (MQ=255)
                                                              |                                                                    
GCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCA  >  minE/2531547‑2531677

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: