Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,253,716 C→T 34.1% A241A (GCG→GCA aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,253,7160CT34.1% 49.2 / 38.4 44A241A (GCG→GCAaspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (9/20);  new base T (8/7);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCT  >  minE/1253651‑1253782
                                                                 |                                                                  
tCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGt                                                               >  1:1033670/1‑71 (MQ=255)
 ccACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTgcggtgcggt                                                               >  1:957286/1‑70 (MQ=255)
 ccACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTgcggtgcggt                                                               >  1:1777133/1‑70 (MQ=255)
     acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGa                                                          <  1:1690905/71‑1 (MQ=255)
     acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGa                                                          <  1:1769098/71‑1 (MQ=255)
     acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGa                                                          <  1:1223927/71‑1 (MQ=255)
        cTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAaag                                                       <  1:1435026/71‑1 (MQ=255)
          tCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGaa                                                         >  1:629705/1‑67 (MQ=255)
          tCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaag                                                       >  1:1091582/1‑69 (MQ=255)
             cTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGtt                                                  <  1:1985570/71‑1 (MQ=255)
              ttCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaagaa                                                     >  1:1566383/1‑67 (MQ=255)
                      acaGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1141743/71‑1 (MQ=255)
                       caGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGa                                        <  1:2095687/71‑1 (MQ=255)
                        aGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGAt                                       >  1:2147974/1‑71 (MQ=255)
                          aTGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaagaag                                                    <  1:1319850/56‑1 (MQ=255)
                               gCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTa                                <  1:1089275/71‑1 (MQ=255)
                               gCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTa                                <  1:1826546/71‑1 (MQ=255)
                                 aCCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaaa                              <  1:1116951/71‑1 (MQ=255)
                                 aCCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaaa                              >  1:2084883/1‑71 (MQ=255)
                                     gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1200770/71‑1 (MQ=255)
                                     gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1460070/71‑1 (MQ=255)
                                     gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:937098/71‑1 (MQ=255)
                                     gcgcTTCCATCACTTCACGAACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAAttta                          <  1:2094978/71‑3 (MQ=255)
                                       gcTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAgg                        >  1:188148/1‑71 (MQ=255)
                                           ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:663105/1‑71 (MQ=255)
                                           ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:350532/1‑71 (MQ=255)
                                           ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:1195897/1‑71 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1808050/45‑1 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1139482/45‑1 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1986370/45‑1 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:2050477/1‑70 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                <  1:505805/70‑1 (MQ=255)
                                                aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                <  1:715687/70‑1 (MQ=255)
                                                  ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:1160249/39‑1 (MQ=255)
                                                  ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:2180078/39‑1 (MQ=255)
                                                  ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:633780/39‑1 (MQ=255)
                                                   tCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAg             >  1:845593/1‑70 (MQ=255)
                                                   tCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAg             >  1:848079/1‑70 (MQ=255)
                                                    cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:48957/1‑66 (MQ=255)
                                                    cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCa           >  1:1699898/1‑71 (MQ=255)
                                                    cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCa           >  1:938408/1‑71 (MQ=255)
                                                    aaCGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:1887692/2‑66 (MQ=255)
                                                       gCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGc        <  1:2011798/71‑1 (MQ=255)
                                                       gCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGc        <  1:107449/71‑1 (MQ=255)
                                                          cTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1041886/50‑1 (MQ=255)
                                                           ttGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaa                               <  1:1439085/44‑1 (MQ=255)
                                                              cggcgCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCt  <  1:2179843/70‑1 (MQ=255)
                                                                 |                                                                  
TCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCT  >  minE/1253651‑1253782

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: