Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,379 T→C 13.0% E13G (GAG→GGG)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3790TC13.0% 89.8 / 6.1 46E13G (GAG→GGG) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (5/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

AAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53242‑53513
                                                                                                                                         |                                                                                                                                      
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)
                                                                                                                                         |                                                                                                                                      
AAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53242‑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: