Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,410,191 C→T 28.5% N215N (AAC→AAT hemX → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,410,1910CT28.5% 171.8 / 68.8 157N215N (AAC→AAThemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base C (60/52);  new base T (28/17);  total (88/69)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 157 total aligned reads displayed.

ATCGCAGGTGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGT  >  minE/2410128‑2410257
                                                               |                                                                  
aTCGCAGGTGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGc                                                                >  1:1193764/1‑68 (MQ=255)
aTCGCAGGTGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCg                                                               >  1:2182056/1‑69 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTg                                                                 <  1:3040144/59‑1 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTg                                                                 <  1:281714/59‑1 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGcc                                                        >  1:1295752/1‑68 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGcc                                                        >  1:675956/1‑68 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGcc                                                        >  1:3183015/1‑68 (MQ=255)
        tGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGc                                                         >  1:381067/1‑67 (MQ=255)
         ggATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCg                                                       >  1:118599/1‑68 (MQ=255)
         ggATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCg                                                       >  1:1713332/1‑68 (MQ=255)
         ggATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCg                                                       >  1:1626989/1‑68 (MQ=255)
         ggATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCg                                                       >  1:3091461/1‑68 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2916649/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:894336/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:82340/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:1844698/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2775035/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:1865819/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2738410/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:27257/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2134001/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2654140/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:731716/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:1772222/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:1755863/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:2999226/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:306550/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:3211409/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:693699/1‑69 (MQ=255)
                      atcatcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGAt                                         >  1:820247/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAAt                                                  <  1:941210/57‑1 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAAt                                                  <  1:3238459/57‑1 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAAt                                                  <  1:620058/57‑1 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAAt                                                  <  1:2231291/57‑1 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATgg                                       >  1:801574/1‑68 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATgg                                       >  1:3417301/1‑68 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:2848342/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:921856/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:2851725/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:2772626/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:1416006/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:1573273/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:1787204/1‑69 (MQ=255)
                         atcCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGt                                      >  1:1279237/1‑69 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3302170/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3371494/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3314986/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:1495491/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3095704/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3166505/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:217680/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:465273/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:2078954/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:2140691/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:1009433/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:713013/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:1181279/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:2434910/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAATCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:876074/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTg                                                                 <  1:1562078/37‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:470480/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:2968877/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:445823/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:311384/69‑1 (MQ=255)
                              tAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGcc                                 <  1:3157587/69‑1 (MQ=255)
                                        cAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTcg                       >  1:2452708/1‑68 (MQ=255)
                                         aGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAg                      >  1:1031344/1‑69 (MQ=255)
                                         aGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAg                      >  1:3356596/1‑69 (MQ=255)
                                         aGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAg                      >  1:3181231/1‑69 (MQ=255)
                                         aGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAg                      >  1:2285046/1‑69 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:3025707/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:167462/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:338801/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:1116834/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:1403577/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:1143776/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:2465707/69‑1 (MQ=255)
                                           cTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGAc                    <  1:707160/69‑1 (MQ=255)
                                            tttCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACg                   >  1:2168483/1‑69 (MQ=255)
                                            tttCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACg                   >  1:423270/1‑69 (MQ=255)
                                                  aTCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAag             <  1:219040/69‑1 (MQ=255)
                                                  aTCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAag             <  1:530699/69‑1 (MQ=255)
                                                    cAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAagag           <  1:549883/69‑1 (MQ=255)
                                                    cAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAagag           <  1:2239732/69‑1 (MQ=255)
                                                    cAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAagag           <  1:758610/69‑1 (MQ=255)
                                                    cAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAagag           <  1:1848345/69‑1 (MQ=255)
                                                    cAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAagag           <  1:1801928/69‑1 (MQ=255)
                                                            tAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTcc     <  1:3282882/67‑1 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGATGAGCTTTCAAGt  >  1:1375590/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:1361419/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:2642152/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:153048/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:1874545/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:2921577/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:2976881/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:906052/1‑69 (MQ=255)
                                                             aaCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGt  >  1:3082727/1‑69 (MQ=255)
                                                              aCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTg                <  1:2635902/54‑1 (MQ=255)
                                                              aCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTg                <  1:241222/54‑1 (MQ=255)
                                                              aCCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTg                <  1:157685/54‑1 (MQ=255)
                                                               |                                                                  
ATCGCAGGTGGATTATGACGGCATCATCCTTAAGCTTAATCAGCTTTCAAATCAGGTAGATAACCTGCGTCTGGCCGATAATGACAGCGATGGTTCGCCGATGGATTCAGACGGTGAAGAGCTTTCCAGT  >  minE/2410128‑2410257

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: