Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,156,754 T→G 40.6% T224T (ACA→ACC astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,156,7540TG40.6% 22.6 / 62.4 64T224T (ACA→ACCastDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (16/22);  new base G (12/14);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTC  >  minE/1156688‑1156822
                                                                  |                                                                    
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:1141960/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:1264815/69‑1 (MQ=255)
                 cGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                                    <  1:897004/68‑1 (MQ=255)
                 cGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTa                                                    <  1:1873004/68‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:2101168/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:341402/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:600942/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:696457/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:1364929/69‑1 (MQ=255)
                   agaATTTTCTCCGGCTGACCGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAAc                                                 <  1:2984110/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:1585686/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:3239420/69‑1 (MQ=255)
                       ttttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                             >  1:1654251/1‑69 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCaa                                            <  1:3255002/69‑1 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                             <  1:2889671/68‑1 (MQ=255)
                        tttCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCa                                             <  1:2838037/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:2724894/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:1367178/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:2544320/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:2339395/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:1557868/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:1355029/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:3125445/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:1927370/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:927418/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:3254600/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:1714821/1‑64 (MQ=255)
                                      cGGAAAGCTGGCGATGCAGCTGGTAGCCGGTATTGGCGCTGCCGGTAAACAGCAAGCCATCgag                                   >  1:3334779/1‑64 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc              >  1:1718775/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc              >  1:335400/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc              >  1:2003896/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc              >  1:1586423/1‑69 (MQ=255)
                                                      cAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAgcgc              >  1:1858909/1‑69 (MQ=255)
                                                              aGCCTGTATTGGCGCTACCGGTAAACAGCAAGCCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCg      >  1:1087134/1‑69 (MQ=255)
                                                              aGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCg      >  1:6242/1‑69 (MQ=255)
                                                                 cTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGttt   <  1:2576649/69‑1 (MQ=255)
                                                                 cTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGttt   >  1:152440/1‑69 (MQ=255)
                                                                 cTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGttt   >  1:2422279/1‑69 (MQ=255)
                                                                  tGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGttt   >  1:2708573/1‑68 (MQ=255)
                                                                  tGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTc  >  1:186711/1‑69 (MQ=255)
                                                                  |                                                                    
ACCGCCCATCTCAAGGGCGAGAATTTTCTCCGGCTGACCGGAGAGCTGGCGATGCAACTGGTAGCCTGTATTGGCGCTACCGGTAAACAGCAAACCGTCGAGATCCTCCAGCGCACTCAGCGCCTGACCCGTTTC  >  minE/1156688‑1156822

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: