Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 983,505 T→C 23.6% intergenic (+417/‑151) ydfC → / → dicB hypothetical protein/cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE983,5050TC23.6% 94.7 / 27.1 55intergenic (+417/‑151)ydfC/dicBhypothetical protein/cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base T (24/18);  new base C (7/6);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCT  >  minE/983435‑983568
                                                                      |                                                               
gTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGAtgt                                                                 >  1:3649234/1‑71 (MQ=255)
  gACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtgg                                                               >  1:1174929/1‑71 (MQ=255)
    cGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtggtg                                                             <  1:3495389/71‑1 (MQ=255)
    cGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtggtg                                                             <  1:1767904/71‑1 (MQ=255)
    cGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtggt                                                              >  1:176120/1‑70 (MQ=255)
      tttGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGaa                                                           <  1:3659306/71‑1 (MQ=255)
       ttGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAAt                                                          >  1:4636548/1‑71 (MQ=255)
        tGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAAt                                                          <  1:4825560/70‑1 (MQ=255)
        tGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATg                                                         >  1:1136924/1‑71 (MQ=255)
          gCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGc                                                        <  1:4054426/70‑1 (MQ=255)
                tCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAg                                                 >  1:3054534/1‑71 (MQ=255)
                 cAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAgc                                                >  1:1048103/1‑71 (MQ=255)
                 cAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAgc                                                >  1:3741604/1‑71 (MQ=255)
                    ttttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtgg                                                               <  1:3132195/53‑1 (MQ=255)
                      ttACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGAtgt                                                                 >  1:598682/1‑49 (MQ=255)
                        aCTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcgc                                          <  1:3643948/70‑1 (MQ=255)
                            cGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa                                      >  1:3093459/1‑70 (MQ=255)
                              tGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtggtg                                                             <  1:974370/45‑1 (MQ=255)
                              tGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATgtggtg                                                             <  1:1066898/45‑1 (MQ=255)
                              tGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAg                                   >  1:3683709/1‑71 (MQ=255)
                                aCTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGt                                  >  1:3041768/1‑70 (MQ=255)
                                  tGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCg                                                       <  1:4819895/47‑1 (MQ=255)
                                    ctctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaaaa                             >  1:44310/1‑71 (MQ=255)
                                      ctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAAc                            >  1:3472261/1‑70 (MQ=255)
                                       tGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGaa                                                           >  1:2422739/1‑38 (MQ=255)
                                       tGCCGCCCTTTTTAAAGTGAATTTTGGGATGTGGTGaa                                                           >  1:5619713/1‑38 (MQ=255)
                                         ccgccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTg                                                   <  1:3661696/44‑1 (MQ=255)
                                         ccgccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaa                        <  1:701072/71‑1 (MQ=255)
                                           gccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                      >  1:2859085/1‑71 (MQ=255)
                                           gccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                      >  1:3810154/1‑71 (MQ=255)
                                            cccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa                      >  1:3801619/1‑70 (MQ=255)
                                            cccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaa                         <  1:2318810/67‑1 (MQ=255)
                                            cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg                                           >  1:209816/1‑49 (MQ=38)
                                            cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg                                           >  1:5018031/1‑49 (MQ=38)
                                            cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg                                           >  1:4553311/1‑49 (MQ=38)
                                             ccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAg                                   >  1:2177498/1‑56 (MQ=255)
                                               tttttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg                                        <  1:1431139/49‑1 (MQ=255)
                                               tttttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg                                        <  1:2351362/49‑1 (MQ=255)
                                               tttttAAAGTGAATTTTGTGATGCGGTGAATGCGGCta                                                   <  1:1015920/38‑2 (MQ=38)
                                                 tttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAg                                                 >  1:5355539/1‑38 (MQ=37)
                                                 tttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTggca                                      >  1:2288969/1‑47 (MQ=37)
                                                  ttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaa                       <  1:1153227/63‑1 (MQ=255)
                                                  ttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAgtgt                <  1:3432727/70‑1 (MQ=255)
                                                  ttAAAGAGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAgtgt                >  1:60501/1‑70 (MQ=255)
                                                   tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa                                      >  1:3522728/1‑47 (MQ=255)
                                                   tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa                                      >  1:4790084/1‑47 (MQ=255)
                                                   tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa                                      >  1:2777883/1‑47 (MQ=255)
                                                   tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt               <  1:4752168/70‑1 (MQ=255)
                                                    aaaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc                     >  1:4377844/1‑63 (MQ=255)
                                                    aaaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc                     >  1:5115125/1‑63 (MQ=255)
                                                     aaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTAt             <  1:1607623/70‑1 (MQ=255)
                                                        tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                           <  1:3577234/37‑1 (MQ=255)
                                                        tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                           <  1:4286088/37‑1 (MQ=255)
                                                        tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg                                           <  1:4858309/37‑1 (MQ=255)
                                                         gAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa                             >  1:4340373/1‑50 (MQ=255)
                                                         gAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa                             >  1:1109611/1‑50 (MQ=255)
                                                          aaTTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTg        >  1:1791464/1‑70 (MQ=255)
                                                            ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa                               <  1:1661474/45‑1 (MQ=25)
                                                            ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa                               <  1:5535740/45‑1 (MQ=25)
                                                            ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa                               <  1:3726905/45‑1 (MQ=25)
                                                              ttGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGAtt    >  1:5099483/1‑70 (MQ=255)
                                                                gtgATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATtct  <  1:3279245/70‑1 (MQ=255)
                                                                  gATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa                               <  1:4210935/39‑1 (MQ=25)
                                                                  gATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa                               <  1:902769/39‑1 (MQ=25)
                                                                      tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt               <  1:1989506/51‑1 (MQ=255)
                                                                      tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt               <  1:4610845/51‑1 (MQ=255)
                                                                      |                                                               
GTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCT  >  minE/983435‑983568

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: