Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,157 C→A 11.5% intergenic (+129/+187) 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,1570CA11.5% 152.1 / 6.0 52intergenic (+129/+187)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (3/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTA  >  minE/337087‑337225
                                                                      |                                                                    
ccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                                                      >  1:1036875/1‑71 (MQ=39)
 cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                                     >  1:300237/1‑71 (MQ=255)
 cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                                                      <  1:286048/70‑1 (MQ=255)
   aGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCgg                                                                   <  1:2826427/71‑1 (MQ=255)
    gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                                                  <  1:4024984/71‑1 (MQ=255)
     cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                                 <  1:1640671/71‑1 (MQ=255)
      aCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                                                      <  1:3707386/65‑1 (MQ=255)
         aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                             >  1:3480069/1‑71 (MQ=255)
         aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                                                              >  1:494883/1‑70 (MQ=255)
          ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                                                  <  1:2518887/65‑1 (MQ=255)
            gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATa                                                                >  1:2111306/1‑65 (MQ=255)
             ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                                         <  1:2578733/71‑1 (MQ=39)
              cTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                                         >  1:123200/1‑70 (MQ=40)
                    cGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                                         <  1:4044008/64‑1 (MQ=40)
                      aCGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                                <  1:2595084/71‑1 (MQ=39)
                           tACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                                           >  1:3004984/1‑71 (MQ=39)
                           tACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                                           >  1:2480968/1‑71 (MQ=39)
                               cGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCgcatccggcagt                                       <  1:3735837/71‑3 (MQ=11)
                                 tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgccg                                                  <  1:4045133/58‑1 (MQ=11)
                                 tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgccg                                                  <  1:541952/58‑1 (MQ=11)
                                 tCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCg                                                                    <  1:1005731/40‑1 (MQ=18)
                                   ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATaa                                                               >  1:3022205/1‑43 (MQ=14)
                                   ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGa                                                                  >  1:129888/1‑40 (MQ=18)
                                   ttATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAg                                   >  1:2326016/1‑71 (MQ=255)
                                    tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGttc                                                        >  1:2213155/1‑48 (MQ=11)
                                      tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgcc                                                   <  1:1979013/52‑1 (MQ=11)
                                      tCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAc                                <  1:2318689/71‑1 (MQ=255)
                                       cAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGc                                  <  1:2027619/68‑1 (MQ=255)
                                        aGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAcc                               <  1:2605495/70‑1 (MQ=255)
                                         ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa                                 >  1:2805773/1‑67 (MQ=255)
                                         ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa                                 >  1:2810483/1‑67 (MQ=255)
                                         ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAt                             <  1:4016339/71‑1 (MQ=255)
                                          gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAt                             >  1:2053727/1‑70 (MQ=255)
                                           ccTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                                                             >  1:3945881/1‑37 (MQ=14)
                                           ccTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTa                                                       >  1:3800090/1‑43 (MQ=14)
                                             tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                                         <  1:1453513/39‑1 (MQ=40)
                                             tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                                         <  1:1825543/39‑1 (MQ=14)
                                             tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCgca                                                <  1:1788757/48‑1 (MQ=11)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                                       >  1:102917/1‑39 (MQ=40)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatc                                              >  1:1808571/1‑48 (MQ=38)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt                        >  1:2976338/1‑70 (MQ=255)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt                        <  1:1752473/70‑1 (MQ=255)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt                        <  1:1419618/70‑1 (MQ=255)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt                        >  1:21914/1‑70 (MQ=255)
                                               cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTCCCTg                       <  1:3968937/71‑1 (MQ=255)
                                                aaaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGa                      <  1:107642/71‑1 (MQ=255)
                                                 aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAt                             >  1:1870416/1‑63 (MQ=255)
                                                  aCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                                  >  1:2299923/1‑41 (MQ=38)
                                                   cTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGc                   <  1:3794814/71‑1 (MQ=255)
                                                     tgtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGa                 >  1:2492461/1‑71 (MQ=255)
                                                       tgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACg               <  1:4001086/71‑1 (MQ=255)
                                                             aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTAcg         >  1:487563/1‑71 (MQ=40)
                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCt    >  1:1394095/1‑71 (MQ=39)
                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa  >  1:26922/1‑71 (MQ=39)
                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa  >  1:2553543/1‑71 (MQ=39)
                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAg                                   <  1:2621777/37‑1 (MQ=40)
                                                                      |                                                                    
CCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTA  >  minE/337087‑337225

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: