Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 638,614 G→T 100% T14K (ACA→AAA)  thiD ← phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730638,6140GT100.0% 110.0 / NA 37T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/18);  total (19/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAT  >  CP000730/638477‑638750
                                                                                                                                         |                                                                                                                                        
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTtttg                                                                                                                                        >  1:260932/1‑140 (MQ=255)
        ggTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                <  1:18991/140‑1 (MQ=255)
             tGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                           >  2:311575/1‑140 (MQ=255)
             tGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                           >  2:71138/1‑140 (MQ=255)
                  tAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGCACGGGCaactagt                                                                                                                      >  2:139227/1‑135 (MQ=255)
                  tAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGAAACGTTTTCAAATCTGCTTGCATACATACGCCAGCACTTTTGTCCGAACCGGCAATTGtt                                                                                                                      >  2:297459/1‑140 (MQ=255)
                            tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTCTTAAAActtcct                                                                                                            >  2:103942/1‑140 (MQ=255)
                                cATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaa                                                                                                        <  2:86732/140‑1 (MQ=255)
                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAACACTACAATCGTGTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCa                                                                                                   <  1:314221/140‑1 (MQ=255)
                                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTACTTCTTTGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTg                                                                                                <  1:151880/140‑1 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                                                              <  2:294674/140‑1 (MQ=255)
                                            tCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCATCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                                                                            >  1:2955/1‑140 (MQ=255)
                                              cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTc                                                                                          <  1:104912/140‑1 (MQ=255)
                                                tAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTCTCAAAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                        >  1:15972/1‑140 (MQ=255)
                                                       gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCGAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAACGCCATTGAGCTTCAATCCTACa                                                                                 <  2:76899/140‑1 (MQ=255)
                                                        aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACat                                                                                >  2:254701/1‑140 (MQ=255)
                                                         tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGACCCGGCAATTGTTACAACTTTCTTTAAAGCCATCGAGCTTCACTcctaccta                                                                               >  1:291597/1‑140 (MQ=255)
                                                             aGTTAAACCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                           <  2:31832/140‑1 (MQ=255)
                                                                 aaaGCGACCAGGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                       <  2:226577/140‑1 (MQ=255)
                                                                            gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCATGCATACATGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                            >  2:200403/1‑140 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                           <  2:315237/140‑1 (MQ=255)
                                                                                 aCGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatcat                                                       <  1:99273/140‑1 (MQ=255)
                                                                                  cGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtt                                                      <  2:226444/140‑1 (MQ=255)
                                                                                          aTTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCt                                              >  1:257075/1‑140 (MQ=255)
                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAAAATTGTATTCATCATATCATTTTTAACCTAATTg                                         >  1:288693/1‑140 (MQ=255)
                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGTAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaa                                       >  1:214612/1‑140 (MQ=255)
                                                                                                   aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaa                                     >  1:147219/1‑140 (MQ=255)
                                                                                                    cGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTCAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTaaaaaa                                    >  1:192737/1‑140 (MQ=255)
                                                                                                    cGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaaa                                    >  1:185077/1‑140 (MQ=255)
                                                                                                     gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAat                                   <  1:279932/140‑1 (MQ=255)
                                                                                                     gTTTTCAAATCTGCGTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATGGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCGTACATAATAATATTGTATTCATCATATCATTTTTAACCTAAGTGAAAAat                                   <  2:236161/140‑1 (MQ=255)
                                                                                                        ttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                <  1:123888/140‑1 (MQ=255)
                                                                                                                    tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatt                    >  2:19690/1‑139 (MQ=255)
                                                                                                                          ccTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAACACTTTCTTTAAAGCCATTCAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgat              <  2:25697/140‑1 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGCtt          >  1:65765/1‑140 (MQ=255)
                                                                                                                                ccAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt        >  2:291558/1‑140 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTGTCGAGCCATTGAGATTCACTCCTACATTAGAATATCCTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaa    <  2:122119/140‑1 (MQ=255)
                                                                                                                                      cTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTACTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAt  <  2:186761/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAT  >  CP000730/638477‑638750

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