Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,600,609	(C)5→6	coding (132/471 nt)	bcp →	peroxiredoxin, thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,600,604	1	.	C	100.0%	90.6 / NA	27	T43P (ACC→CCC)	bcp	peroxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (13/14);  total (13/14)

AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑CCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600469‑2600732                                                                                                                                         |                                                                                                                                 aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc                                                                                                                               >  1:140316/1‑139 (MQ=255)      aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccccccgg                                                                                                                          >  1:371032/1‑136 (MQ=255)         aaTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTg                                                                                                                        <  2:140316/139‑1 (MQ=255)          aTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGt                                                                                                                       >  1:818767/1‑139 (MQ=255)             aaTCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTAcc                                                                                                                    >  2:174778/1‑139 (MQ=255)                 cACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc                                                                                                                               >  2:903708/1‑123 (MQ=255)                 cACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc                                                                                                                               <  1:903708/123‑1 (MQ=255)                             gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTa                                                                                                                 <  1:765373/126‑1 (MQ=255)                             gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTa                                                                                                                 >  2:765373/1‑126 (MQ=255)                              tGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCtt                                                                                                   >  1:1032647/1‑139 (MQ=255)                                   tCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                              <  2:371032/139‑1 (MQ=255)                                        ccGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAAc                                                                                         >  1:823645/1‑139 (MQ=255)                                                   cTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGtt                                                                              >  1:871808/1‑139 (MQ=255)                                                           aTCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAg                                                                      <  2:1032647/139‑1 (MQ=255)                                                               aGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                  <  2:66188/139‑1 (MQ=255)                                                                  cGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGt                                                               >  1:295167/1‑139 (MQ=255)                                                                    gagaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                             <  1:983663/139‑1 (MQ=255)                                                                          aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                       <  1:530473/139‑1 (MQ=255)                                                                          aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccccc                                                                                                                             <  2:279880/68‑1 (MQ=255)                                                                          aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccccc                                                                                                                             >  1:279880/1‑68 (MQ=255)                                                                                 tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                <  2:823645/139‑1 (MQ=255)                                                                                    gACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAg                                             <  2:818767/139‑1 (MQ=255)                                                                                                     aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaa                            >  2:663249/1‑139 (MQ=255)                                                                                                       cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                   <  1:593420/82‑1 (MQ=255)                                                                                                       cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                   >  2:593420/1‑82 (MQ=255)                                                                                                                  ttATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGttt               <  2:871808/139‑1 (MQ=255)                                                                                                                     tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcg                                                                                                >  1:806127/1‑55 (MQ=255)                                                                                                                     tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcg                                                                                                <  2:806127/55‑1 (MQ=255)                                                                                                                          aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGa                                                                                 <  2:135429/65‑1 (MQ=255)                                                                                                                          aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGa                                                                                 >  1:135429/1‑65 (MQ=255)                                                                                                                               aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGa         >  2:425886/1‑132 (MQ=255)                                                                                                                               aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGa         <  1:425886/132‑1 (MQ=255)                                                                                                                               aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag  >  1:888378/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                 AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑CCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600469‑2600732

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