Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,376 G→A 13.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,3760GA13.0% 87.2 / 4.5 41P14L (CCC→CTC) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (2/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGT  >  NC_000913/53235‑53510
                                                                                                                                             |                                                                                                                                      
ggggAATAAGGGCGGAGGGTGAGCGGCAAACCGACCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTTACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCggg                                                                                                                                       <  1:868868/143‑1 (MQ=255)
 gggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGGAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCCCCCACGCGGTGCAAGCTGGACGACTAAGTCCGGGCCAGCACCGGCgggg                                                                                                                                      >  2:91813/1‑142 (MQ=255)
  ggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACGAACCACCGCCTCGACCTGCCACTCACGCCGTGCAAGCTGGACAACTAAGTCCGGGACAAACCCGGCGGGcc                                                                                                                                     >  1:33517/1‑142 (MQ=255)
    aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCg                                                                                                                                   >  1:620029/1‑143 (MQ=255)
          ggCGGTGGGGGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGTACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                                                                                                             <  2:809602/143‑1 (MQ=255)
          ggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                                                                                                             <  2:381639/143‑1 (MQ=255)
              gAGGGTGAGCGGAAAACCGAGCATCGCTGCACGGTTGGTAAGGAGAGTGGCATCGGCAAAAACAACCAGTTCGAACGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGATTgg                                                                                                                         <  1:233259/143‑2 (MQ=255)
                   tGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCCATCCCGGCGCGCACGCCGGGGGTGATccca                                                                                                                    >  1:510384/1‑143 (MQ=255)
                         gCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCGATCCGGGCGGGCTCGCCGGGGGTGGTCACAACAAGt                                                                                                              >  1:195345/1‑143 (MQ=255)
                          cAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGATAGGAGGGAGGAATCACAAAACGtt                                                                                                             >  1:333229/1‑143 (MQ=255)
                           gagccGAGCATCGCGGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                                                                                                            <  2:715018/140‑1 (MQ=255)
                           caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCAGAAACAACCAGTTCGGCCGTGCTCTCGCGGGGTGCAAGCTGGGCAACTAAGTCCGGGCCGATCCCGGCGGGCTCGCCGGGAGTGGTCACAACACGtcg                                                                                                            >  1:650115/2‑141 (MQ=255)
                            aaCCGAGCATCGGTGTTCGGTTGGTAAGGGGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCTGGCTCGCCGGGAGTGATCACAACACGTTgg                                                                                                           <  2:304545/143‑1 (MQ=255)
                               ggggCAGCGCTGCCCGGTTGGTGTGAAGAGTGTCAACGCGGTCGGCAGCCAGTTCGGCCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtt                                                                                                        <  1:567722/140‑1 (MQ=255)
                                         tGCCCGGTGGGGAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTa                                                                                              <  1:361320/143‑1 (MQ=255)
                                            gctgTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGtt                                                                                           <  2:774992/140‑1 (MQ=255)
                                             cGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                          >  1:505811/1‑143 (MQ=255)
                                             cGGTTGGTAAGGAGAGTGGCATCGGAACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCGCGCCGGGAGAGAACACAACACGTTGGGTTTTAACCagtagtgg                                                                                          >  2:154299/1‑141 (MQ=255)
                                                 tGGGAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCa                                                                                      <  1:45413/143‑1 (MQ=255)
                                                  ggTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAg                                                                                     >  1:189446/1‑143 (MQ=255)
                                                    tAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                   <  2:189446/143‑1 (MQ=255)
                                                       ggAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAttt                                                                                >  1:498433/1‑143 (MQ=255)
                                                                     gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCGCGGCGGGCGCGGCGGGAGGGATCACAACACGTTGGGTGTTGACCATTAGTTGCTCAGGATTTTAACGTAGgcggta                                                                  >  2:579339/1‑140 (MQ=255)
                                                                           acaccCCGTTCGACCGGCCACGCACGCTGTGCAAGCGGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgt                                                            <  2:97311/143‑1 (MQ=255)
                                                                                cAGTTCTGCCGGCCTCTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                       <  2:727433/143‑1 (MQ=255)
                                                                                cAGTGCGACCGGCCACTCACGCGGGGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                       <  2:609737/143‑1 (MQ=255)
                                                                                 aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAGTTTAACGTAGGCGCTGGCACGTTGTTcc                                                      >  2:570294/1‑143 (MQ=255)
                                                                                  gTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                                                     <  1:789436/143‑1 (MQ=255)
                                                                                   ttCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCTCTAGCACGTTGCTCCTg                                                    >  2:827338/1‑143 (MQ=255)
                                                                                   gtCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTACTg                                                    >  1:344000/2‑143 (MQ=255)
                                                                                        ccGGCCACTCACGCTGTGCAAGCTGGACAACAAAGACCGGGACAAGCACGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATcc                                               <  1:300432/143‑1 (MQ=255)
                                                                                          ggCCACCCACGCTGCGCAAGCTGGACGACTAAGTCCGGGCGAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCCGGAGTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAg                                             <  2:407034/143‑1 (MQ=255)
                                                                                                        tgCGAGCGGGCCAACTATGTCGGGGCCCATCCCTGCGAGCTCGCCGGGGGTGAGCACAACACGTTGGGTTTTAACCACTAGTTGCTCACTATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCttcttc                               <  2:435663/143‑1 (MQ=255)
                                                                                                        tgCAAGCTGGACAACTACGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCttcttc                               <  1:264253/143‑1 (MQ=255)
                                                                                                        tgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCttcttc                               <  1:487500/143‑1 (MQ=255)
                                                                                                          cAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                             <  2:811786/143‑1 (MQ=255)
                                                                                                          cAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGATGCTTATTccc                             >  1:1182/1‑142 (MQ=255)
                                                                                                              ccggACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCGGCCAGCGTGCTGCTTCTTCCgggga                         >  1:166668/3‑139 (MQ=255)
                                                                                                               tGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGCGAGTGATCACAACACGTTGGGTTTTCACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCCGCTCGCTGCTTCTGCCGAGGAc                        >  2:662979/1‑143 (MQ=255)
                                                                                                                           tCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatc            >  1:697965/1‑143 (MQ=255)
                                                                                                                             cGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAGTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCGAGCTTGCTGCTTCTTCCTATATAGTACGGTTCATCAg          >  1:25876/1‑143 (MQ=255)
                                                                                                                                    ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCCCGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGGTCATCAGCCTGCgg   >  2:492774/2‑143 (MQ=255)
                                                                                                                                     tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATGAGTTGCTCAGGATTTTAACGTAGGCGCTGGCCCGTTGTTCCTGCATCCAGCTTGATGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGt  >  1:769245/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                      
GGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGT  >  NC_000913/53235‑53510

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: ‑