Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,139,593 A→G 62.0% A505A (GCT→GCC infB ← fused protein chain initiation factor 2, IF2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,139,5930AG62.0% 38.8 / 43.8 71A505A (GCT→GCCinfBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base A (19/8);  new base G (42/2);  total (61/10)
Fisher's exact test for biased strand distribution p-value = 5.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTCTTCCGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTGG  >  minE/2139539‑2139659
                                                      |                                                                  
ctctTCTGGCAGGATGCCGTACTGGGATAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:671130/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGTATCGGCTTCTg                                                             >  1:3075920/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:1149434/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:551874/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:414340/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:299941/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2957182/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2567106/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2545741/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2475729/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2388540/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2350847/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2234405/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:2019483/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:1817607/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:1395733/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:1273502/1‑62 (MQ=255)
ctctTCTGGCAGGATGCCCTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:3298029/1‑62 (MQ=255)
ctctTCTGGCAGGATACCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTg                                                             >  1:1579828/1‑62 (MQ=255)
  ctTCCGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGa                                                    >  1:2066641/1‑69 (MQ=255)
  ctTCCGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGa                                                    >  1:2062624/1‑69 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:2476966/69‑1 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:2429523/69‑1 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:2199806/69‑1 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:610696/69‑1 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:2860309/69‑1 (MQ=255)
         cAGGATGCCGTACTGGGAGAGTTCGTTCATAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:3348275/69‑1 (MQ=255)
          aGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGtt                                             <  1:384372/68‑1 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAg                                                                   >  1:432381/1‑41 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTTTCGATCTTGTTCACTGc                                       >  1:2112661/1‑69 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTg                                        >  1:399036/1‑68 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGc                                       >  1:2882583/1‑69 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGc                                       >  1:1064760/1‑69 (MQ=255)
               gCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGc                                       >  1:2928264/1‑69 (MQ=255)
                     cTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAaccac                                 <  1:2781949/69‑1 (MQ=255)
                     cTGGGAGAGTTCGTTCTTAACGCGATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAaccac                                 <  1:1674802/69‑1 (MQ=255)
                                    cTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTgcgc                  >  1:1223523/1‑69 (MQ=255)
                                    cTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTgcg                   >  1:761335/1‑68 (MQ=255)
                                    cTTAACACGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTgcgc                  >  1:1829153/1‑69 (MQ=255)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCacc                                    >  1:1645108/1‑43 (MQ=37)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:3088553/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1198228/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:3289092/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:2858213/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:569458/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:629618/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:619088/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1262587/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:853464/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:2899556/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:883173/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:895180/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1290037/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1290749/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:2498405/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1312317/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1371726/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1447760/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1862369/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1727225/1‑45 (MQ=38)
                                          gcgATCCGGATCGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:978494/1‑45 (MQ=38)
                                          gcgATCCGGATGGCTTCTGGTTTATCGATTTTGTTCACAGCAac                                    >  1:1421889/1‑44 (MQ=25)
                                           cgATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAaccacca                               >  1:2911087/1‑49 (MQ=255)
                                           cgATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTc           >  1:3336187/1‑69 (MQ=255)
                                           cgATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTc           >  1:2818086/1‑69 (MQ=255)
                                           cgATCCGGATCAGCTTCTGGTTTATCGATCTTATTCACTGCAACCACCACCGGTACCTGCGCCGCTTTc           >  1:1175369/1‑69 (MQ=255)
                                               ccGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGt                          <  1:736594/50‑1 (MQ=255)
                                                    tCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTgg  >  1:585032/1‑69 (MQ=255)
                                                    tCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTgg  >  1:2958566/1‑69 (MQ=255)
                                                    tCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTgg  >  1:1662209/1‑69 (MQ=255)
                                                    tCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTg   >  1:2808046/1‑68 (MQ=255)
                                                      |                                                                  
CTCTTCCGGCAGGATGCCGTACTGGGAGAGTTCGTTCTTAACGCGATCCGGATCAGCTTCTGGTTTATCGATCTTGTTCACTGCAACCACCACCGGTACCTGCGCCGCTTTCGCGTGCTGG  >  minE/2139539‑2139659

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: