Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 955,841 T→G 26.6% L84R (CTG→CGG)  phnD1 → putative ABC transporter phosphite bindingprotein PhnD1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122955,8410TG26.6% 41.5 / 2.6 37L84R (CTG→CGG) phnD1putative ABC transporter phosphite bindingprotein PhnD1
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base G (3/8);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAAGTTCAAGCCTTTGGGTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAAGAGCCTGAGCGATC  >  ACB122/955712‑955972
                                                                                                                                 |                                                                                                                                   
aGCGCAAGTTCAAGCCTTTGGGTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTggggg                                                                                                                                <  6:143645/135‑1 (MQ=255)
       gTTCAAGCCTTTGGGTGAGTACATGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGCCCCCCACCGACAGGCCGGGCGGGGGGGCGTCCcg                                                                                                                        >  8:140841/1‑136 (MQ=255)
              ccgttGGGTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACgttcgttc                                                                                                                 <  4:292749/133‑1 (MQ=255)
               cTTTGGGTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGAACATGCCTGGCGGGGGGGCGTGAcgttcgtt                                                                                                                  >  3:208545/1‑134 (MQ=255)
                   gggTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGt                                                                                                             <  3:223044/135‑1 (MQ=255)
                    ggTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGt                                                                                                             <  2:76414/134‑1 (MQ=255)
                     ttGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAATTTCGTACCCGTGGAGGACTACCCGGCGGTGGTCGAGTAGCTGGCCGCAGACCGCCTAGCTCTGGCCTGTCTGGGGGGGGTCACATTCCTTATGGTg                                                                                                            >  5:21223/2‑134 (MQ=255)
                     gTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGCCCGCCTCGCTCTGGCCTGGCTGGGGGGGTTCACCTTCGTTCCGGGGc                                                                                                           >  1:221485/1‑135 (MQ=255)
                           aCCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCCGCCCGCCGCGCTCTGGGCTGGCGGGGGGGCGTGACGGTCGTTAAGGTGCAccg                                                                                                       >  1:190576/1‑132 (MQ=255)
                           aCCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCGTCACGTTCGTTCAGGTGCACCt                                                                                                       >  5:89172/1‑133 (MQ=255)
                           aCCTGGCCAAGCAACTGGGAATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCCCTGGCCCCAGCCCCCCGCGCTCTGGCCTGGCTCGGGGGCGTCCAGGTCGTTAACGTggcactg                                                                                                      >  3:69023/1‑128 (MQ=255)
                            ccTGGCCAAGCAACTGGGCATGGCGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGGTCGTTCAGGTGCACCTGAAg                                                                                                   >  8:77094/1‑136 (MQ=255)
                                 ccACGCAACTGGGCAGGGAGGTGAAGGTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGa                                                                                                 <  5:69218/133‑1 (MQ=255)
                                       aaCTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCCCGTTCGTTCAGGTGCACCTGAAGGACCCGACcgc                                                                                        >  7:258347/1‑136 (MQ=255)
                                       aaCTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACcg                                                                                         >  5:154904/1‑135 (MQ=255)
                                       aaCTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGGCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACcgc                                                                                        <  8:50566/136‑1 (MQ=255)
                                        ccTGGGCAGGGAGGTGAAGTCCGACCCCGTGGCTGACTACCCGGCGGTGGTCGAGGCGCGGGCGGCAGACCGCGTCGATCTGCCCTGGCGGGGGGGCTCCAGGTTCGTTCAGGTGCACGTGAAGGACCCGACcgc                                                                                        <  6:154071/134‑1 (MQ=255)
                                          ggggCATGGCGGTGACGGTCGTCCCCGGGGCGGACTCCCCGGCGGGGGGCGAGTCGCGGGCCGCAGACCGCCTCGAGCGGGCCTGGCGGGGGGGCTGCACGTTCGTTCAGGTGCACCTGAAGGACCCGACcgcca                                                                                      <  1:335228/134‑1 (MQ=255)
                                            ggactGGAGGTGAAGTTCGAACCCGGGGCGGACTACCCGGCGGTGGTCGAGGCGCTGGCCGCAGACCGCCTCGATCGGGCCTGGCGGGGGGGCTTCACGTTCGTTCCGGTGCCCCTGAAGGACCCGACcgccacg                                                                                    <  2:333437/131‑1 (MQ=255)
                                            ggCATGGAGGGGAAGTTCGTACCGGGGGCGGACTACCGGGCGGTGGGCGCGGCGCTGGCCGCAGACCGCCTCGATCTGGCCGGGGTGGGGGGGTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACcgccacg                                                                                    <  5:105804/135‑1 (MQ=255)
                                                     gTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCGGGGGGGCTTCACGTTCGTTCCGGTCCACCTGAAGGACCCGACCGCCACGCCACTGGTAc                                                                          <  7:34743/136‑1 (MQ=255)
                                                         aGTGCGTACCCGGGGCGGACTACCCGGCGGTGGGCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGGGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCg                                                                      <  8:393949/136‑1 (MQ=255)
                                                            ccgcccccGGGGCGCACTACCCGGCGGGGGCCGAGGCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCGGGGGGGCTTCACGTGCGTTCAGGTCCACCTCAAGGCCCCGACCGCCACGCCACTGGTACAGCGTg                                                                    <  1:365352/130‑1 (MQ=255)
                                                                aCCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCtt                                                                                                                            >  5:153092/1‑75 (MQ=255)
                                                                aCCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCtt                                                                                                                            <  6:153092/75‑1 (MQ=255)
                                                                   ggTGGCGGACTACCCGGGGGTGGGCGAGGCGCTGGCGGCAGACCGCGGGGAGCTGGCCTGGCGGGGGGGCTTCACGTTGGTGCCGGTGCACCGGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGAc                                                            <  3:346037/135‑1 (MQ=255)
                                                                   cggggcggcCTACCCGGCGGTGGTCGAGGCGCGGGCCGCAGACCGCCTCGATCGGGCCTGGCGGGGGGGCTTCACGTGCGTTCAGGGGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGAc                                                            <  7:186347/133‑1 (MQ=255)
                                                                    ggggCGGACTACCCGGCGGTGGGCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTGCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGAc                                                            <  5:332804/133‑1 (MQ=255)
                                                                               cccGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCAcc                                                <  4:208545/136‑1 (MQ=255)
                                                                                            gAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCGTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTTCCGCGTGAGCAGGCCGCCCAGTTGACCTCCAAGTtcatc                                    >  2:14425/1‑135 (MQ=255)
                                                                                             aggcGCTGGCGGCAGACCGCCTCGTTCTGGCGTGGCTGGGGGGCTTCACGTTCGTTCAGTTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCAc                                  <  3:44974/133‑1 (MQ=255)
                                                                                                          aGACCGCCTCGATCGGGCCTGGCGGGGGGGCTTCACGTGCGTTCAGGGGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGAc                     <  3:3243/136‑1 (MQ=255)
                                                                                                             ccGCCTGGATCTGGCCTGGTTGGGGGGCTTCACGTTGGTGCCGGTGCACCTGATGCACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGt                   <  2:190576/135‑1 (MQ=255)
                                                                                                             ccGCCTCGAGCGGGCCGGGCTGGGGGGCTTCACGTTCGTTCAGGTGCCCCTGAAGGACCCGACGGCCACGCCACTGGTAGAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGt                   <  6:324346/135‑1 (MQ=255)
                                                                                                                  tCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAaga              >  5:136529/1‑135 (MQ=255)
                                                                                                                           cctggcTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAAGAGCCt          <  2:11075/130‑1 (MQ=255)
                                                                                                                           cctggcTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAAGAGCCt          <  5:330448/130‑1 (MQ=255)
                                                                                                                             tggctggGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAAGAGCCTGAGCGATc  >  4:209688/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                                                   
AGCGCAAGTTCAAGCCTTTGGGTGAGTACCTGGCCAAGCAACTGGGCATGGAGGTGAAGTTCGTACCCGTGGCGGACTACCCGGCGGTGGTCGAGTCGCTGGCCGCAGACCGCCTCGATCTGGCCTGGCTGGGGGGCTTCACGTTCGTTCAGGTGCACCTGAAGGACCCGACCGCCACGCCACTGGTACAGCGTGAGCAGGACGCACAGTTCACCTCCAAGTTCATCACCGCCAACCCCGACGTGAAGAGCCTGAGCGATC  >  ACB122/955712‑955972

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