Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,592	A→G	V173A (GTG→GCG)	btuE ←	thioredoxin/glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,592	0	A	G	100.0%	164.9 / NA	48	V173A (GTG→GCG)	btuE	thioredoxin/glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

TTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGG  >  NC_000913/1793455‑1793727                                                                                                                                          |                                                                                                                                        ttCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATcgca                                                                                                                                       >  1:317264/1‑140 (MQ=255) ttCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAAC‑GCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCaa                                                                                                                                      >  2:652332/1‑140 (MQ=255)      ccccACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATggg                                                                                                                                  >  1:10812/1‑140 (MQ=255)           ccaGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGGAGATTCCGCAACATCTTTTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCtt                                                                                                                             >  1:965604/1‑140 (MQ=255)            caGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAAC‑GCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTcc                                                                                                                           <  1:652332/140‑1 (MQ=255)               gTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCgg                                                                                                                         >  1:21264/1‑140 (MQ=255)                    ggATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCTACATTCTGGAACTGCCACACAATAGAAAGTTATTTACCCACGCCAGTTTAATGCTTTTCAACTCAATGGGAACTTCCGGCGTcc                                                                                                                    >  2:969891/1‑139 (MQ=255)                        ctcCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATc                                                                                                                <  2:21264/140‑1 (MQ=255)                                   gAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAAc                                                                                                     >  1:667536/1‑140 (MQ=255)                                         cGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGa                                                                                               <  2:160351/140‑1 (MQ=255)                                                 aGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACtttt                                                                                       <  1:969891/140‑1 (MQ=255)                                                     gCCCCAGGCGGGGAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGt                                                                                   >  1:981108/1‑140 (MQ=255)                                                          aGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTg                                                                              <  2:317264/140‑1 (MQ=255)                                                          aGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTg                                                                              <  2:667536/140‑1 (MQ=255)                                                            gCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGcc                                                                            <  2:10812/140‑1 (MQ=255)                                                                gTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAAcc                                                                        >  1:987253/1‑140 (MQ=255)                                                                          cAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAAttttt                                                              >  2:588677/1‑140 (MQ=255)                                                                                      aCTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATa                                                  <  2:687983/140‑1 (MQ=255)                                                                                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGTTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAat                                              >  1:489350/1‑140 (MQ=255)                                                                                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATaaagt                                              >  2:998594/1‑138 (MQ=255)                                                                                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAat                                              >  2:202942/1‑140 (MQ=255)                                                                                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAat                                              <  1:372085/140‑1 (MQ=255)                                                                                          cATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAat                                              >  2:372085/1‑140 (MQ=255)                                                                                            tCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAatat                                            <  2:731533/140‑1 (MQ=255)                                                                                            tCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAatat                                            <  2:649317/140‑1 (MQ=255)                                                                                                    gACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATcc                                       >  2:469701/1‑137 (MQ=255)                                                                                                    gACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATcc                                       <  1:469701/137‑1 (MQ=255)                                                                                                      cATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCTGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTa                                  <  2:965604/140‑1 (MQ=255)                                                                                                       attattTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTAc                                 >  1:886864/1‑140 (MQ=255)                                                                                                         tattTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAg                               >  2:524807/1‑140 (MQ=255)                                                                                                         tattTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAg                               >  2:172687/1‑140 (MQ=255)                                                                                                         tattTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAg                               >  2:286818/1‑140 (MQ=255)                                                                                                          attTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                              >  1:21853/1‑140 (MQ=255)                                                                                                          attTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                              >  2:572307/1‑140 (MQ=255)                                                                                                          attTTGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                              >  2:1069122/1‑140 (MQ=255)                                                                                                           ttttGCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCg                             <  1:572307/140‑1 (MQ=255)                                                                                                               gCCAACGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGc                         <  2:21853/140‑1 (MQ=255)                                                                                                                ccAACGCCAGTTTAATGTTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCa                        <  2:489350/140‑1 (MQ=255)                                                                                                                    cGCCAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGc                    <  1:998594/140‑1 (MQ=255)                                                                                                                      ccAGTTTAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCt                  <  1:202942/140‑1 (MQ=255)                                                                                                                          tttAATGCTTTCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGAC‑TTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCt             <  2:205541/140‑1 (MQ=255)                                                                                                                          tttAATGCTTTCCATCGCAATGGGACCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGc              >  2:231260/1‑140 (MQ=255)                                                                                                                                  tttCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATa      <  1:524807/140‑1 (MQ=255)                                                                                                                                   ttCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATa      >  2:571791/1‑139 (MQ=255)                                                                                                                                   ttCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATa      <  1:571791/139‑1 (MQ=255)                                                                                                                                    tCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGc                    >  2:895446/1‑124 (MQ=255)                                                                                                                                    tCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGc                    <  1:895446/124‑1 (MQ=255)                                                                                                                                    tCCATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCAATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACg    >  2:243878/1‑140 (MQ=255)                                                                                                                                     ccATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACgg   <  1:1069122/140‑1 (MQ=255)                                                                                                                                      cATCGCAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACggg  >  2:503174/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        TTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGG  >  NC_000913/1793455‑1793727

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