Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,576,572 T→C 59.0% E509E (GAA→GAG treF ← cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,576,5720TC59.0% 31.4 / 55.3 78E509E (GAA→GAGtreFcytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (25/7);  new base C (33/13);  total (58/20)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCCATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTTCAT  >  minE/2576506‑2576639
                                                                  |                                                                   
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTAttttt                                                                    >  1:2342235/1‑68 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:656791/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:1147399/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:2306083/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:1306049/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:272276/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:1410695/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:2777622/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:1468403/1‑69 (MQ=255)
cccATCCTGCAACGGATACTCGCCACCGCCGCATTCGCGGGGAACACCATCGGCAATATGGTATTTTTc                                                                   >  1:1777096/1‑69 (MQ=255)
      cTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAg                                                             >  1:1060108/1‑69 (MQ=255)
      cTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAg                                                             >  1:1211696/1‑69 (MQ=255)
      cTGCAACGGATACTCGCCACCGCCCCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAg                                                             >  1:3128782/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:1054682/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:1804110/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:686454/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:1097422/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:3224033/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:2453582/1‑69 (MQ=255)
           aCGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:3068340/1‑69 (MQ=255)
           aCGGATACTCGCAACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTtgt                                                        >  1:2683252/1‑69 (MQ=255)
                     gCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCa                                               >  1:1415716/1‑68 (MQ=255)
                     gCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCa                                               >  1:2322736/1‑68 (MQ=255)
                     gCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCa                                               >  1:2788540/1‑68 (MQ=255)
                     gCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCa                                               >  1:2309282/1‑68 (MQ=255)
                        aCCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATa                                           <  1:1431682/69‑1 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:356047/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:2341256/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:3235741/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:2508428/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:3154741/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:371108/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:2020088/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:253375/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:2538597/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:2588974/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:670759/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:1698174/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:696474/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:792734/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:3120535/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:1363868/1‑60 (MQ=255)
                           gccgccTTCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTc                                                 >  1:3071897/1‑60 (MQ=255)
                           gccgccATCGCGAGGAACACCATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGtt                                                  >  1:734927/1‑59 (MQ=255)
                              gccTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTg                                     <  1:1802375/69‑1 (MQ=255)
                                      ggggAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCAcc                              <  1:587513/68‑1 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGt                                                   >  1:1103302/1‑39 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:2523761/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:2990180/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:3245772/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:399126/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1955167/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1941810/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:63510/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1806647/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1521261/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:846406/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1429699/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1207808/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1177298/1‑67 (MQ=255)
                                              ccATCGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTca                       >  1:1109511/1‑67 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:313124/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:472131/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:3128407/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2473682/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2868864/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2820980/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:1460971/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2638272/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2634383/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2499043/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:1963721/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:1832273/69‑2 (MQ=255)
                                                  cGGCAATATGGTATTTCTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTcagccagt                 <  1:2422379/69‑2 (MQ=255)
                                                          tGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTcgcgcg         >  1:1937978/1‑69 (MQ=255)
                                                            gTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       <  1:1911341/69‑1 (MQ=255)
                                                            gTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       <  1:3162746/69‑1 (MQ=255)
                                                            gTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       <  1:1004114/69‑1 (MQ=255)
                                                             tATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       <  1:1681646/68‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGa        >  1:2046663/1‑63 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:1574200/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:1597667/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:172768/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:65297/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:761139/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:2017373/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:597833/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:502776/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGAt       >  1:2773115/1‑64 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:300150/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:2205937/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:346452/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:2464339/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:1499663/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:947387/69‑1 (MQ=255)
                                                                 tttCGATCAGTTTGTGCTGTTCCAGACAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:238643/69‑1 (MQ=255)
                                                                  ttCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTtcat  <  1:3138149/68‑1 (MQ=255)
                                                                  |                                                                   
CCCATCCTGCAACGGATACTCGCCACCGCCGCCTTCGCGGGGAACACCATCGGCAATATGGTATTTTTCGATCAGTTTGTGCTGTTCCAGATAGAACTGATTCACCGTCTTCAGCCAGCTTCGCGCGATTTCAT  >  minE/2576506‑2576639

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: