Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 507,425 C→A 15.1% A71D (GCT→GAT)  ybhH → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE507,4250CA15.1% 388.8 / 31.7 125A71D (GCT→GAT) ybhHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (52/54);  new base A (11/8);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 127 total aligned reads displayed.

ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGC  >  minE/507356‑507495
                                                                     |                                                                      
aCGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:2023637/71‑1 (MQ=255)
 cGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:2533271/70‑1 (MQ=255)
  ggTATTGGCGGCGGTAATCCACTGCCCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGa                                                                     >  1:442034/1‑71 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTg                                                                      <  1:670289/66‑1 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTg                                                                      <  1:429403/66‑1 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTg                                                                      <  1:1556989/66‑1 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTg                                                                      <  1:1525544/66‑1 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATTTCGa                                                               >  1:871863/1‑71 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                <  1:2562027/70‑1 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                <  1:2630038/70‑1 (MQ=255)
            gcggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTAt                                                           <  1:2131832/71‑1 (MQ=255)
             cggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATc                                                          <  1:3935792/71‑1 (MQ=255)
             cggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATc                                                          <  1:3217989/71‑1 (MQ=255)
             cggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATc                                                          <  1:3737020/71‑1 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTATCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgcg                                                  >  1:331150/1‑71 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                                  >  1:1574295/1‑71 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                                  >  1:3805851/1‑71 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                                  >  1:1238991/1‑71 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgcg                                                  >  1:916036/1‑71 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgcg                                                  >  1:1260924/1‑71 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                               <  1:162413/56‑1 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                               <  1:2964827/56‑1 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                               <  1:3166086/56‑1 (MQ=255)
                        tGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTg                                                        >  1:890070/1‑62 (MQ=255)
                        tGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTg                                                        >  1:3465042/1‑62 (MQ=255)
                               tAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:1521541/40‑1 (MQ=255)
                               tAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:2139029/40‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                       <  1:3571978/70‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGtt                                      >  1:3366486/1‑71 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGtt                                      >  1:767836/1‑71 (MQ=255)
                                   gTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:3753835/36‑1 (MQ=255)
                                   gTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                       <  1:4029550/36‑1 (MQ=255)
                                    tCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                   >  1:128410/1‑71 (MQ=255)
                                    tCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTAACTGTTTGCGCAGGTAATCGTTCAt                                   >  1:1675538/1‑71 (MQ=255)
                                     cGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                   >  1:2375792/1‑70 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                               <  1:3348520/41‑1 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:2243138/71‑1 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:4113740/71‑1 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:2156502/71‑1 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:17674/71‑1 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:493326/71‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTa                                            <  1:2174606/57‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTa                                            <  1:1427106/57‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               <  1:4113236/70‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                               >  1:202426/1‑70 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:707194/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:1988523/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:29839/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:701981/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:1244530/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:3916810/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:4095671/40‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                            <  1:672147/40‑1 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAAt                                          >  1:3793048/1‑57 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAAt                                          >  1:4082154/1‑57 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAAt                                          >  1:1857778/1‑57 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGtt                                                      >  1:3312995/1‑44 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGtt                                                      >  1:3171849/1‑44 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGtt                                                      >  1:624931/1‑44 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                       <  1:2280849/71‑1 (MQ=255)
                                                  cGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGa                     <  1:2343483/71‑1 (MQ=255)
                                                  cGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGa                     <  1:2434325/71‑1 (MQ=255)
                                                    ttCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                       <  1:3286748/51‑1 (MQ=255)
                                                    ttCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                       <  1:3525883/51‑1 (MQ=255)
                                                    ttCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                       <  1:853930/51‑1 (MQ=255)
                                                    ttCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                       <  1:1185196/51‑1 (MQ=255)
                                                    ttCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                   <  1:564898/55‑1 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                 >  1:1284586/1‑55 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                 >  1:2942900/1‑55 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                 >  1:366821/1‑55 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                 >  1:520648/1‑55 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:2667772/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:1359107/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:3784684/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:3120224/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:1228897/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:2722436/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:1103808/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:2352911/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                  >  1:832473/1‑70 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:1055382/1‑71 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:3761747/1‑71 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:1007480/1‑71 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:1215797/1‑71 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:1643015/1‑71 (MQ=255)
                                                       cAGCGATCAGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAAt                                          <  1:1629887/45‑1 (MQ=37)
                                                       cAGCGATCAGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAAt                                          <  1:1805776/45‑1 (MQ=37)
                                                        aGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                     <  1:2228346/49‑1 (MQ=255)
                                                        aGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                     <  1:2079674/49‑1 (MQ=255)
                                                        aGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                     <  1:3806404/49‑1 (MQ=255)
                                                          cGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAcc                 >  1:502257/1‑67 (MQ=255)
                                                          cGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAcaa                >  1:1581502/1‑66 (MQ=255)
                                                          cGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                >  1:3631187/1‑68 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAAc                            >  1:2500787/1‑55 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACg                           >  1:3557071/1‑56 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACg                           >  1:1298333/1‑56 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACg                           >  1:370173/1‑56 (MQ=255)
                                                                   tGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTaa         >  1:2281652/1‑66 (MQ=255)
                                                                   tGCTGAAGTCGATTATCTGTTTGCGGAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTaa         >  1:942281/1‑66 (MQ=255)
                                                                   tGCAGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTaa         >  1:2584446/1‑66 (MQ=255)
                                                                     |                                                                      
ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGC  >  minE/507356‑507495

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: