Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,044,955 A→G 13.4% intergenic (‑42/+87) tyrS ← / ← pdxH tyrosyl‑tRNA synthetase/pyridoxine 5'‑phosphate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,044,9550AG13.4% 175.7 / 17.8 67intergenic (‑42/+87)tyrS/pdxHtyrosyl‑tRNA synthetase/pyridoxine 5'‑phosphate oxidase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (4/5);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGC  >  minE/1044887‑1045024
                                                                    |                                                                     
ttGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                                                                      >  1:2146391/1‑70 (MQ=255)
ttGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                                                                      >  1:109237/1‑70 (MQ=255)
ttGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                                                                      >  1:2248546/1‑70 (MQ=255)
ttGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                                                                      >  1:1576237/1‑70 (MQ=255)
ttGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                                                                      >  1:3564189/1‑70 (MQ=255)
  gTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGAc                                                                   >  1:3159018/1‑71 (MQ=255)
  gTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGAc                                                                   >  1:636725/1‑71 (MQ=255)
      aaTCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                                                                <  1:2383362/70‑1 (MQ=255)
      aaTCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                                                                <  1:2205716/70‑1 (MQ=255)
      aaTCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                                                                <  1:461035/70‑1 (MQ=255)
        tCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                             <  1:3543091/71‑1 (MQ=255)
        tCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                             >  1:927420/1‑71 (MQ=255)
        tCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                             <  1:1983021/71‑1 (MQ=255)
         cAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc                                                            >  1:3123536/1‑71 (MQ=255)
          aaGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc                                                            >  1:887728/1‑70 (MQ=255)
            gTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAg                                                         <  1:195349/71‑1 (MQ=255)
            gTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCCTTGCCGGTacac                                                                      <  1:2643041/58‑1 (MQ=255)
                 tGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACGCGACTTTTCGCCAGATGc                                                     >  1:542440/1‑70 (MQ=255)
                  gCTTGCCATCAAAATCTCCTTGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATg                                                      <  1:1134741/68‑1 (MQ=255)
                  gCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGa                                                                    >  1:1858713/1‑54 (MQ=255)
                  gCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGa                                                                    >  1:1854171/1‑54 (MQ=255)
                      gCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc                                                            >  1:3217698/1‑58 (MQ=255)
                        cATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACGCGACTTTTCGCCAGATGCGAAAgaga                                             >  1:1135932/1‑71 (MQ=255)
                        cATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaaa                                                 <  1:3509878/67‑1 (MQ=255)
                         aTCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTc                                                              <  1:1013372/53‑1 (MQ=255)
                          tCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACa                                           <  1:1030295/71‑1 (MQ=255)
                           cAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc                                                            >  1:1627586/1‑53 (MQ=255)
                             aaaTCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                             >  1:714954/1‑50 (MQ=255)
                             aaaTCTCCATGTATATAAAGACTGCACCTTTGGCGGTACACGACTTTTCGCCAGATGCGAAAGAGAc                                            >  1:488988/1‑67 (MQ=255)
                              aaTCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAg                                        >  1:425697/1‑70 (MQ=255)
                              aaTCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAg                                        >  1:395220/1‑70 (MQ=255)
                                   ccATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                                                                <  1:899482/41‑1 (MQ=255)
                                    cATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaaa                                  <  1:3263740/70‑1 (MQ=255)
                                       gTATATAACGACTGCACCTTTGCCGGTACGCGACTTTTCGCCAGATg                                                      <  1:778165/47‑1 (MQ=255)
                                        tatataACGACTGCACCTTTGCCGGTACGCGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGc                              <  1:2463067/70‑1 (MQ=255)
                                        tatataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTg                               >  1:333542/1‑69 (MQ=255)
                                          tataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGc                              >  1:2990815/1‑68 (MQ=255)
                                          tataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGc                              >  1:3004517/1‑68 (MQ=255)
                                          tataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGc                              >  1:781958/1‑68 (MQ=255)
                                          tataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCa                            <  1:578645/70‑1 (MQ=255)
                                           ataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGa                                                   <  1:1984664/46‑1 (MQ=255)
                                           ataACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATa                                         >  1:777057/1‑56 (MQ=255)
                                             aaCGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGt                                <  1:3044900/63‑1 (MQ=255)
                                               cGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGAc                                            <  1:30832/49‑1 (MQ=255)
                                               cGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGAc                                            <  1:2321680/49‑1 (MQ=255)
                                                   tGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaaa                                                 >  1:1551072/1‑40 (MQ=255)
                                                       ccTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaaa                                  <  1:466845/51‑1 (MQ=255)
                                                       ccTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaaa                                  <  1:401103/51‑1 (MQ=255)
                                                        cTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaaa                                  <  1:1558348/50‑1 (MQ=255)
                                                        cTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAg                     <  1:2242428/63‑1 (MQ=255)
                                                        cTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTAc               >  1:2154898/1‑69 (MQ=255)
                                                         tttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGa                          <  1:904143/57‑1 (MQ=255)
                                                         tttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGa                          <  1:2269614/57‑1 (MQ=255)
                                                          ttGCCGGTACGCGACTTTTCGCCAGATGCGAAAGAGAc                                            >  1:2809931/1‑38 (MQ=255)
                                                          ttGCCGGTACGCGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAg                  >  1:858025/1‑64 (MQ=255)
                                                          ttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATa                                    >  1:629416/1‑46 (MQ=255)
                                                          ttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTAcc              >  1:2904423/1‑68 (MQ=255)
                                                          ttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGc           >  1:2574006/1‑71 (MQ=255)
                                                          ttGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGc           >  1:2569628/1‑71 (MQ=255)
                                                           tGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGa                                       >  1:8356/1‑42 (MQ=255)
                                                           tGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGa                                       >  1:3029059/1‑42 (MQ=255)
                                                            gCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGt                                >  1:1392438/1‑48 (MQ=255)
                                                             ccGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAt        <  1:3045747/71‑1 (MQ=255)
                                                              cGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAtt       <  1:482717/71‑1 (MQ=255)
                                                              cGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAtt       <  1:998587/71‑1 (MQ=255)
                                                              cGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAt        >  1:3195147/1‑70 (MQ=255)
                                                               ggTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAtt       >  1:25456/1‑70 (MQ=255)
                                                               ggTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAtt       <  1:1297389/70‑1 (MQ=255)
                                                               ggTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTa      <  1:3632005/71‑1 (MQ=255)
                                                               ggTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTa      <  1:289119/71‑1 (MQ=255)
                                                                gTACGCGACTTTTCGCCAGATGCGAAAGAGACATAGAATa                                    <  1:2635200/40‑1 (MQ=255)
                                                                gTACGCGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCa             <  1:735818/63‑1 (MQ=255)
                                                                gTACGCGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCa             <  1:1539058/63‑1 (MQ=255)
                                                                  acacGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGa         >  1:473501/1‑65 (MQ=255)
                                                                   cacGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAg   <  1:2662712/70‑1 (MQ=255)
                                                                   cacGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGc  <  1:2314705/71‑1 (MQ=255)
                                                                    |                                                                     
TTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGC  >  minE/1044887‑1045024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: