Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 235,640 C→T 22.2% intergenic (+1/‑257) yaiA → / → aroM hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE235,6400CT22.2% 122.8 / 30.4 63intergenic (+1/‑257)yaiA/aroMhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (10/4);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

AGTGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTC  >  minE/235571‑235709
                                                                     |                                                                     
aGTGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGAcc                                                                      <  1:3744687/71‑1 (MQ=255)
aGTGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGAcc                                                                      >  1:2783783/1‑71 (MQ=255)
aGTAAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGCCCCTGACAAAGAGTGAcc                                                                      >  1:3864124/1‑71 (MQ=255)
  tGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGATCg                                                                     >  1:2252620/1‑70 (MQ=255)
   gAACTTCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGAc                                                                       <  1:3025510/67‑1 (MQ=255)
   gAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCa                                                                   >  1:1912191/1‑71 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGATCGCAt                                                                  <  1:453684/71‑1 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGATCGCAt                                                                  >  1:3645465/1‑71 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGATCGCAt                                                                  <  1:2713091/71‑1 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCAt                                                                  >  1:1218267/1‑71 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCAt                                                                  >  1:2549298/1‑71 (MQ=255)
    aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCATGACAAAGAGTGATCGCAt                                                                  >  1:731925/1‑71 (MQ=255)
      cATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCa                                                                <  1:159270/71‑1 (MQ=255)
      cATCCGACCGCTCAGGAAGCGATGCATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCa                                                                <  1:609080/71‑1 (MQ=255)
       aTCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCAt                                                                  >  1:1442778/1‑68 (MQ=255)
         ccGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGATCGCATCAGa                                                              >  1:2816861/1‑70 (MQ=255)
         ccGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGa                                                              >  1:2539946/1‑70 (MQ=255)
         ccGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGa                                                              >  1:1968006/1‑70 (MQ=255)
          cGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGAc                                                             <  1:939108/70‑1 (MQ=255)
          cGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGAc                                                             <  1:2146438/70‑1 (MQ=255)
            aCCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTg                                                           >  1:3357464/1‑70 (MQ=255)
            aCCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGc                                                          <  1:1737221/71‑1 (MQ=255)
                 tCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCgg                                                      >  1:3417312/1‑70 (MQ=255)
                     gAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAgg                                                  >  1:461194/1‑70 (MQ=255)
                       aGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGAt                                               <  1:1179547/71‑1 (MQ=255)
                         cGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCa                                                                >  1:3470960/1‑52 (MQ=255)
                         cGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGa                                                     <  1:3640923/63‑1 (MQ=255)
                          gatggatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCt                                            <  1:3004496/71‑1 (MQ=255)
                            tggatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCt                                            >  1:1989044/1‑69 (MQ=255)
                             ggatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTg                                                           <  1:131867/53‑1 (MQ=255)
                              gatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAg                                         <  1:439423/70‑1 (MQ=255)
                                 gCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGc                                                          <  1:1015727/50‑1 (MQ=255)
                                    aaaaaaCGCTATGAGGACCCTGACAAAGAGTGACCg                                                                     <  1:1575130/36‑1 (MQ=255)
                                    aaaaaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGa                                                     <  1:2994161/52‑1 (MQ=255)
                                    aaaaaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGa                                                     <  1:2982175/52‑1 (MQ=255)
                                    aaaaaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCAc                                   >  1:246547/1‑70 (MQ=255)
                                     aaaaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACt                                  >  1:3717563/1‑70 (MQ=255)
                                      aaaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTa                                 >  1:2521069/1‑70 (MQ=255)
                                       aaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTAc                                >  1:1850174/1‑70 (MQ=255)
                                       aaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACa                               <  1:2386010/71‑1 (MQ=255)
                                        aaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGaa                                                    >  1:3760614/1‑49 (MQ=255)
                                        aaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAggg                                                 >  1:2325660/1‑52 (MQ=255)
                                           gCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAggg                           <  1:1114413/71‑1 (MQ=255)
                                             tATGAGGACCCTGACAAAGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGAt                         >  1:3589625/1‑71 (MQ=255)
                                              aTGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTAc                                <  1:495762/63‑1 (MQ=255)
                                              aTGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATc                        >  1:660995/1‑71 (MQ=255)
                                               tGAGGACCCTGACAAAGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCt                       >  1:2315477/1‑71 (MQ=255)
                                                 aGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGc                     <  1:2563403/71‑1 (MQ=255)
                                                  ggACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCg                    >  1:3920446/1‑71 (MQ=255)
                                                      ccTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTca                 >  1:1636256/1‑70 (MQ=255)
                                                       cTGACAAAGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGt                                        >  1:3081451/1‑46 (MQ=255)
                                                        tGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTg                                           <  1:2378516/42‑1 (MQ=255)
                                                        tGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAt              >  1:225366/1‑71 (MQ=255)
                                                        tGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAt              >  1:283553/1‑71 (MQ=255)
                                                          aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAttt            >  1:2578247/1‑71 (MQ=255)
                                                          aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAttt            >  1:12004/1‑71 (MQ=255)
                                                          aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAttt            >  1:3892992/1‑71 (MQ=255)
                                                          aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAtt             >  1:413117/1‑70 (MQ=255)
                                                            aaaGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGt                                        >  1:1952776/1‑41 (MQ=255)
                                                            aaaGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGt                                        >  1:1989920/1‑41 (MQ=255)
                                                             aaGAGTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTc         <  1:1704170/71‑1 (MQ=255)
                                                             aaGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTc         <  1:2378578/71‑1 (MQ=255)
                                                              agagTGATCGCATCAGACTGCTCGGAAGGGATTCt                                            <  1:1240868/35‑1 (MQ=255)
                                                              agagTGATCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCa        >  1:15301/1‑71 (MQ=255)
                                                              agagTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGTTCTGCGTCACATTTTTCa        >  1:2883096/1‑71 (MQ=255)
                                                              agagTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCa        >  1:1664643/1‑71 (MQ=255)
                                                                    aCCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTc  <  1:1105004/71‑1 (MQ=255)
                                                                     |                                                                     
AGTGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTC  >  minE/235571‑235709

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: