Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 412,564 A→C 75.0% 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,5640AC75.0% 99.0 / 21.0 48intergenic (‑209/+80)leuS/lntleucyl‑tRNA synthetase/apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (9/3);  new base C (20/16);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGA  >  minE/412499‑412633
                                                                 |                                                                     
aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                                   >  1:865419/1‑70 (MQ=255)
aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTAAAGTCAGtt                                                                   >  1:4312707/1‑70 (MQ=255)
  tAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAAc                                                                >  1:2356132/1‑71 (MQ=255)
   aaTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAAc                                                                >  1:1339761/1‑70 (MQ=255)
   aaTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAAc                                                                <  1:5446104/70‑1 (MQ=255)
    aTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     <  1:1552750/64‑1 (MQ=255)
    aTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     <  1:2528471/64‑1 (MQ=255)
      ttACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                                   <  1:3149465/64‑1 (MQ=255)
       tACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:3154624/1‑61 (MQ=255)
       tACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:4392056/1‑61 (MQ=255)
        aCAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAccc                                                           <  1:1582794/70‑1 (MQ=255)
           gagTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGga                                                       <  1:2790775/71‑1 (MQ=255)
            agTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAcc                                                            >  1:3607345/1‑65 (MQ=255)
            agTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGAcc                                                            >  1:500724/1‑65 (MQ=255)
                 gtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGa                                                  <  1:4998281/70‑1 (MQ=255)
                   gggCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAAc                                               >  1:2246144/1‑71 (MQ=255)
                        tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACg                                              <  1:1099357/67‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:930309/1‑42 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:3078737/1‑42 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTa                                                                  <  1:265076/45‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTa                                                                  <  1:4826649/45‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTa                                                                  <  1:2974816/45‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTa                                                                  <  1:3069223/45‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTa                                                                  <  1:2144831/45‑1 (MQ=255)
                          ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATc                                         <  1:2033952/70‑1 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:5129627/1‑40 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     >  1:329951/1‑40 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGAt                                                    >  1:1860628/1‑57 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGAt                                                    >  1:228747/1‑57 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACgat                                            >  1:3233202/1‑63 (MQ=255)
                            agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACgat                                            >  1:2018026/1‑63 (MQ=255)
                              caaagcaaGGAACAAAGAACGTCTATTATTATAGTCAg                                                                     <  1:418194/38‑1 (MQ=255)
                                  gcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTaa                                                                 >  1:59228/1‑38 (MQ=255)
                                  gcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTaa                                                                 >  1:2529939/1‑38 (MQ=255)
                                  gcaaGGAACAAAGAACGTCTATTATTATAGTCAGTTaa                                                                 >  1:4884259/1‑38 (MQ=255)
                                    aaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGga                                                       <  1:3936509/46‑1 (MQ=255)
                                    aaGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGga                                                       <  1:3943331/46‑1 (MQ=255)
                                     aGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGc                             >  1:3087334/1‑71 (MQ=255)
                                      ggAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGcc                            >  1:2898740/1‑71 (MQ=255)
                                           aaaGAACGTCTATTATTATAGTCAGTTAACGACCCggg                                                        <  1:5499309/38‑1 (MQ=255)
                                           aaaGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGaaa                        >  1:2321233/1‑70 (MQ=255)
                                           aaaGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGaaa                        >  1:2247054/1‑70 (MQ=255)
                                               aaCGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATgg                   <  1:3641835/71‑1 (MQ=255)
                                                aCGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGaaa                                                >  1:4270683/1‑41 (MQ=255)
                                                   tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg               >  1:4603936/1‑71 (MQ=255)
                                                            ataGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCgg       >  1:3816750/1‑70 (MQ=255)
                                                                tAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGa  >  1:3182209/1‑71 (MQ=255)
                                                                tAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGa  >  1:286313/1‑71 (MQ=255)
                                                                 |                                                                     
ATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGA  >  minE/412499‑412633

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: