Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,573,339 T→C 28.8% intergenic (+29/‑66) ysaB → / → glyQ hypothetical protein/glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,573,3390TC28.8% 91.2 / 47.5 66intergenic (+29/‑66)ysaB/glyQhypothetical protein/glycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (8/11);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAAT  >  minE/2573271‑2573404
                                                                    |                                                                 
tGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtata                                                                  <  1:1281872/70‑1 (MQ=255)
tGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtata                                                                  <  1:586279/70‑1 (MQ=255)
  aCGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCt                                                               <  1:3199686/71‑1 (MQ=255)
     cgGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtat                                                                   <  1:2002985/64‑1 (MQ=255)
       gATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtcca                                                           <  1:3121106/70‑1 (MQ=255)
         tGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtccatc                                                         >  1:3821690/1‑70 (MQ=255)
          ggCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtccatc                                                         >  1:3856533/1‑69 (MQ=255)
           gCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                                                      <  1:1282650/71‑1 (MQ=255)
            ccATTTTTTGCATCTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:3545743/71‑1 (MQ=255)
             cATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCg                                                    >  1:635819/1‑71 (MQ=255)
              aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtccat                                                          >  1:3222381/1‑64 (MQ=255)
              aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCgg                                                   <  1:161784/71‑1 (MQ=255)
               ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                                                  <  1:2590114/71‑1 (MQ=255)
               ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCggg                                                  <  1:23949/71‑1 (MQ=255)
               ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCggg                                                  <  1:2362479/71‑1 (MQ=255)
               ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCgg                                                   >  1:1774731/1‑70 (MQ=255)
                 ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:2320038/1‑71 (MQ=255)
                 ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:1418642/1‑71 (MQ=255)
                 ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCggg                                                  >  1:3625832/1‑69 (MQ=255)
                  tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                <  1:599776/70‑1 (MQ=255)
                   ttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtat                                               <  1:1988829/70‑1 (MQ=255)
                       aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                <  1:1627846/65‑1 (MQ=255)
                       aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtc                                                             <  1:2236741/52‑1 (MQ=255)
                       aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtc                                                             <  1:2168231/52‑1 (MQ=255)
                       aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtc                                                             <  1:492052/52‑1 (MQ=255)
                       aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGa                                          <  1:32536/71‑1 (MQ=255)
                        tCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCt                                                               <  1:1540638/49‑1 (MQ=255)
                         cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATc                                        >  1:2935690/1‑71 (MQ=255)
                         cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATc                                        >  1:3881590/1‑71 (MQ=255)
                          tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCtt                                                              <  1:1466346/48‑1 (MQ=255)
                           ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       >  1:1991702/1‑70 (MQ=255)
                             ccATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTATTGTATATCTTCCATCCAGCGGGTATACTGATCcct                                     <  1:2096044/70‑1 (MQ=255)
                             ccATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCcctt                                    >  1:375782/1‑71 (MQ=255)
                              cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtccatcc                                                        <  1:3199737/50‑1 (MQ=255)
                              cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGTGTATACTGATCcctt                                    <  1:3880849/70‑1 (MQ=255)
                              cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttc                                   >  1:1691747/1‑71 (MQ=255)
                               aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:2984662/1‑59 (MQ=255)
                                tGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCccttcct                                 <  1:568359/71‑1 (MQ=255)
                                tGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCccttcct                                 >  1:43537/1‑71 (MQ=255)
                                 gCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttcctt                                >  1:3261944/1‑71 (MQ=255)
                                 gCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttcctt                                >  1:1423844/1‑71 (MQ=255)
                                   gTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              <  1:2204223/71‑1 (MQ=255)
                                    ttAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                                                  <  1:2783029/50‑1 (MQ=255)
                                    ttAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              >  1:2930578/1‑70 (MQ=255)
                                    ttAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              <  1:3241368/70‑1 (MQ=255)
                                    ttAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              <  1:2981567/70‑1 (MQ=255)
                                     tAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       >  1:349569/1‑60 (MQ=255)
                                     tAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       >  1:745736/1‑60 (MQ=255)
                                     tAAGCCATGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATc                                        >  1:6196/1‑59 (MQ=255)
                                       aGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          <  1:3666062/71‑1 (MQ=255)
                                        gCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATcc                         <  1:1479558/71‑1 (MQ=255)
                                          ccTGCTTTTTCCCGTTTCGTACTGTATATCTtcca                                                           <  1:1497950/35‑1 (MQ=255)
                                          ccTGCTTTTTCCCGTTTCGTACTGTATATCTtcca                                                           <  1:2617077/35‑1 (MQ=255)
                                            tGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCacac                     >  1:1670708/1‑71 (MQ=255)
                                             gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACg                    >  1:111509/1‑71 (MQ=255)
                                                  ttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATc                <  1:1653431/70‑1 (MQ=255)
                                                   tCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          >  1:3804492/1‑59 (MQ=255)
                                                   tCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          >  1:2162060/1‑59 (MQ=255)
                                                       gTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCAc          >  1:3254568/1‑71 (MQ=255)
                                                        tttCGTACTGTATATCTTCCATCCAGCGGGTATACTGAt                                         >  1:2344292/1‑39 (MQ=255)
                                                        tttCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCAc          >  1:257095/1‑70 (MQ=255)
                                                          tCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGa        <  1:655384/70‑1 (MQ=255)
                                                            gTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaat     <  1:948701/71‑1 (MQ=255)
                                                             tACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaat     >  1:3096497/1‑70 (MQ=255)
                                                               cTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTa                  >  1:2059966/1‑55 (MQ=255)
                                                               cTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTa                  >  1:2920224/1‑55 (MQ=255)
                                                               cTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaataa   >  1:293267/1‑70 (MQ=255)
                                                               cTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaataat  <  1:1937055/71‑1 (MQ=255)
                                                               cTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaataa   >  1:796905/1‑70 (MQ=255)
                                                                    |                                                                 
TGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAAT  >  minE/2573271‑2573404

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