Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,719,183 A→G 41.4% V331V (GTA→GTG argD → bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,719,1830AG41.4% 28.4 / 64.1 87V331V (GTA→GTGargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (29/22);  new base G (9/27);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GGCATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCAGAGCTGAAACCACAGTACAAAGG  >  minE/2719118‑2719245
                                                                 |                                                              
ggCATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTAttt                                                             <  1:1258015/69‑1 (MQ=255)
ggCATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTAttt                                                             <  1:1795277/69‑1 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:1472509/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:585490/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:462291/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:2294225/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:1245543/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:3009700/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:2224604/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:2686671/1‑68 (MQ=255)
  cATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTa                                                            >  1:363336/1‑68 (MQ=255)
      cAGGCGAAACGCCATCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:1447415/69‑1 (MQ=255)
      cAGGCGAAACGCCAGTGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       >  1:3319883/1‑69 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:3423900/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:1871459/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:32611/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:3166246/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:146098/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:101748/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:660942/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:1313199/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:2865005/69‑1 (MQ=255)
      cAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:824282/69‑1 (MQ=255)
        ggCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGata                                                      <  1:241361/68‑1 (MQ=255)
         gCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATAtt                                                    <  1:1610964/69‑1 (MQ=255)
               cGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGt                                              <  1:1321098/69‑1 (MQ=255)
               cGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGt                                              <  1:18344/69‑1 (MQ=255)
                  cAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATAtt                                                    <  1:39965/60‑1 (MQ=255)
                   aGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGat                                                       <  1:1446402/56‑1 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:1811646/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:2994885/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:2530055/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:283299/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:1861572/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:1179653/1‑68 (MQ=255)
                         ttgttgACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:421032/1‑68 (MQ=255)
                         ttgttgACCATCTGAAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTg                                     >  1:3224629/1‑68 (MQ=255)
                           cgagACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  <  1:2726330/66‑1 (MQ=255)
                           cgagACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  <  1:2479535/66‑1 (MQ=255)
                               atcATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:3281118/3‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATGAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:1057078/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:1207066/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:365247/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:1049349/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:301510/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:1337053/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:1550217/1‑69 (MQ=255)
                               aCCATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              >  1:450210/1‑69 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATTGGGCTGTTGATTg                              <  1:512294/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:3102638/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:3264155/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:437554/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2379263/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:674604/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:703700/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2748342/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2142972/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:1022636/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:1118733/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:135702/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:1671233/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:1696448/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:179297/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:267962/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:228335/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2296776/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2303215/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2309643/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2603394/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2379355/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:241884/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:2566897/68‑1 (MQ=255)
                                ccATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:252489/68‑1 (MQ=255)
                                 cATTTGCAGAAGATCGATCAGCAGTACGATGTGTTTAGCGATATTCGCGGTATGGGGCTGTTGATTg                              <  1:1719668/67‑1 (MQ=255)
                                  ttttGCAGAAGATCGATCAGCAGTACGATGTGTTTTGCGATTTTGGCGGTATGGGGCTGTTGATTg                              <  1:2033035/63‑1 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:230219/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:787549/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:1114450/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:493637/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:1783921/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTg                                  >  1:2277780/1‑59 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCag                        >  1:1675475/1‑69 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCag                        >  1:2834391/1‑69 (MQ=255)
                                     tGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCag                        >  1:2311850/1‑69 (MQ=255)
                                                  tcAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCAGAGCTGAAACcac            <  1:2465753/68‑1 (MQ=255)
                                                           cGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCAGAGCTGAAACCACAGTACAAAgg  >  1:157440/1‑69 (MQ=255)
                                                           cGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCAGAGCTGAAACCACAGTACAAAgg  >  1:441150/1‑69 (MQ=255)
                                                                 |                                                              
GGCATTCAGGCGAAACGCCAGCGTTTTGTTGACCATCTGCAGAAGATCGATCAGCAGTACGATGTATTTAGCGATATTCGCGGTATGGGGCTGTTGATTGGCGCAGAGCTGAAACCACAGTACAAAGG  >  minE/2719118‑2719245

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: