Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,437,591	T→C	A87A (GCA→GCG)	pdxB ←	erythronate‑4‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,437,591	0	T	C	100.0%	122.7 / NA	41	A87A (GCA→GCG)	pdxB	erythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/20);  total (21/20)

GTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGTGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACTTTCGTGACCGAACGCAC  >  NC_000913/2437474‑2437725                                                                                                                      |                                                                                                                                        gTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCg                                                                                                                  <  1:68205/140‑1 (MQ=255)  ttACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAACCCAATT‑CCCGc                                                                                                                 <  1:82845/140‑1 (MQ=255)   tACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGcc                                                                                                                >  1:828873/1‑140 (MQ=255)           cccACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑c                                                                                                                     <  2:825754/127‑1 (MQ=255)           cccACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑c                                                                                                                     >  1:825754/1‑127 (MQ=255)           cccACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCaa                                                                                                        <  1:813310/140‑1 (MQ=255)             cACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAAcc                                                                                                      >  1:559445/1‑140 (MQ=255)               cgtcccccccCGTACGGTCGTAAAGTGAAAATCCATCGCGTCCGGCAAGCATCAGCAGGGAGGAGAAAACATATCCCCCCCCCGCAACCGCTTTACGGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCAt                                                                                                    <  2:828873/136‑1 (MQ=255)                    cccACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCCCCGCAATCGCATTACAGCCAGGCGCAGCGGAAAACCCAATT‑CCCGCCTGCTTCAACCATGCTTc                                                                                               <  1:817667/140‑1 (MQ=255)                         cGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCGGCGGAAAAACCAATG‑CCCGCCTGCTTCAACCATGCTTCATCGa                                                                                          <  2:555574/140‑1 (MQ=255)                          gTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGa                                                                                          <  2:364906/139‑1 (MQ=255)                          gTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGa                                                                                          >  1:364906/1‑139 (MQ=255)                          gTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGAc                                                                                         >  2:653992/1‑140 (MQ=255)                          gTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGAc                                                                                         >  1:472290/1‑140 (MQ=255)                           tACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACa                                                                                        <  2:559445/140‑1 (MQ=255)                             cGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCAGCCTGCTTCAACCATGCTTCATCGACATg                                                                                      >  2:93959/1‑140 (MQ=255)                              ggTCGTACAGTAAAAACCCACCGCTTTCGGCAAGCATCAGCAGGGAGGAAAAAACATATCCCCCCCCCGAAATCCCATTACACCCAGGCGCAG‑GGAAAAACCCATTCCCCGCCTGCTTCAACCATGCTTCATCGACATgg                                                                                     <  2:526245/140‑1 (MQ=255)                                tCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCAGCCTGCTTCAACCATGCTTCATCGACATggtc                                                                                   <  1:93959/140‑1 (MQ=255)                                  gTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATggtcgg                                                                                 <  2:664418/140‑1 (MQ=255)                                   tACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCa                                                                                              >  1:284588/1‑126 (MQ=255)                                   tACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCa                                                                                              <  2:284588/126‑1 (MQ=255)                                      aGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCCCCCCCGCAATCGCATTACAGCCAGGCGCAGCGGAAAACCCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTccc                                                                             <  1:591301/140‑1 (MQ=255)                                      aGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTccc                                                                             <  2:326765/140‑1 (MQ=255)                                        tGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCg                                                                           >  1:470838/1‑140 (MQ=255)                                          aaaaTCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCgca                                                                         >  1:504275/1‑140 (MQ=255)                                          aaaaTCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCgca                                                                         >  2:609852/1‑140 (MQ=255)                                          aaaaTCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCg                                                                           <  1:784225/138‑1 (MQ=255)                                          aaaaTCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCg                                                                           >  2:784225/1‑138 (MQ=255)                                                  tCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCgcagtggcagt                                                                 >  1:24497/1‑140 (MQ=255)                                                  ccgcgTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCgcagtggcagt                                                                 <  1:649335/139‑1 (MQ=255)                                                         cGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCCCCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCcaaca                                                          <  2:505650/140‑1 (MQ=255)                                                         cGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCcaaca                                                          <  1:653992/140‑1 (MQ=255)                                                               gCATCAGCAGGGAGGAGAAAACATATTCCACCCCCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTa                                                    <  2:472290/140‑1 (MQ=255)                                                                    agcagGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATggg                                               >  1:599003/1‑140 (MQ=255)                                                                          gaggagAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTcc                                         >  2:680616/1‑140 (MQ=255)                                                                                   aCATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCaa                                >  2:390677/1‑140 (MQ=255)                                                                                   aCATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCaa                                <  1:680616/140‑1 (MQ=255)                                                                                   aCATATTCCACCACCGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCaa                                >  1:805479/1‑140 (MQ=255)                                                                                             ccaccGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGattcattc                      >  2:363969/1‑140 (MQ=255)                                                                                                 cGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACtt                  >  2:213575/1‑140 (MQ=255)                                                                                                 cGCAATCGCATTACAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACtt                  >  1:215459/1‑140 (MQ=255)                                                                                                             aCAGCCAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACTTTCGTGACCGAAc      <  2:506208/140‑1 (MQ=255)                                                                                                                 ccAGGCGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACTTTCGTGACCGAACGCAc  >  2:403757/1‑140 (MQ=255)                                                                                                                      |                                                                                                                                        GTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGCAGGGAGGAGAAAACATATTCCACCACCGCAATCGCATTACAGCCAGGTGCAGCGGAAAAACCAATT‑CCCGCCTGCTTCAACCATGCTTCATCGACATGGTCGGTCCCCGCAGTGGCAGTGCCAACAAATTTAATGGGTTTTCCTGCCAGCAAAGATTCATTCACTTTCGTGACCGAACGCAC  >  NC_000913/2437474‑2437725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑