Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,255,649 G→A 17.6% A113A (GCG→GCA yecE → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,255,6490GA17.6% 116.2 / 16.7 51A113A (GCG→GCAyecEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/27);  new base A (3/6);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTT  >  minE/1255581‑1255719
                                                                    |                                                                      
ttGGCTCCGCGCATTGGTCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCt                                                                      <  1:1357052/71‑1 (MQ=255)
  ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCttt                                                                    >  1:1694014/1‑71 (MQ=255)
  ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCttt                                                                    >  1:2882525/1‑71 (MQ=255)
  ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt                                                                     <  1:2893194/70‑1 (MQ=255)
  ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt                                                                     <  1:3923784/70‑1 (MQ=255)
  ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt                                                                     <  1:2577941/70‑1 (MQ=255)
    cTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTg                                                                   >  1:2270399/1‑70 (MQ=255)
           cATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtc                                                           <  1:1581803/71‑1 (MQ=255)
             ttGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc                                                         <  1:274454/71‑1 (MQ=255)
                  cAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtc                                                    <  1:2282760/71‑1 (MQ=255)
                   aaTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtct                                                   >  1:2282387/1‑71 (MQ=255)
                    aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAttt                                                             <  1:1014249/60‑1 (MQ=255)
                    aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc                                                  >  1:3716368/1‑71 (MQ=255)
                    aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc                                                  >  1:3745330/1‑71 (MQ=255)
                     tACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct                                                 <  1:789457/71‑1 (MQ=255)
                     tACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct                                                 <  1:2739988/71‑1 (MQ=255)
                       cTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTc                                               <  1:1274923/71‑1 (MQ=255)
                       cTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCg                                                        >  1:3328318/1‑62 (MQ=255)
                        tGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc                                                         <  1:1511582/60‑1 (MQ=255)
                          gCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGa                                                       >  1:3441020/1‑60 (MQ=255)
                          gCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtct                                                   >  1:59007/1‑64 (MQ=255)
                            tGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc                                                         <  1:1032105/56‑1 (MQ=255)
                            tGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTg                                            <  1:2447934/69‑1 (MQ=255)
                             gCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc                                                  >  1:2165122/1‑62 (MQ=255)
                             gCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACtt                                                                     <  1:617042/43‑1 (MQ=255)
                              cAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGa                                        <  1:806358/71‑1 (MQ=255)
                               aaCTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGaa                                       >  1:482184/1‑71 (MQ=255)
                                    cctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtct                                                          >  1:3874622/1‑47 (MQ=255)
                                    cctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTa                                   <  1:734921/70‑1 (MQ=255)
                                      tgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt                                    >  1:1180285/1‑67 (MQ=255)
                                       gccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTa                               >  1:852514/1‑71 (MQ=255)
                                        ccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAtttac                                                           <  1:246899/42‑3 (MQ=255)
                                        ccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTTGATTCTCTTCCTGGCGAATTTAATTa                               <  1:3889856/70‑1 (MQ=255)
                                        ccACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAtt                                <  1:3209061/69‑1 (MQ=255)
                                          acaTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATgg                            >  1:1936153/1‑71 (MQ=255)
                                            aTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCg                          <  1:1036221/71‑1 (MQ=255)
                                            aTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCg                          <  1:619601/71‑1 (MQ=255)
                                             ttCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTg                                            <  1:1631846/52‑1 (MQ=255)
                                             ttCGGCCACGGGAGCAGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTa                                   >  1:3249420/1‑60 (MQ=255)
                                                ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg                       <  1:3198055/70‑1 (MQ=255)
                                                ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg                       <  1:368080/70‑1 (MQ=255)
                                                ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg                       <  1:1984661/70‑1 (MQ=255)
                                                   ccACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct                                                 <  1:1502424/41‑1 (MQ=255)
                                                    cACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt                                    <  1:3150780/53‑1 (MQ=255)
                                                    cACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt                                    <  1:3621757/53‑1 (MQ=255)
                                                     aCGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTcc                 >  1:763412/1‑71 (MQ=255)
                                                      cGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAtt                                <  1:1903074/55‑1 (MQ=255)
                                                      cGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCg                >  1:3720488/1‑71 (MQ=255)
                                                       gggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAt                                 <  1:1637750/53‑1 (MQ=255)
                                                       gggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGc               <  1:1480323/71‑1 (MQ=255)
                                                           gCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATc           <  1:3272407/71‑1 (MQ=255)
                                                                  gcgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtt    >  1:1371676/1‑71 (MQ=255)
                                                                   cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTaa                                  >  1:570281/1‑40 (MQ=255)
                                                                   cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGttt   <  1:409982/71‑1 (MQ=255)
                                                                   cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGttt   <  1:396966/71‑1 (MQ=255)
                                                                    gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt  >  1:3711623/1‑71 (MQ=255)
                                                                    gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt  <  1:769697/71‑1 (MQ=255)
                                                                    gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt  <  1:1457525/71‑1 (MQ=255)
                                                                    gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt  >  1:3224484/1‑71 (MQ=255)
                                                                    |                                                                      
TTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTT  >  minE/1255581‑1255719

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: