Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,128 G→T 29.4% intergenic (+100/+216) ybaQ → / ← copA predicted DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE337,1280GT29.4% 144.2 / 39.8 68intergenic (+100/+216)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base T (4/16);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAT  >  minE/337058‑337187
                                                                      |                                                           
ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAg                                                             <  1:617656/71‑1 (MQ=39)
ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAg                                                             >  1:5361954/1‑71 (MQ=39)
ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAg                                                             <  1:255182/71‑1 (MQ=39)
 gATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAgg                                                            <  1:2850960/71‑1 (MQ=39)
  aTAAGGCGTTCACGCCGCATCCGGCATCCGCGCCCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGc                                                           >  1:4705287/1‑71 (MQ=255)
   tAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGcc                                                          >  1:5426725/1‑71 (MQ=39)
   tAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGcc                                                          >  1:5234925/1‑71 (MQ=39)
      ggCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTa                                                        <  1:579626/70‑1 (MQ=40)
      ggCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTa                                                        <  1:3936434/70‑1 (MQ=40)
        cGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACTCGTCTTATCAGGCCTACa                                                      >  1:222679/1‑70 (MQ=25)
          ttCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAAc                                                   <  1:268372/71‑1 (MQ=39)
           tCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                  >  1:4098414/1‑71 (MQ=39)
           tCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                  >  1:4206352/1‑71 (MQ=39)
            cACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACtt                                                 >  1:3328428/1‑71 (MQ=39)
                ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGc                                             >  1:201713/1‑71 (MQ=39)
                  gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCt                                           >  1:5332384/1‑71 (MQ=39)
                   catccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                                          <  1:5002858/71‑1 (MQ=39)
                     tccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaa                                                     <  1:4128991/58‑1 (MQ=35)
                     tccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGaa                                        <  1:4903546/71‑1 (MQ=39)
                      ccggcatccgAGCTCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACtt                                                 <  1:4656307/61‑1 (MQ=255)
                       cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAcc                                      >  1:1167226/1‑71 (MQ=39)
                       cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAcc                                      >  1:1160387/1‑71 (MQ=39)
                          catccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                   <  1:4643591/71‑1 (MQ=39)
                            tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGcc                                                          <  1:3464144/46‑1 (MQ=38)
                            tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                  <  1:2581297/54‑1 (MQ=38)
                            tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAgg                                 <  1:4124728/71‑1 (MQ=39)
                              cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                               >  1:3645992/1‑71 (MQ=39)
                                 gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                            >  1:2203401/1‑71 (MQ=255)
                                 gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                            >  1:2769352/1‑71 (MQ=255)
                                  cACCGTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                           <  1:841036/71‑1 (MQ=255)
                                  cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                           <  1:4611801/71‑1 (MQ=255)
                                  cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                           <  1:3458981/71‑1 (MQ=255)
                                      aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                           <  1:1673968/67‑1 (MQ=255)
                                      aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                       >  1:3156300/1‑71 (MQ=255)
                                       ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgtg                                              >  1:2478334/1‑47 (MQ=255)
                                       ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                                          >  1:2808927/1‑51 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCt                                           >  1:3006852/1‑48 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                         >  1:1919618/1‑66 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                        >  1:3132075/1‑67 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                        >  1:5474537/1‑67 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCg                     <  1:1366499/70‑1 (MQ=40)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                    >  1:1283872/1‑71 (MQ=39)
                                            tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgc                                               <  1:1898891/41‑2 (MQ=255)
                                            tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                   >  1:2726734/1‑53 (MQ=255)
                                            tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                   >  1:5067364/1‑53 (MQ=255)
                                              aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                           <  1:31433/59‑1 (MQ=255)
                                                 cGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg            <  1:2408332/71‑1 (MQ=39)
                                                    cGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcat         <  1:1080292/71‑1 (MQ=39)
                                                      cTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatcc       >  1:4314238/1‑71 (MQ=39)
                                                      cTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatcc       <  1:277113/71‑1 (MQ=39)
                                                         gcgcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg     <  1:709523/67‑1 (MQ=11)
                                                         aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg     >  1:2097000/1‑70 (MQ=40)
                                                           gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAg                                  <  1:3841883/38‑1 (MQ=11)
                                                             gTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                >  1:2300586/1‑38 (MQ=14)
                                                              tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag  <  1:5343446/68‑2 (MQ=11)
                                                              tCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGttt                  <  1:2583956/52‑1 (MQ=11)
                                                              tCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca          <  1:4141989/60‑1 (MQ=35)
                                                               cTTATCATGCCTACAAACTTGTGCCTGAACGGTAg                                  <  1:2725590/35‑1 (MQ=14)
                                                               cTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCg                     >  1:1829399/1‑48 (MQ=11)
                                                                ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCg                     <  1:3592969/47‑1 (MQ=11)
                                                                ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg     <  1:4532308/63‑1 (MQ=11)
                                                                ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg     <  1:3746486/63‑1 (MQ=11)
                                                                ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg     <  1:712778/63‑1 (MQ=11)
                                                                 tATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                >  1:2986557/1‑34 (MQ=14)
                                                                 tATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                >  1:3834808/1‑34 (MQ=14)
                                                                 tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGAt                           <  1:2194573/40‑1 (MQ=14)
                                                                  aTCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                       <  1:1529795/43‑1 (MQ=14)
                                                                  aTCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                       <  1:3865459/43‑1 (MQ=14)
                                                                   tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                   <  1:2434942/46‑1 (MQ=11)
                                                                   tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                   <  1:2441412/46‑1 (MQ=11)
                                                                   tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                   <  1:3298669/46‑1 (MQ=11)
                                                                   tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgca          <  1:2503124/55‑1 (MQ=11)
                                                                      |                                                           
GGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAT  >  minE/337058‑337187

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