Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,511,849 T→C 58.5% W373R (TGG→CGG)  pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,511,8490TC58.5% 23.7 / 47.3 41W373R (TGG→CGG) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base C (14/10);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGT  >  minE/1511780‑1511911
                                                                     |                                                              
gTTCCGGTTGACGATCACGAACGTATCAAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTg                                                               <  1:1134539/71‑1 (MQ=255)
    cGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATc                                                           <  1:905826/71‑1 (MQ=255)
    cGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATc                                                           <  1:357536/71‑1 (MQ=255)
      gTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGa                                                         <  1:989728/71‑1 (MQ=255)
                  gAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGc                                             >  1:1108751/1‑71 (MQ=255)
                        aTCAAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTc                                       >  1:1515521/1‑71 (MQ=255)
                         tCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCg                                                          <  1:1108309/51‑1 (MQ=255)
                           gagaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGaa                                                        >  1:1764828/1‑51 (MQ=255)
                           gagaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGaa                                                        >  1:390607/1‑51 (MQ=255)
                           aagaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAtc                                                      >  1:567096/2‑53 (MQ=255)
                             gaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGa                                                         >  1:1622170/1‑48 (MQ=255)
                             gaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGa                                                 >  1:345036/1‑56 (MQ=255)
                                aGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAgcg                                 >  1:1301573/1‑69 (MQ=255)
                                    cAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGa                                    >  1:615446/1‑62 (MQ=255)
                                    cAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGcc                            <  1:1049995/70‑1 (MQ=255)
                                    cAGGAATACGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGcc                            <  1:340664/70‑1 (MQ=255)
                                      ggAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCcgtc                         <  1:1242886/71‑1 (MQ=255)
                                         aTACGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCgtcgtc                      >  1:397299/1‑71 (MQ=255)
                                          tACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCAc                                          <  1:339185/50‑1 (MQ=255)
                                          tACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCAc                                          <  1:244260/50‑1 (MQ=255)
                                          tACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCAc                                          <  1:150249/50‑1 (MQ=255)
                                           aCGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtct                     >  1:1799430/1‑70 (MQ=255)
                                            cGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtctg                    <  1:676223/70‑1 (MQ=255)
                                            cGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGGCACTTCTGAGCGCAGc                             >  1:1663248/1‑61 (MQ=255)
                                            cGTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtctgt                   >  1:1435808/1‑71 (MQ=255)
                                             gTTGCTAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtctgtc                  <  1:1762847/71‑1 (MQ=255)
                                               tGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCa                                           >  1:1683849/1‑44 (MQ=255)
                                                  tAACTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCgtcgtc                      <  1:1468889/62‑1 (MQ=255)
                                                    aCTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGc             <  1:402149/69‑1 (MQ=255)
                                                    aCTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCt           >  1:492348/1‑71 (MQ=255)
                                                    aCTACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCt           >  1:621203/1‑71 (MQ=255)
                                                      tACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTg                                     >  1:461450/1‑43 (MQ=255)
                                                      tACATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGc         >  1:1556781/1‑71 (MQ=255)
                                                       aCATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAg                                   >  1:1020126/1‑44 (MQ=255)
                                                        cATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCcgtc                         <  1:775474/53‑1 (MQ=255)
                                                        cATCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCg        <  1:1468181/70‑1 (MQ=255)
                                                          tCAACGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtct                     <  1:1859303/55‑1 (MQ=255)
                                                            aaCGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtctgt                   >  1:107310/1‑55 (MQ=255)
                                                            aaCGCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCg   <  1:1376931/71‑1 (MQ=255)
                                                             aCGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGt  <  1:1317038/71‑1 (MQ=255)
                                                               gCTGACCGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCgtct                     >  1:1859785/1‑50 (MQ=255)
                                                                     |                                                              
GTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGT  >  minE/1511780‑1511911

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: