Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,376 G→A 12.0% P14L (CCC→CTC)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3760GA12.0% 92.8 / 6.2 45P14L (CCC→CTC) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (3/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

GGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53234‑53513
                                                                                                                                              |                                                                                                                                         
gggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATGGCTGCTGGGGAGGTAAGGAGAGTGGCATCGGCACACACAACCAGTTCGACCGGCCCCTCACGCTGGGCAAGCTGGACAACTACGACCGGGCCCATCAcggcgg                                                                                                                                           >  2:272264/1‑143 (MQ=255)
 ggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTGGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCcgggcggg                                                                                                                                          >  2:649323/1‑143 (MQ=255)
        gaGGGCGGGGGGTGAGCGGCAAACCGAGCAGCGCTGCCCGGTTGGTAAGGAGAGTGGCGTCGGCGCAAACAGCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                                                                   <  2:219925/142‑1 (MQ=255)
           ggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGGGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                                                                                                                <  1:455807/143‑1 (MQ=255)
              ggAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTGCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGATTg                                                                                                                             <  1:409648/143‑1 (MQ=255)
                aGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAGTCCCGGCGGGCTCGCCGGGAGTGAt                                                                                                                           >  1:627763/1‑143 (MQ=255)
                   gTGAGCGCCAAACCGAGCACCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCACGCTGGAAAACTAAGTCCGGTCCAATCCCGGCGGGCTCGCCGGGAGTGATcac                                                                                                                        <  1:169869/143‑1 (MQ=255)
                      aGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCCGCGGCCCCGCCGAGGGTGagcgggac                                                                                                                     >  2:772478/1‑136 (MQ=255)
                         ggCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTGCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTAAGAGCCAATCCCGGCGGGCTCGCGGGAGGGGGACACAACACg                                                                                                                  >  2:690473/1‑143 (MQ=255)
                          gCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGGGATGACAACACGt                                                                                                                 >  1:195981/1‑143 (MQ=255)
                              aCCGAGCATCGCTGCCCGGTTCGTAAACCGAATAAGGTTTTGACGAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                                                             <  1:842939/143‑1 (MQ=255)
                               ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                                                                                                            <  1:516943/143‑1 (MQ=255)
                                cgcgCAGCGCGGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCCCTCACGCTGTGCAGGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtt                                                                                                           <  1:399488/140‑1 (MQ=255)
                                    cATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                                                       <  1:353342/143‑1 (MQ=255)
                                     aTCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGCGATCCCAACCCGTTGGGTCTTCAc                                                                                                      >  2:11086/1‑143 (MQ=255)
                                      tCGCTGCCCGGGTAGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc                                                                                                     <  1:316690/143‑1 (MQ=255)
                                            cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCGATCCCGGCGGGCTCGCCGGGAGTGATCACGACACGTTGGGTTTTAACCATTAGt                                                                                               >  1:32559/1‑143 (MQ=255)
                                            cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGATCTCTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCACTAGt                                                                                               >  2:636476/1‑143 (MQ=255)
                                             ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGATGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGGCGGGAGTGATCACAACACGTTGGGTTTTAGCCATTAGtt                                                                                              >  2:82166/1‑143 (MQ=255)
                                                 ttcgTAAGGAGAGTGGCATCGCCACAAACACCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTc                                                                                          <  1:269705/140‑1 (MQ=255)
                                                  cggcaaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCa                                                                                         <  1:248928/139‑1 (MQ=255)
                                                     tAAGTAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                      <  1:373434/143‑1 (MQ=255)
                                                        ggAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGCGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCCGGCttt                                                                                   >  1:340252/1‑143 (MQ=255)
                                                         gagaGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGGCAACTAAGTCCGGGCCAATCCCGGCGGGCCCGCCGGGAGGGAACACAACACGTCGGGCCTTAACCGCGAGTTGCTCAggcatct                                                                                  >  2:3305/1‑138 (MQ=255)
                                                                            acaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTGGGCGCTGGCACgt                                                               >  1:399192/1‑143 (MQ=255)
                                                                                 taGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                          <  1:855864/142‑1 (MQ=255)
                                                                                   gTTCGACCGGCCACTCAGGCGGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                        <  2:859231/143‑1 (MQ=255)
                                                                                         ccGGCCACTCACGCTGTGCAAGCTGGGCAACTAAGTCCGGGCCAGTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCAAGCATcc                                                  <  2:834687/143‑1 (MQ=255)
                                                                                         ccGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCCCAACACGTTGGGTGTGAACCATGAGTTGCTCACGATATTAACGGAGGCGCGGCCACGTTGTTccgcgctca                                                  >  2:55740/1‑136 (MQ=255)
                                                                                           ggCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTCTAACCATTAGTTGGCCAGGAGTTTGACGTACGCGCCGGCACGTTGTTCGTGCATcccg                                                >  2:99483/1‑141 (MQ=255)
                                                                                            gCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCCTTCGTTGCGTAGGATTTTAACGTAGGCGCCAGCACGTTGTTCCTGCATCCAGc                                               >  2:307736/1‑143 (MQ=255)
                                                                                               aCTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTGCTGCCTGCAGCTtg                                            >  1:116504/1‑143 (MQ=255)
                                                                                                 tCACGCTGTGCAAGCTGGACAACTAAGTCCGTGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgct                                          <  1:504253/143‑1 (MQ=255)
                                                                                                     gCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGctt                                      >  1:803634/1‑143 (MQ=255)
                                                                                                            aaGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCga                               <  1:101636/143‑1 (MQ=255)
                                                                                                            aaGATGGACAACTAAGTCCGGGCCAATCACGGCGCGCTCGCCGGGAGTGATCACAACACGTTGGGTCTTAACCATTAGTTGACCAGGATTTTAACGTCGGCGCTGGCACGCTGTTCCTCCATGCAGCTTACCGCTACTTACgc                               >  2:594948/1‑142 (MQ=255)
                                                                                                                 ggACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGCTGGGTTTTAACCATTGGTTGATCAGGCTTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCGTCCGAGAACt                          >  2:239309/1‑143 (MQ=255)
                                                                                                                   gcaggtAAGTCCGGGCCAACCCCGGCGGGCTCGCCGGGAGTGCTCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTa                        <  1:502168/138‑1 (MQ=255)
                                                                                                                       cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCCCCACACGTTGAGTTTTAACCAGTAAGCGCCCAGGATCGTAACGAAGGCGCTGGCACGTTGGTCCTGCCTCCAGCTTGCTGCTTCTTACGCGAACTATCGGt                    >  2:159175/1‑143 (MQ=255)
                                                                                                                       cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGt                    >  2:644196/1‑143 (MQ=255)
                                                                                                                             ccGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGACCTTACGGTtcatca              <  1:750924/143‑1 (MQ=255)
                                                                                                                               gagCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGc            <  1:150252/141‑1 (MQ=255)
                                                                                                                                    actCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg       >  2:569601/3‑143 (MQ=255)
                                                                                                                                    actCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCACCTTGCTGCTTCTTCCGAGAACTTGCGGTTCATCCGCAGGCg       >  1:807448/3‑143 (MQ=255)
                                                                                                                                     aTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCCTGcag      <  2:364431/143‑3 (MQ=255)
                                                                                                                                      tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACGTACGGTTCATCAGCATGCGGt     >  2:390182/1‑143 (MQ=255)
                                                                                                                                        ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCACTAGTTGCTCAGGATTTTAACGCAGGCGCTGGCACGTTGTTCCTCCATCCAGCTTGCTGCTTCTTCCGAGCACGTACGGTTCATCAGCGTGCgggaa   >  2:674079/1‑140 (MQ=255)
                                                                                                                                         cggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCGTCAGCATGCGGtgtg  >  1:539245/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
GGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53234‑53513

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