Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,897,223 A→G 29.4% P315P (CCA→CCG argA → 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,897,2230AG29.4% 92.0 / 43.5 95P315P (CCA→CCGargAfused acetylglutamate kinase homolog (inactive) and amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/41);  new base G (13/15);  total (39/56)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAA  >  minE/1897155‑1897286
                                                                    |                                                               
gAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCa                                                                 >  1:800519/1‑69 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:1479170/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:1461671/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:1401542/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2668550/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:1347279/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:1345687/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2763032/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:862214/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2206488/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2119929/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:3184359/69‑1 (MQ=255)
       ccGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:818673/69‑1 (MQ=255)
        cGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:91300/68‑1 (MQ=255)
        cGAGCAGATTCGTCGCGCAACAATCAACGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2259412/68‑1 (MQ=255)
         gAGCAGATTCGTCGCGCAACAATCATCGATATCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGagc                                                          <  1:2484970/67‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:1873488/69‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:73782/69‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:1245847/69‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:1200757/69‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:351389/69‑1 (MQ=255)
              gATTCGTCGCGCAACAATCAACGATATTGGCGGAATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAgg                                                   <  1:3142079/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAAGGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGt                                                  >  1:2788910/1‑68 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGt                                                  >  1:2629138/1‑68 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 >  1:558509/1‑69 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 >  1:3378884/1‑69 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:824599/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 >  1:3179443/1‑69 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:1767693/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:1190648/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:925718/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:1854208/69‑1 (MQ=255)
                ttCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTa                                                 <  1:936442/69‑1 (MQ=255)
                      gcgcAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTgg                                           <  1:1728433/69‑1 (MQ=255)
                      gcgcAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTgg                                           <  1:167166/69‑1 (MQ=255)
                      gcgcAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTgg                                           <  1:2739634/69‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:1120087/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:2847593/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:1221564/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:1451846/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:1453242/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:1466533/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:160171/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:3040000/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:849257/60‑1 (MQ=255)
                             aaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCt                                             <  1:3076280/60‑1 (MQ=255)
                              aTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGt                                   <  1:2681626/69‑1 (MQ=255)
                              aTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGt                                   <  1:1603264/69‑1 (MQ=255)
                              aTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGt                                   <  1:115095/69‑1 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:3315917/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:2938241/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:3372361/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:275683/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:89444/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtctc                               >  1:2463467/1‑68 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtct                                >  1:3368779/1‑67 (MQ=255)
                                   cGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTtct                                >  1:1029741/1‑67 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGt                                          >  1:1718001/1‑55 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:3376157/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:346360/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:675606/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:679556/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:793765/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:2773093/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:3038475/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:2076207/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:2455724/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:1011712/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:2790564/1‑69 (MQ=255)
                                     atatCGGCGGTATTCTGGAGTTGATTCGCCCGCTGGAGCAGCAAGGTATTCTGGTGCGCCGTTCTcgcg                            >  1:2838834/1‑69 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGgag                 <  1:2858235/69‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGgag                 <  1:673468/69‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGgag                 <  1:2692554/69‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:949178/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:1065018/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:2973233/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:65387/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:1667677/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:2066903/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGga                  <  1:2150018/68‑1 (MQ=255)
                                                aTTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTAGCGAGCAGCTGga                  <  1:553223/68‑1 (MQ=255)
                                                 ttCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGgaga                >  1:1998028/1‑69 (MQ=255)
                                                   cTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGAt               >  1:2835708/1‑68 (MQ=255)
                                                   cTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGAt               >  1:1342446/1‑68 (MQ=255)
                                                   cTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGAt               >  1:490910/1‑68 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTTCGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:3153124/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:790086/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:2205898/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:884348/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:894222/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:1398087/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:1224956/1‑62 (MQ=255)
                                                           gATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg             >  1:418447/1‑62 (MQ=255)
                                                               cGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACaaa  >  1:2977410/1‑69 (MQ=255)
                                                               cGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACaaa  <  1:15792/69‑1 (MQ=255)
                                                               cGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACaaa  <  1:278471/69‑1 (MQ=255)
                                                                    |                                                               
GAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAA  >  minE/1897155‑1897286

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