Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 412,564 A→C 41.1% intergenic (‑209/+80) leuS ← / ← lnt leucyl‑tRNA synthetase/apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE412,5640AC41.1% ‑0.5 / 49.5 39intergenic (‑209/+80)leuS/lntleucyl‑tRNA synthetase/apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base C (7/9);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCAC  >  minE/412496‑412631
                                                                    |                                                                   
tACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                   <  1:2121686/71‑1 (MQ=255)
tACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                   <  1:59915/71‑1 (MQ=255)
  cATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                                 <  1:1297888/71‑1 (MQ=255)
  cATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                                 <  1:1648033/71‑1 (MQ=255)
   aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTa                                                                <  1:1559964/71‑1 (MQ=255)
      aaTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAAc                                                              >  1:2042025/1‑70 (MQ=255)
        tttACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAc                                                           <  1:1926830/71‑1 (MQ=255)
          tACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAcc                                                          <  1:1116043/70‑1 (MQ=255)
          tACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGAccc                                                         >  1:3683701/1‑71 (MQ=255)
            cAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCg                                                        <  1:1129852/70‑1 (MQ=255)
             agagTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                                 >  1:1049356/1‑60 (MQ=255)
             agagTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAc                                                           <  1:1500717/66‑1 (MQ=255)
             agagTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCgg                                                       >  1:3031692/1‑70 (MQ=255)
               agTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCgg                                                       >  1:2257325/1‑68 (MQ=255)
                gTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCgg                                                       >  1:516471/1‑67 (MQ=255)
                   tgtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATg                                                 >  1:682087/1‑70 (MQ=255)
                         cGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACga                                           >  1:2365952/1‑69 (MQ=255)
                         cGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCgg                                                       <  1:923621/58‑1 (MQ=255)
                          gTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACgat                                          >  1:854334/1‑68 (MQ=255)
                           tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTaa                                                               <  1:2551560/48‑1 (MQ=255)
                           tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACgatg                                         >  1:309852/1‑67 (MQ=255)
                            aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTaa                                                               >  1:1498896/1‑47 (MQ=255)
                            aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTaa                                                               >  1:3300776/1‑47 (MQ=255)
                             ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                                 <  1:1682568/44‑1 (MQ=255)
                               agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTaa                                                               >  1:2053636/1‑44 (MQ=255)
                                gcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                                 <  1:1210406/41‑1 (MQ=255)
                                  aaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGt                                   >  1:2434753/1‑69 (MQ=255)
                                    agcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGc                               <  1:3169370/71‑1 (MQ=255)
                                          gAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGa                                 <  1:2689536/63‑1 (MQ=255)
                                             cAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGCGAGATGAAACGGCATCCGGTGAGCGTGCCGaaa                      <  1:3426331/71‑1 (MQ=255)
                                                      tCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAAc                                             <  1:1500354/39‑1 (MQ=255)
                                                      tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAAc                                             <  1:1721425/39‑1 (MQ=255)
                                                      tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg             <  1:393686/71‑1 (MQ=255)
                                                       cTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg             >  1:2281166/1‑70 (MQ=255)
                                                       cTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg             >  1:3618396/1‑70 (MQ=255)
                                                               ataGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCgg     <  1:2297796/70‑1 (MQ=255)
                                                               ataGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCgg     <  1:30054/70‑1 (MQ=255)
                                                                taGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGc    <  1:820353/70‑1 (MQ=255)
                                                                  gTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCAc  >  1:1089075/1‑70 (MQ=255)
                                                                    |                                                                   
TACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCAC  >  minE/412496‑412631

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: