Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,624 C→T 22.9% intergenic (+94/+59) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,6240CT22.9% 178.2 / 48.3 96intergenic (+94/+59)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/27);  new base T (10/12);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTG  >  minE/5554‑5691
                                                                      |                                                                   
cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATc                                                                     >  1:4198018/1‑71 (MQ=255)
  gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                   >  1:1033395/1‑71 (MQ=255)
   cATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                   >  1:3429588/1‑70 (MQ=255)
     ttCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt                                                                    <  1:4489936/67‑1 (MQ=255)
      tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCAt                                                                >  1:1566045/1‑70 (MQ=255)
        aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                                                             >  1:626255/1‑71 (MQ=255)
               gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                                                      >  1:2140600/1‑71 (MQ=255)
                aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                     >  1:2157516/1‑71 (MQ=255)
                 tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                     >  1:3454739/1‑70 (MQ=255)
                   cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                  >  1:4258947/1‑71 (MQ=255)
                   cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                  >  1:4075143/1‑71 (MQ=255)
                   cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                  <  1:5461251/71‑1 (MQ=255)
                    gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                  >  1:3075199/1‑70 (MQ=255)
                      cGCTGGCGCGTCTTATCAGTCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                               <  1:2522751/71‑1 (MQ=255)
                         tGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                   <  1:4933170/48‑1 (MQ=255)
                           gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                                                                 >  1:4682794/1‑48 (MQ=255)
                           gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                           <  1:2819361/70‑1 (MQ=255)
                              cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                       <  1:2253435/71‑1 (MQ=255)
                                tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                                              <  1:1069708/62‑1 (MQ=255)
                                  ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                         >  1:467451/1‑49 (MQ=255)
                                   tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGt                                                        >  1:5188918/1‑49 (MQ=255)
                                   tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                  >  1:3893917/1‑71 (MQ=255)
                                       aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc                               <  1:2775898/70‑1 (MQ=255)
                                       aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc                               <  1:2229800/70‑1 (MQ=255)
                                       aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCa                              >  1:4362773/1‑71 (MQ=255)
                                         gCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                     <  1:3846967/62‑1 (MQ=255)
                                           cTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc                           <  1:3393367/70‑1 (MQ=255)
                                            tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                                              <  1:5540013/50‑1 (MQ=255)
                                            tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                         >  1:5297533‑M1/1‑71 (MQ=255)
                                            tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                         >  1:3953771‑M1/1‑71 (MQ=255)
                                             aCGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                         >  1:3297847‑M1/1‑70 (MQ=255)
                                              cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                       >  1:3059320‑M1/1‑71 (MQ=255)
                                              cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                       >  1:153691‑M1/1‑71 (MQ=255)
                                              cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                       >  1:2783007‑M1/1‑71 (MQ=255)
                                               gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                         >  1:1939379/1‑36 (MQ=255)
                                               gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                         >  1:2783382/1‑36 (MQ=255)
                                               gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                         >  1:3526172/1‑36 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                         <  1:3792198/35‑1 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                       >  1:2543164/1‑37 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     >  1:3747856/1‑71 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     >  1:2146278/1‑71 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     <  1:1288116/71‑1 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     <  1:4113015/71‑1 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     <  1:3994012/71‑1 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     >  1:2323434/1‑71 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     >  1:1654681/1‑71 (MQ=255)
                                                ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     <  1:1825340/71‑1 (MQ=255)
                                                 tAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                     >  1:1791558/1‑70 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                       >  1:154/1‑35 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                       >  1:2031531/1‑35 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAgg                                                 >  1:344469/1‑41 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAgg                                                 >  1:3636197/1‑41 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg                    >  1:5380938/1‑70 (MQ=255)
                                                  aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGa                   >  1:5387744/1‑71 (MQ=255)
                                                   aTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGC                        <  1:1566591‑M1/65‑1 (MQ=255)
                                                    ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTc                                      >  1:4311763/1‑50 (MQ=255)
                                                    ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgcc                                 <  1:79045/55‑1 (MQ=255)
                                                    ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgcc                                 <  1:2981545/55‑1 (MQ=255)
                                                    ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                            >  1:427967/1‑60 (MQ=255)
                                                      cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                             >  1:3338352/1‑41 (MQ=255)
                                                      cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                             >  1:5620756/1‑41 (MQ=255)
                                                      cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                             >  1:691059/1‑41 (MQ=255)
                                                       tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataa                <  1:162167‑M1/69‑7 (MQ=255)
                                                       tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataa                <  1:5640649‑M1/69‑7 (MQ=255)
                                                       tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaat                  <  1:1857813‑M1/67‑5 (MQ=255)
                                                        gCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                               <  1:388539/37‑1 (MQ=255)
                                                         cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcg            <  1:3860884‑M1/71‑11 (MQ=255)
                                                         cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcg            >  1:5612008‑M1/1‑61 (MQ=255)
                                                         cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaa               <  1:2088042‑M1/68‑8 (MQ=255)
                                                         cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaa               <  1:5534453‑M1/68‑8 (MQ=255)
                                                          aaTATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgc           <  1:4347374‑M1/71‑12 (MQ=255)
                                                           atatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgca          <  1:4381488‑M1/71‑13 (MQ=255)
                                                           atatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaat                  <  1:3884668‑M1/63‑5 (MQ=255)
                                                           atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                       <  1:1848515/42‑1 (MQ=255)
                                                           atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                       <  1:3812902/42‑1 (MQ=255)
                                                            tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATga                    >  1:2398733‑M1/1‑58 (MQ=255)
                                                            tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATga                    >  1:5604673‑M1/1‑58 (MQ=255)
                                                            tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa              <  1:110977/66‑1 (MQ=255)
                                                            tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa              <  1:4012164/66‑1 (MQ=255)
                                                            tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACt          >  1:284536/1‑70 (MQ=255)
                                                             atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                       <  1:1638927/40‑1 (MQ=255)
                                                             atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                       <  1:2689296/40‑1 (MQ=255)
                                                             atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                     >  1:1252321/1‑42 (MQ=255)
                                                             atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                     >  1:3930461/1‑42 (MQ=255)
                                                              tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactc       <  1:2571599‑M1/71‑16 (MQ=255)
                                                              tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact        >  1:1579168‑M1/1‑56 (MQ=255)
                                                              tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact        >  1:1537235‑M1/1‑56 (MQ=255)
                                                              tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact        >  1:637586‑M1/1‑56 (MQ=255)
                                                              tatTGAATCTGCATGCTTTTGTAGGAAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTaa       >  1:5696364/1‑71 (MQ=255)
                                                                ttGAATTTGCGTGCTTTTGTAGGCAGGTTAAGGCGTTCACGCCGCATCCGGCATgaataaagcgca          >  1:579422‑M1/1‑54 (MQ=255)
                                                                ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAAc      >  1:1645072/1‑70 (MQ=255)
                                                                 tGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaa     >  1:526587‑M1/1‑53 (MQ=255)
                                                                 tGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaa     >  1:1531146‑M1/1‑53 (MQ=255)
                                                                   aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaa  >  1:4185550‑M1/1‑51 (MQ=255)
                                                                   aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaata                 <  1:4928258‑M1/56‑6 (MQ=255)
                                                                   aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaata                 <  1:1233799‑M1/56‑6 (MQ=255)
                                                                   aaTCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctg  >  1:3503473/1‑71 (MQ=255)
                                                                      |                                                                   
CGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTG  >  minE/5554‑5691

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: