Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,847,097 G→C 24.2% G284G (GGG→GGC iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,847,0970GC24.2% 148.9 / 45.7 93G284G (GGG→GGCiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (46/24);  new base C (17/5);  total (64/29)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GACGCAGAAATATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATG  >  minE/1847029‑1847161
                                                                    |                                                                
gACGCAGAAATATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTggg                                                                  <  1:2303369/69‑1 (MQ=255)
gACGCAGAAATATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTggg                                                                  <  1:2326197/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTggg                                                                  >  1:2508298/1‑61 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          <  1:664246/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:1490725/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:2735002/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:1309765/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          <  1:3191119/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:3158885/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:2891536/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:880349/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          <  1:2976031/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          <  1:964497/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          <  1:2937816/69‑1 (MQ=255)
        aaTATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:2625187/1‑69 (MQ=255)
        aaTATTCGACAAAGCGGCCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAgta                                                          >  1:2788615/1‑67 (MQ=255)
        aaTATTAGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGa                                                          >  1:435768/1‑69 (MQ=255)
            ttCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGt                                                      <  1:2787652/69‑1 (MQ=255)
            ttCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGt                                                      <  1:1020807/69‑1 (MQ=255)
                             tGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     <  1:764097/69‑1 (MQ=255)
                             tGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     <  1:59053/69‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:2696747/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:1923686/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:1889212/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:2049563/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:2177948/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:433229/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:390882/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAg                                                            <  1:1804345/45‑1 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:408165/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:2923471/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:2596339/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:2804015/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:546189/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:1577381/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:1801056/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:49146/1‑68 (MQ=255)
                              gCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCaaa                                     >  1:511328/1‑68 (MQ=255)
                                  tACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCGACGCGCAAAAacac                                >  1:1834892/1‑69 (MQ=255)
                                  tACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAacac                                >  1:1334256/1‑69 (MQ=255)
                                  tACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAacac                                >  1:2862906/1‑69 (MQ=255)
                                  tACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAacac                                >  1:2992193/1‑69 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCaa                                                                >  1:165401/1‑35 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:473367/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:2687683/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:2978215/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1013365/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:3142062/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:981266/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1166899/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1220322/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1265141/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1306806/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1541559/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:1859064/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:3173617/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacagc                               >  1:301212/1‑66 (MQ=255)
                                    cTTTCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacag                                >  1:172988/1‑66 (MQ=255)
                                    atttCGGTGGAAGCGACTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAacag                                >  1:2401904/2‑66 (MQ=255)
                                      ttCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGc                            >  1:869376/1‑69 (MQ=255)
                                            ggAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTccc                      >  1:1133200/1‑69 (MQ=255)
                                            ggAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTccc                      >  1:998719/1‑69 (MQ=255)
                                                   aCTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGa               <  1:85623/69‑1 (MQ=255)
                                                   aCTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGa               <  1:2624990/69‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGa               <  1:2385756/68‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGa               <  1:1353011/68‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGa               <  1:3414888/68‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAACAGCCGCATTccc                      <  1:2901969/61‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAACAGCCGCATTccc                      <  1:891852/61‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAACAGCCGCATTccc                      <  1:1135598/61‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAACAGCCGCATTccc                      <  1:1260132/61‑1 (MQ=255)
                                                    cTAACTGGAATCTTGGCAATAAGGATGGTTATCAGCAACGCGCAAAAACAGCCGCATTccc                      <  1:2153827/61‑1 (MQ=255)
                                                      aaCTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAAt            <  1:3424269/69‑1 (MQ=255)
                                                           gAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCt        >  1:1057058/1‑68 (MQ=255)
                                                           gAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCt        >  1:1275830/1‑68 (MQ=255)
                                                           gAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCt        >  1:1480864/1‑68 (MQ=255)
                                                           gAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCt        >  1:546355/1‑68 (MQ=255)
                                                           gAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCt        >  1:1572346/1‑68 (MQ=255)
                                                             aTCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTgg      <  1:2816045/68‑1 (MQ=255)
                                                                ttGGTAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:1901695/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCaaaaa                                   >  1:3406616/1‑36 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAaca                                 >  1:2216808/1‑38 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAaca                                 >  1:2229212/1‑38 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAaca                                 >  1:1551938/1‑38 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAaca                                 >  1:669323/1‑38 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAaca                                 >  1:237949/1‑38 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCTGCGGGAAATAGCTGGc     >  1:961398/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:1931192/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:1019121/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:2579758/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:1104819/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:2465367/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:3060521/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:2441115/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:2320801/1‑66 (MQ=255)
                                                                ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGc     >  1:2917167/1‑66 (MQ=255)
                                                                  gggAATAAGGATGGTTATCAACAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCAt   >  1:1890293/1‑66 (MQ=255)
                                                                    gAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATg  <  1:2399258/65‑1 (MQ=255)
                                                                    gAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATg  <  1:699817/65‑1 (MQ=255)
                                                                    gAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATg  <  1:538321/65‑1 (MQ=255)
                                                                    |                                                                
GACGCAGAAATATTCGACAAAGCGGGCATTGCTGTACTTTCGGTGGAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATG  >  minE/1847029‑1847161

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: