Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,003,911 G→A 28.6% R113C (CGC→TGC)  pobA ← p‑hydroxybenzoate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,003,9110GA28.6% 52.2 / 23.9 35R113C (CGC→TGC) pobAp‑hydroxybenzoate hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base A (5/5);  total (16/19)
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.94e-01

AGCCGGCGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTT  >  ACB143/4003776‑4004045
                                                                                                                                       |                                                                                                                                      
aGCCGGCGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTcgcg                                                                                                                                        <  7:88341/136‑1 (MQ=255)
aGCCGGCGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTcgcg                                                                                                                                        <  6:248740/136‑1 (MQ=255)
      cGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTg                                                                                                                                  <  1:146070/136‑1 (MQ=255)
      cGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTg                                                                                                                                  <  6:255800/136‑1 (MQ=255)
      cGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTg                                                                                                                                  <  3:18888/136‑1 (MQ=255)
                 cAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCAc                                                                                                                       >  5:311631/1‑136 (MQ=255)
                   gTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTg                                                                                                                                  >  1:74200/1‑123 (MQ=255)
                   gTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTg                                                                                                                                  <  2:74200/123‑1 (MQ=255)
                      aaTCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCCCTGGCTTCGCAGGCCTGCATCAGGTCCCGGGTg                                                                                                                  >  7:273073/1‑136 (MQ=255)
                           gCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTc                                                                                                             <  7:294607/136‑1 (MQ=255)
                           gCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTc                                                                                                             <  4:186950/136‑1 (MQ=255)
                            cTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCg                                                                                                            <  8:65836/136‑1 (MQ=255)
                               cACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCCCCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTc                                                                                                         >  7:151137/1‑136 (MQ=255)
                               cACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTc                                                                                                         >  2:206752/1‑136 (MQ=255)
                                 cccGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTg                                                                                                       >  2:29782/1‑136 (MQ=255)
                                  ccGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCATGGGTGACTTCGGTCTgg                                                                                                      <  3:263987/136‑1 (MQ=255)
                                   cGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGTTCACGGGTGACTTCGGTCTGGc                                                                                                     >  2:229358/1‑136 (MQ=255)
                                                                  tctcGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAg                                                                      >  8:101637/1‑136 (MQ=255)
                                                                      gCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCt                                                                  >  7:303898/1‑136 (MQ=255)
                                                                      gCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCt                                                                  >  7:114656/1‑136 (MQ=255)
                                                                        cTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGAGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGCTCTGGACGTAGACAATCACAGTCTTGCCGCCGGTCAGGGCttt                                                                >  6:209479/1‑136 (MQ=255)
                                                                           tCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCttt                                                                >  2:221894/1‑133 (MQ=255)
                                                                                       gCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTg                                                 >  1:222201/1‑136 (MQ=255)
                                                                                       gCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTg                                                 >  8:198685/1‑136 (MQ=255)
                                                                                           aaCGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTggc                                               <  8:220422/134‑1 (MQ=255)
                                                                                              gTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCGTTGAGATCCAGACGCTGGCGgcgc                                          <  8:262019/136‑1 (MQ=255)
                                                                                                 ggcggcTGAATAAATGCTCGGCGCACCACTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGcgccc                                        <  4:91772/135‑1 (MQ=255)
                                                                                                     gcTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTg                                                                                                                                  <  7:295909/41‑1 (MQ=255)
                                                                                                     gcTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTg                                                                                                                                  >  8:295909/1‑41 (MQ=255)
                                                                                                         aaTAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAAcc                                 <  7:172962/134‑1 (MQ=255)
                                                                                                           tAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACcagca                             <  6:72696/136‑1 (MQ=255)
                                                                                                                 aTCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCa                        <  6:311631/135‑1 (MQ=255)
                                                                                                                    ggCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAAccc                    <  4:16603/136‑1 (MQ=255)
                                                                                                                      cgcACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGCCCATCCCCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCCCCCAGCAGGTCAAccccc                  >  7:59832/1‑136 (MQ=255)
                                                                                                                           cacTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTg             <  1:90540/136‑1 (MQ=255)
                                                                                                                             cTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCa           <  1:44628/136‑1 (MQ=255)
                                                                                                                             cTGGCTTCGCAGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCa           <  1:29782/136‑1 (MQ=255)
                                                                                                                                     gcgGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCt   <  7:101637/136‑1 (MQ=255)
                                                                                                                                      cgGGCCTGCATCAGGTCACGGGTGACTTCGGGCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGGGCCCGCGAACCAGCAGTTCAACCCCCTCGTGCACCAGGCGtt  <  6:33997/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
AGCCGGCGATATAGTCGCAGTCAATCCGCTGCACCCGGCCATCCTTTTCGAACGTCAGGTAGGGCTTCTCGCCTTTCAATTCATGCGGCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTT  >  ACB143/4003776‑4004045

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