Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 412,564 A→C 12.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,5640AC12.1% 121.0 / 17.4 58intergenic (‑209/+80)leuS/lntleucyl‑tRNA synthetase/apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/14);  new base C (3/4);  total (40/18)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

ACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGG  >  minE/412497‑412628
                                                                   |                                                                
aCATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGt                                                               >  1:1070044/1‑71 (MQ=255)
aCATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGt                                                               >  1:2805333/1‑71 (MQ=255)
 cATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGt                                                               >  1:3735519/1‑70 (MQ=255)
  aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                              >  1:3658326/1‑70 (MQ=255)
  aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                              >  1:3663132/1‑70 (MQ=255)
  aTTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                              >  1:2865936/1‑70 (MQ=255)
         tACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGAccc                                                      >  1:2419853/1‑71 (MQ=255)
                  tgtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTacc                                                           >  1:616835/1‑55 (MQ=255)
                   gtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGa                                             <  1:3199598/70‑1 (MQ=255)
                   gtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGa                                             <  1:2780015/70‑1 (MQ=255)
                        cGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACgg                                        >  1:3572496/1‑70 (MQ=255)
                        cGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACgg                                        >  1:2252221/1‑70 (MQ=255)
                          tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTa                                                                  <  1:1065062/42‑1 (MQ=255)
                          tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTa                                                                  <  1:4062031/42‑1 (MQ=255)
                          tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGtt                                                              <  1:39446/46‑1 (MQ=255)
                          tATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGtt                                                              <  1:2104347/46‑1 (MQ=255)
                           aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATc                                    >  1:3166437/1‑71 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGa                                                         >  1:2828395/1‑49 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTaa                                                            >  1:3787570/1‑46 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGaaa                                           >  1:2032578/1‑63 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATcc                                   >  1:2221765/1‑71 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATcc                                   >  1:3932766/1‑71 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATcc                                   >  1:3967515/1‑71 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATcc                                   >  1:344865/1‑71 (MQ=255)
                            ttAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATc                                    >  1:1827329/1‑70 (MQ=255)
                              agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCgg                                 >  1:3526532/1‑71 (MQ=255)
                               gcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGt                                >  1:2418179/1‑71 (MQ=255)
                                 aaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGgag                                                 <  1:1503642/52‑1 (MQ=255)
                                 aaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGt                                >  1:1296239/1‑69 (MQ=255)
                                    gcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCg                           >  1:511175/1‑71 (MQ=255)
                                     caaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGt                          <  1:2430500/71‑1 (MQ=255)
                                        ggAACAAAGAACGTCTATTATTATAGTCAGTTAAc                                                           <  1:2198958/35‑1 (MQ=255)
                                        ggAACAAAGAACGTCTATTATTATAGTCAGTTAAc                                                           <  1:701295/35‑1 (MQ=255)
                                        ggAACAAAGAACGTCTATTATTATAGTAAGTTAACGa                                                         <  1:1023061/37‑1 (MQ=255)
                                         gAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGaaa                                           >  1:3395739/1‑50 (MQ=255)
                                          aaCAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCg                      >  1:734113/1‑70 (MQ=255)
                                           aCAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCg                      >  1:3628298/1‑69 (MQ=255)
                                             aaaGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGaaa                   <  1:316210/70‑1 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGaa                                            >  1:2662559/1‑43 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCa                                      >  1:741354/1‑49 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:709592/1‑71 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:1208454/1‑71 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:267737/1‑71 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:1296862/1‑71 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:1397924/1‑71 (MQ=255)
                                               aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                <  1:1451127/71‑1 (MQ=255)
                                               aGAACGACTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAatat                >  1:7789/1‑71 (MQ=255)
                                                gAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATg               <  1:3283324/71‑1 (MQ=255)
                                                 aaCGTCTATTATTATAGTAAGTTAACGACCCGGGAGATg                                              >  1:2878790/1‑39 (MQ=255)
                                                     tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAat                  >  1:2680980/1‑63 (MQ=255)
                                                     tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg          <  1:2044833/71‑1 (MQ=255)
                                                     tCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATg          <  1:4007374/71‑1 (MQ=255)
                                                       tattattatAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTa        >  1:2300909/1‑71 (MQ=255)
                                                       tattattatAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTa        >  1:1219822/1‑71 (MQ=255)
                                                       tattattatAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTa        >  1:2852720/1‑71 (MQ=255)
                                                           attatAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTcc    >  1:1021251/1‑71 (MQ=255)
                                                               taGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCgg  <  1:1293042/69‑1 (MQ=255)
                                                                  tAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATgg              <  1:2750102/54‑1 (MQ=255)
                                                                  tAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTcc    <  1:1921080/64‑1 (MQ=255)
                                                                  tAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTcc    <  1:699184/64‑1 (MQ=255)
                                                                   |                                                                
ACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGG  >  minE/412497‑412628

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: