Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 62,472 T→G 5.2% A329A (GCA→GCC)  polB ← DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE62,4720TG 5.2% 217.1 / 6.4 77A329A (GCA→GCC) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (48/25);  new base G (2/2);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

AATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGA  >  minE/62407‑62539
                                                                 |                                                                   
aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCa                                                                 >  1:5404202/1‑70 (MQ=255)
aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:2615766/1‑71 (MQ=255)
aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:2338859/1‑71 (MQ=255)
aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:6787259/1‑71 (MQ=255)
  tGGCATGATTTCATTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:5620236/71‑1 (MQ=255)
  tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:2180690/71‑1 (MQ=255)
  tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:6758609/71‑1 (MQ=255)
  tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:694694/71‑1 (MQ=255)
       tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         >  1:1933420/1‑71 (MQ=255)
       tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:4466893/71‑1 (MQ=255)
       tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:848712/71‑1 (MQ=255)
       tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         >  1:1031841/1‑71 (MQ=255)
       tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:252925/71‑1 (MQ=255)
        gATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGtt                                                        >  1:5307856/1‑71 (MQ=255)
              aGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCt                                                   <  1:6394136/70‑1 (MQ=255)
                  ttGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:2461105/1‑71 (MQ=255)
                   tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:1799094/1‑63 (MQ=255)
                   tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:808554/1‑70 (MQ=255)
                   tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              <  1:4014366/70‑1 (MQ=255)
                     tgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            <  1:2167280/70‑1 (MQ=255)
                      ggAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            >  1:2489917/1‑69 (MQ=255)
                       gAAGATGCGCGTGACCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGaaca                                         >  1:241191/1‑69 (MQ=255)
                          gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGaaa                                          <  1:5216757/67‑1 (MQ=255)
                          gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggc                                      >  1:6362388/1‑71 (MQ=255)
                          gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggc                                      >  1:1253820/1‑71 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:1082083/1‑54 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:5266919/1‑54 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:1019963/1‑71 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:3210811/1‑71 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:5958180/1‑71 (MQ=255)
                            tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:6150460/1‑71 (MQ=255)
                             cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCCGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:4736683/1‑71 (MQ=255)
                             cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:607650/1‑71 (MQ=255)
                             cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:932236/1‑71 (MQ=255)
                             cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCCGt                                   >  1:6365210/1‑71 (MQ=255)
                                cGTGACCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTcga                                <  1:1547846/71‑3 (MQ=255)
                                    aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:3229676/1‑54 (MQ=255)
                                    aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:4784393/1‑54 (MQ=255)
                                    aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:4776020/1‑54 (MQ=255)
                                    aCCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTCGAtctc                            >  1:2422970/1‑71 (MQ=255)
                                     ccAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            <  1:2914662/54‑1 (MQ=255)
                                     ccAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTCGATCTCg                           <  1:2096503/71‑1 (MQ=255)
                                       aGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:2484461/1‑50 (MQ=255)
                                        gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATc                                                    <  1:191269/43‑1 (MQ=255)
                                        gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:6141302/1‑61 (MQ=255)
                                        gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:3597422/1‑61 (MQ=255)
                                        gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:6314886/1‑61 (MQ=255)
                                        gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTc                         >  1:5993921/1‑70 (MQ=255)
                                         cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGTCGGTCAATTTCGTCCa                       <  1:2241917/71‑1 (MQ=255)
                                         cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCTGTCAATTTCGTCCa                       <  1:2330194/71‑1 (MQ=255)
                                         cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCa                       <  1:3561457/71‑1 (MQ=255)
                                           cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGaaa                                          >  1:4517680/1‑50 (MQ=255)
                                           cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCaa                                >  1:5567953/1‑60 (MQ=255)
                                           cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCaa                                >  1:6584173/1‑60 (MQ=255)
                                            gCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      <  1:5604015/69‑1 (MQ=255)
                                             cAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCg                   <  1:23250/71‑1 (MQ=255)
                                               aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcg                                        >  1:3756049/1‑48 (MQ=255)
                                               aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTcc                        >  1:5704933/1‑64 (MQ=255)
                                               aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTcc                        >  1:601340/1‑64 (MQ=255)
                                                tCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATc                >  1:5338803/1‑71 (MQ=255)
                                                  tttCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTc                         <  1:6052000/60‑1 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcg                                        >  1:2513290/1‑43 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      >  1:4784659/1‑61 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      >  1:6651069/1‑61 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATcc               <  1:5823581/68‑1 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCa             <  1:2528479/70‑1 (MQ=255)
                                                    tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCa             <  1:2566889/70‑1 (MQ=255)
                                                       ggTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACggg         <  1:198289/71‑1 (MQ=255)
                                                           ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      <  1:2125439/70‑1 (MQ=255)
                                                           ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      >  1:3036434/1‑70 (MQ=255)
                                                           ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      >  1:178466/1‑70 (MQ=255)
                                                           ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTAt     >  1:1854288/1‑71 (MQ=255)
                                                           ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTAt     >  1:2905641/1‑71 (MQ=255)
                                                            taAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    <  1:5432507/71‑1 (MQ=255)
                                                             aaGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    >  1:2200664/1‑70 (MQ=255)
                                                             aaGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    >  1:1007033/1‑70 (MQ=255)
                                                              aGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGa  >  1:6806464/1‑71 (MQ=255)
                                                                 |                                                                   
AATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGA  >  minE/62407‑62539

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: ‑