Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 275,944 C→A 34.6% V77V (GTG→GTT cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE275,9440CA34.6% 134.0 / 88.1 124V77V (GTG→GTTcyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/40);  new base A (21/22);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 125 total aligned reads displayed.

AAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAA  >  minE/275882‑276009
                                                              |                                                                 
aaGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGc                                                             >  1:723014/1‑69 (MQ=255)
aaGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGc                                                             >  1:815682/1‑69 (MQ=255)
aaGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGc                                                             >  1:2453875/1‑69 (MQ=255)
       cAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAt                                                      <  1:1230236/69‑1 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCTAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttct                                                   >  1:457939/1‑68 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:3260077/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:2076752/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:2227579/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:1710575/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:1606177/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:10604/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:1140092/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:426841/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:2995920/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttctt                                                  >  1:2932444/1‑69 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttct                                                   >  1:1551563/1‑68 (MQ=255)
           aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGAttct                                                   >  1:2382119/1‑68 (MQ=255)
                          gTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCtt                                   <  1:1330045/69‑1 (MQ=255)
                            acacCATCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:437058/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGAttcttc                                                 >  1:2553422/1‑53 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:352597/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:3294605/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:2712374/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:674453/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:3094386/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:2170427/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:968930/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:2033165/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:965975/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:719298/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:1363680/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:1425761/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:148465/1‑69 (MQ=255)
                            acacCAGCGAGACAGCAGGAGAGATCAGGCCTTTAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTc                                 >  1:86345/1‑69 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCtt                                   >  1:1832878/1‑62 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGt                              >  1:2382650/1‑67 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:2989590/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:1369030/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:2377120/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:2444952/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:1487076/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTc                             >  1:150835/1‑68 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCg                            <  1:3292048/69‑1 (MQ=255)
                                 aGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCAATCTTTCTGTc                             >  1:2026037/1‑68 (MQ=255)
                                  gCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGa                           <  1:1560890/69‑1 (MQ=255)
                                  gCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGa                           <  1:1985835/69‑1 (MQ=255)
                                  gCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGa                           <  1:1760781/69‑1 (MQ=255)
                                      gaCAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATc                       <  1:2477280/69‑1 (MQ=255)
                                      gaCAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATc                       <  1:2397993/69‑1 (MQ=255)
                                      gaCAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATc                       <  1:902147/69‑1 (MQ=255)
                                       aCAGCAGCAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATc                       <  1:3130947/68‑1 (MQ=255)
                                          gcagGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTg                   >  1:3291737/1‑69 (MQ=255)
                                          gcagGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTg                   >  1:368330/1‑69 (MQ=255)
                                                   aTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           >  1:327293/1‑68 (MQ=255)
                                                   aTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAg          >  1:3423218/1‑69 (MQ=255)
                                                   aTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAg          >  1:1224923/1‑69 (MQ=255)
                                                   aTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAg          >  1:2543191/1‑69 (MQ=255)
                                                   aTCAGGCCTTTCACCAGCACCCGATTATTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAg          >  1:2547395/1‑69 (MQ=255)
                                                     cAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAgtt        >  1:3297887/1‑69 (MQ=255)
                                                         ccTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTa    <  1:926368/69‑1 (MQ=255)
                                                         ccTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTa    <  1:1660363/69‑1 (MQ=255)
                                                         ccTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTa    <  1:2613164/69‑1 (MQ=255)
                                                          cTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTa    <  1:692278/68‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:85853/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:665789/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:484322/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:2613318/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:465849/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:323430/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:1079268/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:1131402/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:1639901/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:1730767/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:1756577/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:2008850/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:2553867/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:262968/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:3274877/69‑1 (MQ=255)
                                                           tttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:2843204/69‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:362078/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:3005519/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:957470/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:1311687/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:1421666/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:2976927/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:857579/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:2816539/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:1825350/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:674864/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:2199246/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:222913/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:542280/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:24807/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:2600148/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:408730/60‑1 (MQ=255)
                                                           tttAACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTa           <  1:427114/60‑1 (MQ=255)
                                                            ttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:183716/68‑1 (MQ=255)
                                                            ttCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTaaa  <  1:613302/68‑1 (MQ=255)
                                                             tCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCg                            <  1:2803604/41‑1 (MQ=255)
                                                              cACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATc                       <  1:2728203/45‑1 (MQ=255)
                                                              |                                                                 
AAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAA  >  minE/275882‑276009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: