Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 168,776 A→G 16.5% K26K (AAA→AAG rpsB → 30S ribosomal subunit protein S2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE168,7760AG16.5% 199.9 / 29.8 85K26K (AAA→AAGrpsB30S ribosomal subunit protein S2
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (6/8);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTACCGATG  >  minE/168707‑168845
                                                                     |                                                                     
tGTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGaaaa                                                                      >  1:570678/1‑71 (MQ=255)
 gTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAAt                                                                     >  1:2691611/1‑71 (MQ=255)
 gTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAAt                                                                     >  1:3082617/1‑71 (MQ=255)
 gTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAAt                                                                     >  1:1477803/1‑71 (MQ=255)
  tttCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAAt                                                                     <  1:4877715/70‑1 (MQ=255)
         gcgaGACATGCTCAAGGCTGGTGTTAACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGt                                                             >  1:1938177/1‑71 (MQ=255)
          cgcgACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGtt                                                            >  1:2080191/1‑71 (MQ=255)
          cgcgACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGtt                                                            >  1:331490/1‑71 (MQ=255)
            cgACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCa                                                          >  1:2351522/1‑71 (MQ=255)
            cgACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCa                                                          >  1:3721297/1‑71 (MQ=255)
                aTGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCtt                                                      <  1:3924032/71‑1 (MQ=255)
                 tGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTc                                                     <  1:4992802/71‑1 (MQ=255)
                   cTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCgg                                                   <  1:4024852/71‑1 (MQ=255)
                   cTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCgg                                                   <  1:3843249/71‑1 (MQ=255)
                   cTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCgg                                                   <  1:181749/71‑1 (MQ=255)
                        ggCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCg                                                    >  1:4026613/1‑65 (MQ=255)
                              tgtTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAAct                                        <  1:2321346/71‑2 (MQ=255)
                              tgtTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACa                                        <  1:1047860/71‑1 (MQ=255)
                                  cACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGt                                    >  1:2150660/1‑71 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:3635253/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:2762509/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:4640291/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:3984658/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:4006017/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:2827670/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:925167/71‑1 (MQ=255)
                                   aCTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:4969061/71‑1 (MQ=255)
                                       cGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCtt                                                      >  1:1425125/1‑48 (MQ=255)
                                          tCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGtt                                   <  1:3640484/64‑1 (MQ=255)
                                          ccACCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCGGtc                                                 <  1:2411383/49‑2 (MQ=255)
                                          ccACCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCGGtc                                                 <  1:3183768/49‑2 (MQ=255)
                                          ccACCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCGGtc                                                 <  1:3278474/49‑2 (MQ=255)
                                          ccACCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCGGtc                                                 <  1:3691794/49‑2 (MQ=255)
                                          ccACCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCGGtc                                                 <  1:4473186/49‑2 (MQ=255)
                                            aCCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCg                                                    <  1:4117064/45‑1 (MQ=255)
                                            aCCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCg                                                    <  1:345409/45‑1 (MQ=255)
                                            aCCAGACCCGTTTCTGGAACCCAAAGATGGCGCCGTTCATCTTCg                                                    <  1:4524942/45‑1 (MQ=255)
                                            aCCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACatcatc                          <  1:1405846/71‑1 (MQ=255)
                                               aGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAAc                       >  1:4709798/1‑71 (MQ=255)
                                               aGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAAc                       >  1:3284444/1‑71 (MQ=255)
                                                      ttACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatca                            <  1:5164079/59‑1 (MQ=255)
                                                      ttACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTg                   <  1:3305117/68‑1 (MQ=255)
                                                      ttACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTg                   <  1:1815441/68‑1 (MQ=255)
                                                      ttACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTg                   <  1:619016/68‑1 (MQ=255)
                                                      ttACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTg                   <  1:949097/68‑1 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:3345596/3‑68 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:1957858/3‑68 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:3153511/3‑68 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:613621/3‑68 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:3981477/3‑68 (MQ=255)
                                                       ttcTGGAACCCAAAGATGGCGCCGTTCATCTTCGGTCATCGCAACAAGATCCACATCATCAACCTGgaa                 >  1:3975713/3‑68 (MQ=255)
                                                        aCTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAAc                       >  1:4794908/1‑62 (MQ=255)
                                                        aCTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGaa               >  1:4375234/1‑70 (MQ=255)
                                                        aCTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGaa               >  1:4675061/1‑70 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:994371/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:2002149/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:3781457/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:3776523/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:2006176/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:2231271/69‑1 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTg          >  1:2662749/1‑71 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTg          >  1:3133820/1‑71 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTg          >  1:3621608/1‑71 (MQ=255)
                                                            gAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTg          >  1:3371208/1‑71 (MQ=255)
                                                             aaCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGaaa              <  1:3869276/66‑1 (MQ=255)
                                                             aaCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGaaa              <  1:4648313/66‑1 (MQ=255)
                                                             aaCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGt         <  1:3921620/71‑1 (MQ=255)
                                                             aaCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGt         <  1:4215555/71‑1 (MQ=255)
                                                             aaCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGt         <  1:591818/71‑1 (MQ=255)
                                                              aCCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGt         <  1:2549488/70‑1 (MQ=255)
                                                               cccGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTa        >  1:2329402/1‑70 (MQ=255)
                                                               cccGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTa        >  1:1064246/1‑70 (MQ=255)
                                                               cccGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTa        >  1:3796152/1‑70 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:2604352/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:4916655/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:433530/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:5056753/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:3793428/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:1653585/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:1522768/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:1159104/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCAcatc                             >  1:3954565/1‑46 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACatcatc                          <  1:757374/49‑1 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACatcatc                          >  1:945975/1‑49 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTACCGa    >  1:3355882/1‑71 (MQ=255)
                                                                  gAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTACCGa    >  1:1157110/1‑71 (MQ=255)
                                                                    aaaTGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTACCGATg  >  1:4179055/1‑71 (MQ=255)
                                                                     aaTGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAAc            <  1:1896171/60‑1 (MQ=255)
                                                                     |                                                                     
TGTTTCCATGCGCGACATGCTCAAGGCTGGTGTTCACTTCGGTCACCAGACCCGTTACTGGAACCCGAAAATGAAGCCGTTCATCTTCGGTGCGCGTAACAAAGTTCACATCATCAACCTTGAGAAAACTGTACCGATG  >  minE/168707‑168845

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: