Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 412,492 C→T 18.7% intergenic (‑137/+152) 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,4920CT18.7% 146.6 / 26.2 64intergenic (‑137/+152)leuS/lntleucyl‑tRNA synthetase/apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/22);  new base T (5/7);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

ATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTA  >  minE/412422‑412559
                                                                      |                                                                   
atagcatagcCCAAACGTCCGCGTCAAAACAGCCTTTCTCGCACTCGACGTTGAAATGATGCCCGGATTAc                                                                     <  1:2388802/71‑1 (MQ=255)
 tagcatagcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACt                                                                    >  1:3453857/1‑71 (MQ=255)
     atagcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATa                                                                <  1:1404499/71‑1 (MQ=255)
        gcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACAt                                                             >  1:898685/1‑71 (MQ=255)
        gcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACAt                                                             >  1:1677529/1‑71 (MQ=255)
          ccAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTc                                                                   >  1:2829458/1‑63 (MQ=255)
          ccAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTc                                                                   >  1:317797/1‑63 (MQ=255)
            aaaCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAAt                                                         >  1:2946677/1‑71 (MQ=255)
            aaaCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAAt                                                         >  1:1917871/1‑71 (MQ=255)
             aaCGTCCGCGTCAAAACAGCCTTTCTCGCACTCGACGTTGAAATGATGGCGGATTACTCa                                                                  >  1:4021438/1‑60 (MQ=255)
             aaCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCa                                                                  >  1:763870/1‑61 (MQ=255)
             aaCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAAtt                                                        <  1:529523/71‑1 (MQ=255)
              aCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTAATACATTAAtt                                                        >  1:2155204/1‑70 (MQ=255)
               cGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTa                                                      <  1:1618528/71‑1 (MQ=255)
                 tCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTAc                                                                     <  1:3053609/54‑1 (MQ=255)
                     cgTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACagag                                                 <  1:136302/70‑1 (MQ=255)
                      gTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGtt                                               <  1:36780/71‑1 (MQ=255)
                       tCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGttt                                              <  1:3674367/71‑1 (MQ=255)
                         aaaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACt                                                                    >  1:711619/1‑47 (MQ=255)
                         aaaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTc                                                                   >  1:2268428/1‑48 (MQ=255)
                         aaaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACag                                                   >  1:1495423/1‑64 (MQ=255)
                         aaaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtg                                             <  1:648889/70‑1 (MQ=255)
                          aaaCAGCCTTTCGCGCACTCGACGTTGATATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtg                                             <  1:1190062/69‑1 (MQ=255)
                          aaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtgtg                                           >  1:1312853/1‑71 (MQ=255)
                           aaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATa                                                                >  1:2850461/1‑49 (MQ=255)
                           aaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATa                                                                >  1:1083890/1‑49 (MQ=255)
                                ccTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTa                                     <  1:3730663/71‑1 (MQ=255)
                                 cTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATTCATTAATTTACAGAGTTTGTGGGCGTAt                                    <  1:2070649/71‑1 (MQ=255)
                                 cTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAAtt                                                        >  1:691051/1‑51 (MQ=255)
                                  tttCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTAt                                    >  1:3837424/1‑70 (MQ=255)
                                  tttCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTAt                                    >  1:3164205/1‑70 (MQ=255)
                                    tCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtgtg                                           >  1:3797935/1‑61 (MQ=255)
                                    tCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTag                                 >  1:402850/1‑71 (MQ=255)
                                        gcACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGttt                                              <  1:11162/54‑1 (MQ=255)
                                            tCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTAt                                    <  1:1573833/60‑1 (MQ=255)
                                             cGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAg                        >  1:4042868/1‑71 (MQ=255)
                                             cGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAg                        >  1:2296633/1‑71 (MQ=255)
                                              gACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtgtg                                           >  1:1033498/1‑51 (MQ=255)
                                              gACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtgtg                                           >  1:1657858/1‑51 (MQ=255)
                                              gACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtgt                                            >  1:2403935/1‑50 (MQ=255)
                                              gACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAg                        >  1:253562/1‑70 (MQ=255)
                                                 gTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGaa                     >  1:722316/1‑70 (MQ=255)
                                                 gTTGAAATGATGCACGGATTACTCATACATTAATTTACAGAGttt                                              >  1:327981/1‑45 (MQ=255)
                                                  ttGAAATGATGCCCGGATTACTCATACATTAATTTACaga                                                  <  1:2524501/40‑1 (MQ=255)
                                                  ttGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGttt                                              <  1:1246278/44‑1 (MQ=255)
                                                    gAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTtg                                             >  1:947702/1‑43 (MQ=255)
                                                    gAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACa                   >  1:4017472/1‑69 (MQ=255)
                                                     aaaTGATGCCCGGATTATTCATACATTAATTTACag                                                   <  1:479715/36‑1 (MQ=255)
                                                     aaaTGATGCCCGGATTATTCATACATTAATTTACag                                                   <  1:2701531/36‑1 (MQ=255)
                                                     aaaTGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACaaa                 <  1:985994/70‑1 (MQ=255)
                                                     aaaTGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACaaa                 <  1:2207209/70‑1 (MQ=255)
                                                     aaaTGATGCCCGGATTACTCATACATTAATTTACAGAGTTtg                                             <  1:3027604/42‑1 (MQ=255)
                                                      aaTGATGCCCGGATTACTCATACATTAATTTACAGAGTTtg                                             <  1:667239/41‑1 (MQ=255)
                                                      aaTGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACaaa                 >  1:1268263/1‑69 (MQ=255)
                                                        tgatgCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTagcaaagc                           <  1:972631/57‑1 (MQ=255)
                                                        tgatgCCCGGATTACTCATACATTAATTTACAGAGTTtgtg                                           >  1:2095050/1‑41 (MQ=255)
                                                         gatgCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACg            >  1:1807049/1‑71 (MQ=255)
                                                          atgCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTagcaaagca                          <  1:1270661/56‑1 (MQ=255)
                                                           tgCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTc          <  1:1253275/71‑1 (MQ=255)
                                                           tgCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTc          <  1:33949/71‑1 (MQ=255)
                                                             cccGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAATAACGTCt         <  1:4008617/70‑1 (MQ=255)
                                                                ggATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCtatta     >  1:1861141/1‑71 (MQ=255)
                                                                ggATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAg                        <  1:3067001/52‑1 (MQ=255)
                                                                  aTTACTCATACATTAATTTACAGAGTTTGTGGGCg                                       <  1:2426891/35‑1 (MQ=255)
                                                                   ttACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCtattatta  <  1:712463/71‑1 (MQ=255)
                                                                     aCTCATACATTAATTTACAGAGTTTGTGGGCGTATTagca                               <  1:3681650/40‑1 (MQ=255)
                                                                      |                                                                   
ATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTA  >  minE/412422‑412559

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: