Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,156,744 A→G 44.1% L228L (TTG→CTG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,156,7440AG44.1% 9.3 / 45.5 59L228L (TTG→CTG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (12/21);  new base G (12/14);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGC  >  minE/1156688‑1156810
                                                        |                                                                  
aCCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGt                                                        >  1:947890/1‑69 (MQ=255)
aCCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGt                                                        >  1:2422717/1‑69 (MQ=255)
          tcAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTg                                                         >  1:2070120/1‑58 (MQ=255)
          tcAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTa                                              >  1:3044311/1‑69 (MQ=255)
          tcAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTa                                              >  1:1846618/1‑69 (MQ=255)
               ggCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCt                                               <  1:874323/63‑1 (MQ=255)
               ggCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCt                                               <  1:1546206/63‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:1264815/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:3429695/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:29931/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:2976295/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:2816779/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:218615/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:1979074/69‑1 (MQ=255)
                gCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:1141960/69‑1 (MQ=255)
                 cGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:1873004/68‑1 (MQ=255)
                 cGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                        <  1:897004/68‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2101168/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2984110/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:1364929/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:3239420/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:1585686/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:341402/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2182781/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2149691/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2123020/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2056916/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:2053038/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:1341658/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:600942/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                     <  1:696457/69‑1 (MQ=255)
                       ttttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                 >  1:1654251/1‑69 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCaa                                <  1:3255002/69‑1 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                 <  1:2838037/68‑1 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                 <  1:2889671/68‑1 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                 <  1:262556/68‑1 (MQ=255)
                         ttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCaaa                               >  1:1303754/1‑69 (MQ=255)
                                  tGACCGGAGAGCTGGCGATGCAACTGGTAGCCTTTATTGGCGCTACCGGTAAACAGCAAACCGTCgaga                      <  1:1942759/69‑1 (MQ=255)
                                  tGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCgaga                      <  1:3375486/69‑1 (MQ=255)
                                  tGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCgaga                      <  1:3145888/69‑1 (MQ=255)
                                  tGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCgaga                      <  1:2493957/69‑1 (MQ=255)
                                   gACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCTCTACCGGTAAACAGCAAACCGTCgaga                      <  1:780244/68‑1 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:3254600/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:1714821/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:1927370/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:927418/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:3334779/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:2339395/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:3125445/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:1355029/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:2544320/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:1367178/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:2724894/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                       >  1:1557868/1‑64 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc  >  1:1858909/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc  >  1:2003896/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc  >  1:1586423/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc  >  1:335400/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc  >  1:1718775/1‑69 (MQ=255)
                                                        |                                                                  
ACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGC  >  minE/1156688‑1156810

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: