Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,474,315	(C)6→5	coding (95/861 nt)	pdxI →	pyridoxine 4‑dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,474,310	0	C	.	100.0%	110.2 / NA	27	coding (90/861 nt)	pdxI	pyridoxine 4‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/16);  total (11/16)

ATTTTTCTTAATCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGAAGCGCTTTATCC  >  NC_000913/1474184‑1474441                                                                                                                               |                                                                                                                                    aTTTTTCTTAATCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGcccccac                                                                                                                               >  2:102259/1‑131 (MQ=255)   ttttCTTAATCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGcccccacga                                                                                                                             >  2:379414/1‑129 (MQ=255)      tCTTAATCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGAt                                                                                                                         >  1:332765/1‑133 (MQ=255)            tCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCAc                                                                                                                   <  2:111591/133‑1 (MQ=255)            tCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCAc                                                                                                                   <  2:101826/133‑1 (MQ=255)             ccGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACg                                                                                                                  <  2:387642/133‑1 (MQ=255)                 tttATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGc                                                                                                              <  1:379414/133‑1 (MQ=255)                       aCAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGATAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAAc                                                                                                        >  2:452491/1‑133 (MQ=255)                           ggcagggTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGTCATGTCGCTTTAACCGtg                                                                                                    >  1:244716/1‑133 (MQ=255)                                ggTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCg                                                                                               >  1:642647/1‑133 (MQ=255)                                     gATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAgg                                                                                          >  2:620195/1‑133 (MQ=255)                                         agcagcAATACATTTGCTGTCGATACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCt                                                                                      <  1:293836/133‑1 (MQ=255)                                              caATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAt                                                                                 <  2:642647/133‑1 (MQ=255)                                                   cATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                                                                            <  1:33975/133‑1 (MQ=255)                                                   cATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                                                                            >  2:33975/1‑133 (MQ=255)                                                   cATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                                                                            <  2:329685/133‑1 (MQ=255)                                                         ctctCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAAt                                                                      <  1:590860/133‑1 (MQ=255)                                                          tctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATc                                                                     <  2:332765/133‑1 (MQ=255)                                                              ggTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCatat                                                                 >  1:362020/1‑133 (MQ=255)                                                                        ccGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc                                                       <  1:620195/133‑1 (MQ=255)                                                                        ccGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc                                                       <  1:102259/133‑1 (MQ=255)                                                                          gTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa                                                     >  1:635558/1‑133 (MQ=255)                                                                                            ggCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGt                                   >  1:98852/1‑133 (MQ=255)                                                                                                               gTCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCcgc                >  1:384494/1‑133 (MQ=255)                                                                                                                tCCTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCcgcg               >  2:509902/1‑133 (MQ=255)                                                                                                                 ccTGGAGTTTTTGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGa              <  2:11331/133‑1 (MQ=255)                                                                                                                  cTGGAGTTTTTGG‑CCCCCACGAGACCGCCACGACGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGaa             <  2:190514/133‑1 (MQ=255)                                                                                                                  cTGGAGTTTTCGG‑CCTCCACGAGATCGCTACGTCGCTATAAACGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGaa             <  2:244716/133‑1 (MQ=255)                                                                                                                        tttttGG‑CCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGAAGCGCtt       <  2:384494/133‑1 (MQ=255)                                                                                                                          tttGG‑CCCCCACGAGATCGCCA‑‑‑CGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGAAGCGCTTTATcc  >  1:539702/1‑133 (MQ=255)                                                                                                                               |                                                                                                                                    ATTTTTCTTAATCCGTTTTATTACAGGGCAGGGTGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTATGGTCCGCACGTCACCAATCAGATTATCCGCGAAGCGCTTTATCC  >  NC_000913/1474184‑1474441

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