Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,370 C→G 11.3% intergenic (+27/+16) crr → / ← pdxK Enzyme IIA(Glc)/pyridoxal kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,3700CG11.3% 114.4 / 5.2 48intergenic (+27/+16)crr/pdxKEnzyme IIA(Glc)/pyridoxal kinase I
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base G (3/3);  total (22/26)
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.82e-01

CCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAG  >  NC_000913/2536237‑2536512
                                                                                                                                     |                                                                                                                                              
cccTGACTCTGGTTTTTATTTTCAAAAAGGACGAAATCAAAGAACGGAACAAACGGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGcc                                                                                                                                       <  1:86084/143‑1 (MQ=255)
  cTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAGGGAGTAATTCTGGCCGCAGTGAAAAAGGGCGCCCAGCGGCGGCag                                                                                                                                     >  1:325764/1‑142 (MQ=255)
    gACTCCGGTTGTTAGCTCCAACGTGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttt                                                                                                                                   <  2:892706/143‑1 (MQ=255)
       tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAATCGTGGGTGAAACCCCGGTTATCCGCATCACGCCGTAATTCTTGCCGCAGTGAACAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                >  1:179966/1‑143 (MQ=255)
       tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGACGTAATTCTTGCCGCAGTGAAAAATTGCGCCCATCGGCGCCAtttttt                                                                                                                                >  1:392912/1‑143 (MQ=255)
       tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGGAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                <  1:935823/143‑1 (MQ=255)
          tgttgttATCTCCCACATGGACGAAATCAAAGACCGGATCAAACTGTCCGGCAGCGTAACCGTGGGTGAATCCCCGGTTAGCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCAATCGGCGCCATTTTTTTAt                                                                                                                             <  1:752400/142‑1 (MQ=255)
               ttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATAGCTCCCATCGGAGCCATTTTTTTATGCTTc                                                                                                                        >  2:419293/1‑143 (MQ=255)
                  gctcCAACATGGACGAAATCAAAGAACTGATCAGACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGc                                                                                                                     <  2:161537/142‑1 (MQ=255)
                      cAACATGGACGAAATCAAAGAACTGATCAAACTGTACGGTAGCGTAACCGTGGGTGAAACCCCAGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGATCCCGCCAgc                                                                                                                 >  2:997085/1‑143 (MQ=255)
                          acggACTAAATCAAATAACTGATCAGACTGTCCGGTAGCGCAACCGTGGGTGAAACCCCGGTTATCCGCATCAATAAGTCATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcggcg                                                                                                             <  2:599971/141‑1 (MQ=255)
                                  aaTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCCTCGGCCCCATTTTTTTATGCTTCCGGCAGCgacggccacaac                                                                                                     >  2:486989/1‑137 (MQ=255)
                                     tcgagagCTGATGATGCTGTCCTGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTACTCCTTTCCGTAGGGTAGAATGGCGCCCATCGGCGCGATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAAt                                                                                                  <  2:698366/140‑1 (MQ=255)
                                     cAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGACACCCCGGTTCTCCGCATCATGAAGTAATTCTTGCCGCAGTGAAAAATGACGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCACAATCact                                                                                                  >  1:741112/1‑141 (MQ=255)
                                      aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGACACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAGATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCCGGAGAATCAAtt                                                                                                 >  1:640603/1‑143 (MQ=255)
                                      aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAAtt                                                                                                 <  1:131877/143‑1 (MQ=255)
                                         gAACTGATCTAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGACAAAATCAATTCAt                                                                                              <  2:119293/143‑1 (MQ=255)
                                             tGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCCGCGGCGGCAAAATCCATTCATCGct                                                                                          >  1:956181/1‑143 (MQ=255)
                                               aTCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCCTTCTTTTATGCACACGCCAGCGGCGCCACAACCAATCCATCGctct                                                                                        >  2:23514/1‑143 (MQ=255)
                                                 cAAACTGTCCGCTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCa                                                                                      <  1:17844/143‑1 (MQ=255)
                                                   aaCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtg                                                                                    <  1:986795/143‑1 (MQ=255)
                                                      tGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgctg                                                                                 >  2:478715/1‑143 (MQ=255)
                                                        tCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATgctgct                                                                               >  1:555054/1‑143 (MQ=255)
                                                         ccGGTGGCGGAACCGTGTGTGAAACCGCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGctgctg                                                                              <  2:628592/143‑1 (MQ=255)
                                                            gTAGGGTAACCGTGGGTGAAACTGCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt                                                                           <  2:855947/143‑1 (MQ=255)
                                                                 gTAGCCGTGGGTGAGCCCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgc                                                                      <  2:374814/143‑1 (MQ=255)
                                                                    aCCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATCGGCGCCGTTTTTTTATGCTTCCGCCAGCGGCGGCACAATCAGTTCATCGCTCTCCTGCCGCTGTGTGTAACGCa                                                                   >  1:564234/1‑143 (MQ=255)
                                                                            tGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTcc                                                           <  1:407035/143‑1 (MQ=255)
                                                                               aaCCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGt                                                        <  1:901862/143‑1 (MQ=255)
                                                                                   ccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTCcgc                                                    >  2:123537/1‑143 (MQ=255)
                                                                                              cATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCCTAGCTGTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCt                                         >  1:772445/1‑143 (MQ=255)
                                                                                              cATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTCGCGCCTCCCTTCCAGTACGCGCAACCCCGCt                                         >  2:547129/1‑143 (MQ=255)
                                                                                                   cgaagTAATTCTTGCCGCAGTGCAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTGGCGCATCACTTCCAGTACGCGCAACCTCGCTTGGTg                                    <  2:991394/142‑1 (MQ=255)
                                                                                                       gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAACATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCGCt                                >  1:584485/1‑143 (MQ=255)
                                                                                                          aTTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTCCTGGGTGTAGCGCATCACTTCCAGCACACGCAACCCCTCTCGGTggtctacc                             >  2:329060/1‑136 (MQ=255)
                                                                                                           ttCTTGCCGCAGGGAAAAATGGCGCCGATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCCATGCATCGCTCGCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCAGCGCAt                            <  2:632099/143‑1 (MQ=255)
                                                                                                              ttGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCGAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                         >  1:571537/1‑143 (MQ=255)
                                                                                                               tGCCGTAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGt                        <  2:392912/143‑1 (MQ=255)
                                                                                                                 ccGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTa                      >  1:799964/1‑143 (MQ=255)
                                                                                                                        aaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCtt               >  2:776114/2‑143 (MQ=255)
                                                                                                                          aaaaTGGCGCCCATCGGCGCCATTTTTTGATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTcc             <  2:131568/143‑1 (MQ=255)
                                                                                                                            aaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCt           >  1:996269/1‑143 (MQ=255)
                                                                                                                              tGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCttt         <  2:612683/143‑1 (MQ=255)
                                                                                                                              tGGCGCCCATCGGCGCCATTGTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATGCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCttt         <  1:209154/143‑1 (MQ=255)
                                                                                                                               ggCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc        <  1:151397/143‑1 (MQ=255)
                                                                                                                                gcgcCGATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCACTTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCa       <  2:741112/143‑1 (MQ=255)
                                                                                                                                gcgcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCGTCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCa       <  1:850320/143‑1 (MQ=255)
                                                                                                                                gcgcCCATCGGCGCCATTTGTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCa       <  2:253499/143‑1 (MQ=255)
                                                                                                                                    ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCGTTGCAGCaa   >  1:898271/1‑143 (MQ=255)
                                                                                                                                     cATCGGCGCCATTTTTTTATGCGTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAg  <  1:784414/143‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                              
CCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAG  >  NC_000913/2536237‑2536512

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