Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,896,276 G→A 54.0% intergenic (‑229/‑3) amiC ← / → argA N‑acetylmuramoyl‑L‑alanine amidase/fused acetylglutamate kinase homolog (inactive) and amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,896,2760GA54.0% 3.3 / 52.3 63intergenic (‑229/‑3)amiC/argAN‑acetylmuramoyl‑L‑alanine amidase/fused acetylglutamate kinase homolog (inactive) and amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (8/21);  new base A (21/13);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATCAA  >  minE/1896211‑1896340
                                                                 |                                                                
acTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:2017986/68‑1 (MQ=255)
acTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:185149/68‑1 (MQ=255)
acTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:901229/68‑1 (MQ=255)
acTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:729464/68‑1 (MQ=255)
acTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:422390/68‑1 (MQ=255)
 cTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGcc                                                                <  1:2951470/67‑1 (MQ=255)
   tGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTg                                                             >  1:2623022/1‑68 (MQ=255)
         cTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAgg                                                      <  1:994668/69‑1 (MQ=255)
         cTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAgg                                                      <  1:1260932/69‑1 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAgg                                                      >  1:2880593/1‑68 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1658093/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:2179720/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1060929/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1989445/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1875264/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1108904/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1818777/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1020622/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:2325587/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1376175/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:814523/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1601468/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:2788222/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1596636/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:2915688/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:467111/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:3150829/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGa                                                     >  1:1382469/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGAAAAGGa                                                     >  1:2223957/1‑69 (MQ=255)
          tGCAAAAAAACAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGAAAAGGa                                                     >  1:2220071/1‑69 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAATGAACGt                                                 <  1:2733689/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:2981242/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:103329/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:606795/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:3050808/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:2986707/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:1460670/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:151664/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:2766286/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:1977891/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:2533620/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:1695043/69‑1 (MQ=255)
              aaaaaaCAGAATAAAAATACAATAATTTCGAATAATCATGCAAAGAGGTGTACCGTGGTAAAGGAACGt                                                 <  1:2372540/69‑1 (MQ=255)
                     aGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCg                                          <  1:1657568/69‑1 (MQ=255)
                     aGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCg                                          <  1:934953/69‑1 (MQ=255)
                     aGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCg                                          <  1:123818/69‑1 (MQ=255)
                        aTAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGt                                       <  1:3268196/69‑1 (MQ=255)
                        aTAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGt                                       <  1:1335489/69‑1 (MQ=255)
                           aaaaTACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTgg                                    >  1:3192545/1‑69 (MQ=255)
                                   cTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGa                            <  1:444332/69‑1 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:1436626/1‑69 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:1413286/1‑69 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:3422293/1‑69 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:676382/1‑69 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:128622/1‑69 (MQ=255)
                                                atCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCgg               >  1:246802/1‑69 (MQ=255)
                                                   aTGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTc            <  1:2950816/69‑1 (MQ=255)
                                                           agGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATc    <  1:2917731/69‑1 (MQ=255)
                                                           agGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATc    <  1:1297137/69‑1 (MQ=255)
                                                           agGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATc    <  1:160062/69‑1 (MQ=255)
                                                           agGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATc    <  1:1230087/69‑1 (MQ=255)
                                                           agGTGTGCCGTGGTAAAAGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATc    <  1:2741917/69‑1 (MQ=255)
                                                             gtgtGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATCaa  <  1:1761059/69‑1 (MQ=255)
                                                               gtgCCGTGGTAAAGGTACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTAGGTTCCCTATATCaa  <  1:2201641/67‑1 (MQ=255)
                                                                 |                                                                
ACTTGCACCCTGCGAAAAAACAGAATAAAAATACACTAATTTCGAATAATCATGCAAAGAGGTGTGCCGTGGTAAAGGAACGTAAAACCGAGTTGGTCGAGGGATTCCGCCATTCGGTTCCCTATATCAA  >  minE/1896211‑1896340

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