Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,283 A→G 28.7% N66D (AAT→GAT)  ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,2830AG28.7% 85.0 / 38.4 87N66D (AAT→GAT) ycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base G (20/5);  total (56/31)
Fisher's exact test for biased strand distribution p-value = 8.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTGGGCG  >  minE/626221‑626345
                                                              |                                                              
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTAt                                                           >  1:1746876/1‑68 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:516147/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2577171/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:1731285/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2997019/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:484176/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:1523200/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:594451/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:3017598/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:1393987/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2307125/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:1819781/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2237045/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:131070/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:3362860/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:3420740/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2738872/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:1093865/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          >  1:2234190/1‑69 (MQ=255)
aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTCTCAGCCGTCTGGGAGATTATc                                                          >  1:2102096/1‑69 (MQ=255)
 aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGTGAGATTATc                                                          <  1:3004942/68‑1 (MQ=255)
 aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          <  1:5181/68‑1 (MQ=255)
 aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          <  1:533575/68‑1 (MQ=255)
 aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          <  1:2498651/68‑1 (MQ=255)
 aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                                                          <  1:885545/68‑1 (MQ=255)
   tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                                                       >  1:1427952/1‑69 (MQ=255)
   tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                                                       >  1:867053/1‑69 (MQ=255)
   tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                                                       >  1:3360344/1‑69 (MQ=255)
              aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg                                            >  1:2216587/1‑69 (MQ=255)
              aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg                                            >  1:3272872/1‑69 (MQ=255)
              aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg                                            >  1:757102/1‑69 (MQ=255)
              aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg                                            >  1:349397/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:2730976/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:2808472/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:3293259/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:3324505/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:1066277/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:1306717/1‑69 (MQ=255)
                  ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg                                        >  1:1653515/1‑69 (MQ=255)
                    gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa                                       <  1:2149632/68‑1 (MQ=255)
                    gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa                                       <  1:2574915/68‑1 (MQ=255)
                    gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa                                       <  1:1999733/68‑1 (MQ=255)
                       gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc                                    <  1:2382793/68‑1 (MQ=255)
                       gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc                                    <  1:1357729/68‑1 (MQ=255)
                       gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc                                    <  1:3325009/68‑1 (MQ=255)
                       gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc                                    <  1:3352669/68‑1 (MQ=255)
                          cGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTctct                                 <  1:1269484/68‑1 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGt                          >  1:1410443/1‑68 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:3215624/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:3068812/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:1334499/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:1731875/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:1770403/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:2634035/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:2120053/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:412140/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTg                         >  1:461938/1‑69 (MQ=255)
                                 cATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGAGTg                         >  1:1914437/1‑69 (MQ=255)
                                   tatCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTgg                        >  1:1988987/1‑68 (MQ=255)
                                   tatCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTgg                        >  1:3243071/1‑68 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:1276587/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:1435371/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:3215951/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:1544706/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:1634396/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:269730/1‑69 (MQ=255)
                                     tCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGa                     >  1:236781/1‑69 (MQ=255)
                                           gCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAAtt                  <  1:3135630/66‑1 (MQ=255)
                                           gCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAAtt                  <  1:2618148/66‑1 (MQ=255)
                                           gCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAAtt                  <  1:2320527/66‑1 (MQ=255)
                                           gCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAAtt                  <  1:2152586/66‑1 (MQ=255)
                                              tttCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATgg             <  1:2193294/68‑1 (MQ=255)
                                              tttCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATgg             <  1:3340818/68‑1 (MQ=255)
                                              tttCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATgg             <  1:2035335/68‑1 (MQ=255)
                                                 cAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAAt          <  1:3118812/68‑1 (MQ=255)
                                                 cAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAAt          <  1:402954/68‑1 (MQ=255)
                                                   gtcgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATAc        <  1:956212/68‑1 (MQ=255)
                                                   gtcgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATAc        <  1:1876373/68‑1 (MQ=255)
                                                   gtcgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATAc        <  1:1376260/68‑1 (MQ=255)
                                                   gtcgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATAc        <  1:1348143/68‑1 (MQ=255)
                                                   gtcgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATAc        <  1:1066632/68‑1 (MQ=255)
                                                     cgtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAAtt                  <  1:529690/56‑1 (MQ=255)
                                                      gtgtGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTggg    >  1:1521844/4‑69 (MQ=255)
                                                      gtcTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTggg    >  1:2626610/1‑69 (MQ=255)
                                                        cTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTGGgcg  <  1:300914/69‑1 (MQ=255)
                                                        cTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTGGgcg  <  1:2775429/69‑1 (MQ=255)
                                                        cTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTGGgcg  <  1:1193479/69‑1 (MQ=255)
                                                              |                                                              
AAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCTGGCGCGTGGCGAATTAATGGAATACTGGGCG  >  minE/626221‑626345

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: