Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,531,609 A→G 7.3% 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,6090AG 7.3% 140.2 / 4.2 55I14I (ATT→ATCcoaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (2/2);  total (29/26)
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.18e-01

CGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAG  >  minE/2531546‑2531678
                                                               |                                                                     
cGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGAt                                                                >  1:610997/1‑71 (MQ=255)
 gCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATc                                                               >  1:1157188/1‑71 (MQ=255)
   cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACAATTGGTGATGGGATCGa                                                             <  1:205095/71‑1 (MQ=255)
            gTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCg                                                    >  1:603254/1‑71 (MQ=255)
                  gAACATCTGGGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                             <  1:1742911/56‑1 (MQ=255)
                  gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                             <  1:198460/56‑1 (MQ=255)
                   aaCATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATa                                                <  1:1178738/68‑1 (MQ=255)
                     cATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaaa                                              >  1:1978680/1‑68 (MQ=255)
                       tCTGGGTGGCGCGCGTCACGATATCGATATGACCATTGGTaa                                                                      <  1:252961/42‑1 (MQ=255)
                       tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGa                                                             <  1:1666656/51‑1 (MQ=255)
                         tGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCg                                       >  1:160429/1‑71 (MQ=255)
                           cGTGTCGCGAGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCg                                       <  1:2011464/69‑1 (MQ=255)
                            gTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGttt                                    <  1:1036816/71‑1 (MQ=255)
                             tGGCGCGCGTCACGATATCGATATGACCATTGGTa                                                                       <  1:1527863/35‑1 (MQ=255)
                             tGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                             <  1:137778/45‑1 (MQ=255)
                             tGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtttt                                   <  1:1456016/71‑1 (MQ=255)
                               gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTg                                 >  1:2073313/1‑71 (MQ=255)
                                 gcgcgTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAAt                                             >  1:2107791/1‑57 (MQ=255)
                                  cgcgTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCAt                              <  1:1365663/71‑1 (MQ=255)
                                   gcgTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTAcc                                                      >  1:524569/1‑46 (MQ=255)
                                   gcgTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTAcc                                                      >  1:994376/1‑46 (MQ=255)
                                    cgTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATa                             <  1:1844277/70‑1 (MQ=255)
                                      tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCg                                                    <  1:397179/45‑1 (MQ=255)
                                      tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaaca                          <  1:1823148/71‑1 (MQ=255)
                                       cACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGttttt                                  >  1:1091958/1‑62 (MQ=255)
                                       cACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCAt                              >  1:572944/1‑66 (MQ=255)
                                        aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                        >  1:1920187/1‑71 (MQ=255)
                                        aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                        <  1:926144/71‑1 (MQ=255)
                                        aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                        <  1:871002/71‑1 (MQ=255)
                                        aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                        >  1:521012/1‑71 (MQ=255)
                                        aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaa                            >  1:1011184/1‑67 (MQ=255)
                                          gatatcgatatGACCATTGGTAATGGGATCGAAAg                                                          >  1:1659635/1‑35 (MQ=255)
                                           atatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGAt                                                 >  1:624736/1‑43 (MQ=255)
                                            tatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaac                           <  1:5634/64‑1 (MQ=255)
                                            tatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTc                     >  1:768265/1‑70 (MQ=255)
                                            tatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCa                    <  1:2126304/71‑1 (MQ=255)
                                            tatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCa                    <  1:1213639/71‑1 (MQ=255)
                                            tatagatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCt                      >  1:1756603/1‑69 (MQ=255)
                                             atcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCaa                   >  1:2096652/1‑71 (MQ=255)
                                                gatatGAGCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                     <  1:1261719/50‑1 (MQ=255)
                                                gatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCa                    >  1:975820/1‑67 (MQ=255)
                                                gatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGc                <  1:954322/71‑1 (MQ=255)
                                                   atGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGtt             >  1:839521/1‑71 (MQ=255)
                                                    tGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                               <  1:62130/36‑1 (MQ=255)
                                                      aCCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaaca                          >  1:475426/1‑55 (MQ=255)
                                                         aTTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATa                             <  1:903687/49‑1 (MQ=255)
                                                          ttGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc          <  1:435397/67‑1 (MQ=255)
                                                           tGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaa                         >  1:666936/1‑51 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaa                         >  1:465735/1‑50 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaa                         >  1:268945/1‑50 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaa                         >  1:1244307/1‑50 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaa                         >  1:2184536/1‑50 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAAt                  >  1:356305/1‑57 (MQ=255)
                                                            ggTAATGGGATCGAAAGTACCCCGATAAATCGCCCGTTTTTGCATAACAACCTCAAt                  >  1:2062609/1‑57 (MQ=255)
                                                              tAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGttt            <  1:1488684/61‑1 (MQ=255)
                                                              tAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc          <  1:1285584/63‑1 (MQ=255)
                                                              tAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAg  <  1:157506/71‑1 (MQ=255)
                                                               aaTGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCaa                   >  1:694589/1‑53 (MQ=255)
                                                               |                                                                     
CGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAG  >  minE/2531546‑2531678

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: