Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,531,609 A→G 64.2% 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,6090AG64.2% 70.1 / 68.0 67I14I (ATT→ATCcoaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base G (18/25);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

TCGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAGG  >  minE/2531545‑2531679
                                                                |                                                                      
tCGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGa                                                                  >  1:205679/1‑71 (MQ=255)
 cGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGAt                                                                 >  1:3831027/1‑71 (MQ=255)
    cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                              <  1:751677/71‑1 (MQ=255)
    cAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                              <  1:3062264/71‑1 (MQ=255)
           aCGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTAcc                                                       <  1:501099/71‑1 (MQ=255)
           aCGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTAc                                                        <  1:3327675/70‑1 (MQ=255)
             gTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCg                                                     >  1:2033381/1‑71 (MQ=255)
               gATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCgg                                                    <  1:2571571/70‑1 (MQ=255)
                aTCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                              <  1:1565558/59‑1 (MQ=255)
                aTCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                              <  1:532214/59‑1 (MQ=255)
                aTCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGa                                                              <  1:2231608/59‑1 (MQ=255)
                 tCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 <  1:1768260/71‑1 (MQ=255)
                 tCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 <  1:1773287/71‑1 (MQ=255)
                  cGAAGATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATc                                                                >  1:1087639/1‑55 (MQ=255)
                  cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 >  1:2202748/1‑70 (MQ=255)
                  cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 >  1:3176528/1‑70 (MQ=255)
                  cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 >  1:1865324/1‑70 (MQ=255)
                  cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 >  1:1234865/1‑70 (MQ=255)
                  cGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATa                                                 >  1:3866122/1‑70 (MQ=255)
                   gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATaa                                                <  1:2514677/70‑1 (MQ=255)
                   gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATaa                                                <  1:1179245/70‑1 (MQ=255)
                   gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATaa                                                <  1:1740743/70‑1 (MQ=255)
                   gAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATggg                                                                   <  1:2236288/51‑1 (MQ=255)
                    aaCATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATaa                                                <  1:450704/69‑1 (MQ=255)
                     aCATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAAt                                              >  1:3062011/1‑70 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATGTGACCATTGGTGATGGGATCGAAAGTa                                                         <  1:2188996/56‑1 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGcc                                          >  1:682700/1‑71 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGcc                                          >  1:3840741/1‑71 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGcc                                          >  1:3005436/1‑71 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGcc                                          >  1:2939177/1‑71 (MQ=255)
                        tCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGcc                                          >  1:371847/1‑71 (MQ=255)
                         cTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGccc                                         >  1:1396730/1‑71 (MQ=255)
                          tGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCg                                        >  1:619109/1‑71 (MQ=255)
                          atcGTGACGCGCGCCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCg                                        >  1:3579149/3‑71 (MQ=255)
                          atcGTGACGCGCGCCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCg                                        >  1:1712729/3‑71 (MQ=255)
                                gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGttttt                                   >  1:2219614/1‑70 (MQ=255)
                                gcgcgcGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTg                                  >  1:1896799/1‑71 (MQ=255)
                                   cgcgTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCAt                               <  1:2091960/71‑1 (MQ=255)
                                    gcgTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                      <  1:3744273/63‑1 (MQ=255)
                                    gcgTCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                      <  1:1694755/63‑1 (MQ=255)
                                    gcgACACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                      <  1:2004298/63‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGTTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtttt                                    <  1:3917272/62‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaaca                           <  1:1560531/71‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaac                            <  1:375803/70‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaac                            <  1:2569731/70‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATaa                                                <  1:938122/50‑1 (MQ=255)
                                       tCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGtttt                                    <  1:2196123/62‑1 (MQ=255)
                                         aCGATATCGATATGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                         >  1:387747/1‑71 (MQ=255)
                                         aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATaacaac                         >  1:3585132/1‑71 (MQ=255)
                                         aCGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGGTTTTGCATaacaac                         >  1:2078606/1‑71 (MQ=255)
                                            atatcgatatGACCATTGGTAATTGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTc                      <  1:962699/71‑1 (MQ=255)
                                            atatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTc                      <  1:3380911/71‑1 (MQ=255)
                                            atatcgatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTc                      <  1:1883986/71‑1 (MQ=255)
                                               tcgatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGtt                                      <  1:1392505/52‑1 (MQ=255)
                                                cgatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCa                                >  1:4044222/1‑57 (MQ=255)
                                                 gatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCaa                    >  1:1194817/1‑68 (MQ=255)
                                                 gatatGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCaa                    >  1:1602960/1‑68 (MQ=255)
                                                 gatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGc                 <  1:1666520/71‑1 (MQ=255)
                                                 gatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGc                 >  1:490565/1‑71 (MQ=255)
                                                 gatatGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGc                 >  1:1885637/1‑71 (MQ=255)
                                                     tGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCAt                               >  1:663493/1‑53 (MQ=255)
                                                     tGACCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCAt                               >  1:2593069/1‑53 (MQ=255)
                                                       aCCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc           <  1:1721973/71‑1 (MQ=255)
                                                       aCCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc           <  1:594714/71‑1 (MQ=255)
                                                       aCCATTGGTGATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc           <  1:788481/71‑1 (MQ=255)
                                                       aCCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTc           <  1:1818370/71‑1 (MQ=255)
                                                              gTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTTCATAACAACCTCAATGCGTTTTCGGTGGc     <  1:2386384/70‑1 (MQ=255)
                                                                aaTGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAgg  >  1:2855933/1‑71 (MQ=255)
                                                                |                                                                      
TCGCCAGAATAACGTGATCGAACATCTGCGTGGCGCGCGTCACGATATCGATATGACCATTGGTAATGGGATCGAAAGTACCCGGATAAATCGCCCGTTTTTGCATAACAACCTCAATGCGTTTTCGGTGGCAGG  >  minE/2531545‑2531679

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: