Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,727	G→A	A128V (GCC→GTC)	btuE ←	thioredoxin/glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,727	0	G	A	100.0%	166.1 / NA	47	A128V (GCC→GTC)	btuE	thioredoxin/glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/25);  total (22/25)

AATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTG  >  NC_000913/1793594‑1793858                                                                                                                                      |                                                                                                                                    aaTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGa                                                                                                                               >  1:15937/1‑140 (MQ=255)  aTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGaa                                                                                                                              >  1:515336/1‑140 (MQ=255)   tGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAAt                                                                                                                             >  2:808376/1‑140 (MQ=255)   tGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAAt                                                                                                                             >  1:1147502/1‑140 (MQ=255)          ttCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGctct                                                                                                                      <  1:187402/140‑1 (MQ=255)              ggCGTAATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTc                                                                                                                  >  1:1080299/1‑140 (MQ=255)               gCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTcc                                                                                                                 >  2:1082169/1‑140 (MQ=255)               gCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTcc                                                                                                                 >  1:303235/1‑140 (MQ=255)                   cATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGgcg                                                                                                             >  2:385924/1‑140 (MQ=255)                             gAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGt                                                                                                   <  1:1082169/140‑1 (MQ=255)                               aaaaCGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCg                                                                                                 >  2:40025/1‑140 (MQ=255)                                          tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGca                                                                                      <  2:282606/140‑1 (MQ=255)                                            aCTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaat                                                                                    <  2:515336/140‑1 (MQ=255)                                            aCTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGaaat                                                                                    <  2:303235/140‑1 (MQ=255)                                              ttttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatca                                                                                  <  1:438523/140‑1 (MQ=255)                                                ttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaat                                                                                >  1:518723/1‑140 (MQ=255)                                                   cGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAtttt                                                                             <  1:40025/140‑1 (MQ=255)                                                     tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGAACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTg                                                                           <  2:1080299/140‑1 (MQ=255)                                                       ccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCCGCAATCAATTTTTGAt                                                                         >  2:446927/1‑140 (MQ=255)                                                         tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATAc                                                                       <  2:51388/140‑1 (MQ=255)                                                            cAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGc                                                                    <  1:610529/140‑1 (MQ=255)                                                               ccAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACGATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGa                                                                 >  2:844331/1‑140 (MQ=255)                                                                 aGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                                                               <  1:385924/140‑1 (MQ=255)                                                                    aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCg                                                            <  2:15937/140‑1 (MQ=255)                                                                        tttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCt                                                        >  1:274727/1‑140 (MQ=255)                                                                         ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCtt                                                       <  1:808376/140‑1 (MQ=255)                                                                         ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCtt                                                       <  2:67246/140‑1 (MQ=255)                                                                         ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCCGCAATCAATTTTTGATACAGCGGATGGCGTCCtt                                                       <  1:446927/140‑1 (MQ=255)                                                                          tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc                                                      <  2:971772/140‑1 (MQ=255)                                                                          tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTc                                                      <  2:484006/140‑1 (MQ=255)                                                                               aTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCAt                                                 <  1:844331/140‑1 (MQ=255)                                                                                   cATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAAc                                             >  1:1066855/1‑140 (MQ=255)                                                                                        aaTATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCaa                                        <  2:1147502/140‑1 (MQ=255)                                                                                        aaTATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCaa                                        <  2:518723/140‑1 (MQ=255)                                                                                                          cGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATc                                      <  2:204036/124‑1 (MQ=255)                                                                                                          cGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATc                                      >  1:204036/1‑124 (MQ=255)                                                                                                          cGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGa                      <  2:114006/140‑1 (MQ=255)                                                                                                          cGGTGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTACTTCGCCATTCACGTCAATCTTACTGAACATCGGGa                      <  2:531816/140‑1 (MQ=255)                                                                                                             tGCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACg                   >  1:821281/1‑140 (MQ=255)                                                                                                              gCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATc                          <  2:342059/132‑1 (MQ=255)                                                                                                              gCACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATc                          >  1:342059/1‑132 (MQ=255)                                                                                                               cACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTc                 >  2:791979/1‑140 (MQ=255)                                                                                                               cACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTc                 >  2:975935/1‑140 (MQ=255)                                                                                                               cACGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTc                 >  2:1186369/1‑140 (MQ=255)                                                                                                                 cGGCCTTTGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCAc               <  1:791979/140‑1 (MQ=255)                                                                                                                       ttGCTGACCATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCa         <  1:975935/140‑1 (MQ=255)                                                                                                                              ccATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATgtggtg  >  2:133089/1‑140 (MQ=255)                                                                                                                              ccATACGGACATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATgtggtg  >  1:934668/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                    AATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTG  >  NC_000913/1793594‑1793858

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