Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,433	C→A	A319D (GCC→GAC)	msrP →	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,433	0	C	A	100.0%	129.2 / NA	43	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/22);  total (21/22)

GCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCT  >  NC_000913/2040188‑2040657                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  gCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAgg                                                                                                                                                                                                                              >  2:3467/1‑251 (MQ=255)           gtgCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCt                                                                                                                                                                                                                    <  1:153253/251‑1 (MQ=255)              ccGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGt                                                                                                                                                                                                                  >  2:434182/1‑250 (MQ=255)                           gCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCTCCTGGAAGCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGAt                                                                                                                                                                                                   <  2:96619/251‑1 (MQ=255)                                aaacggatttaaaagATAGTCAGTATTATGCTGACCCGCGAGCGTCCGCCAACCACCTTGAATCTGGCTGCGCATGAGGAATACGGTTTTTACGCCCACGTTAATCCGTATGTTGATCATCCGCGCTGGTCACAGGCTACCGAACGAGTTATTGGTTCAGTCGGCATCCTCGATGTACAGCGCCTGACCACGCTAATGTGGAATGGTTACGCCGACCAGGTGGCATCGC‑GTGATCGTGGCCTGGATTTGc                                                                                                                                                                                               <  2:90471/236‑1 (MQ=255)                                aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTg                                                                                                                                                                                                >  1:135799/1‑249 (MQ=255)                                    ggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGgag                                                                                                                                                                                          <  1:281369/251‑1 (MQ=255)                                             cGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTaa                                                                                                                                                                                 <  1:437065/251‑1 (MQ=255)                                                 aGTCAGTATTAAGCTGACCCGAGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGc                                                                                                                                                                             >  2:23552/1‑251 (MQ=255)                                                     aGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGCTTTGCGGGAGAATTTCTAAATGCGTCt                                                                                                                                                                         >  2:417425/1‑251 (MQ=255)                                                          tAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                     <  1:174321/250‑1 (MQ=255)                                                                  cccGCGAGCGTCCGCCAACCCCCTGGAATCTGGCAGCGCCTGACGAATACGGTTGTTACGCCAACGGTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGTAGAATTTCTAAATGCGTCTGACAGCAAAACAg                                                                                                                                                            <  2:288054/251‑1 (MQ=255)                                                                         gCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACAt                                                                                                                                                     <  2:361905/251‑1 (MQ=255)                                                                               gCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACCGTTTTTACGCCAACCTTAATCCGTATGCTGATCACCCGCGATGGTCACCGGCAACAGAACGATTTATTGGATAAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTg                                                                                                                                                >  1:213044/1‑250 (MQ=255)                                                                                 cAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGCTTTTACGCCAACGCTCAACCGAATGTTGATCAAACGCGCTGGTCACAGGCTACCGAACGATTTATTGGTACAGGCGGCATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGaa                                                                                                                                              >  1:188817/1‑250 (MQ=255)                                                                                       ccTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGc                                                                                                                                       <  1:349888/251‑1 (MQ=255)                                                                                                 ggcagagacTGACGAATACGTTTTTTACGCCAACGTTAAGCCGTTTGTTGATCACCCTCGTTGGTTACTGGCTACCGAACGATTTATTGTTTCAGGCGGCATCCTCGTGGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGG‑ATCGTTGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGTTGACATGGCTGCAAGTTTTCCTGCAt                                                                                                                                 <  2:240304/239‑1 (MQ=255)                                                                                                           gACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttg                                                                                                                     <  1:9542/249‑1 (MQ=255)                                                                                                                aTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACGGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGt                                                                                                              <  2:386260/251‑1 (MQ=255)                                                                                                                     gTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGCTGGTCCAAGGCTACCGAACGATTTCTTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTCCGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTtctc                                                                                                         >  1:55438/1‑251 (MQ=255)                                                                                                                       ttttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTTGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCCGCAGCTTGCCGTATTGTTGCCGTTTctct                                                                                                        <  2:228831/250‑1 (MQ=255)                                                                                                                            cGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGCTTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTctggctgg                                                                                                  >  2:373644/1‑251 (MQ=255)                                                                                                                            cGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTTGTTCAGGAGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCCGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGTGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTTAACGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTctggctgg                                                                                                  >  1:185292/1‑251 (MQ=255)                                                                                                                                   gTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGc                                                                                           <  2:25945/251‑1 (MQ=255)                                                                                                                                   gTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGc                                                                                           >  2:445241/1‑251 (MQ=255)                                                                                                                                   gTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGc                                                                                           <  1:417425/251‑1 (MQ=255)                                                                                                                                    ttAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCg                                                                                          >  1:395544/1‑251 (MQ=255)                                                                                                                                         ccGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCGAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGatca                                                                                      >  1:261983/1‑250 (MQ=255)                                                                                                                                             aTGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGGTGCCGTTTCTCTGGCTGGTCTGGGCGatcaatca                                                                                  <  1:65208/250‑1 (MQ=255)                                                                                                                                                            tcttttCACAGGCTACCGAACGCTTTATTGGTTCAGGCGGCATCCTCGCTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGAAATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGCTCAATCACGGTGGACTGGGTGcc                                                                  >  2:249749/6‑251 (MQ=255)                                                                                                                                                               ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGTTGCCGAt                                                               <  2:195777/251‑1 (MQ=255)                                                                                                                                                                     aGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTg                                                         <  1:39906/251‑1 (MQ=255)                                                                                                                                                                             gAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGatat                                                 >  1:246010/1‑251 (MQ=255)                                                                                                                                                                                  aTTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGCTCCGGTGAAAGATATTCAGc                                            >  2:71752/1‑251 (MQ=255)                                                                                                                                                                                  aTTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                            >  1:240272/1‑251 (MQ=255)                                                                                                                                                                                     tATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATACTCAGCAt                                          <  2:90855/250‑1 (MQ=255)                                                                                                                                                                                        tGGTTCCGGCGGCATCCTCGATGTACAGCGCACGCCCACGCTACTGTTTAATGGTCACGTCGACCAGGAGCCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGCGACATGGCTGACAGTTTGCATGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCAtttt                                       >  1:170208/1‑250 (MQ=255)                                                                                                                                                                                               ggcggcATCCTCGATGTACAGCGCCAGCCAACGCTGCTGTTTAATGGTTACGCCGACCAGGTTGCAGCGCTGT‑ATCGGGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCGGAATGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGGTGGTCTGGGCGATCAATCACGGTGTACTCGGTGCCTATCCGGTGAAAGATATTCAGCATTTTACTGGTCg                               <  1:306270/251‑1 (MQ=255)                                                                                                                                                                                                       ccTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcac                       >  1:423951/1‑251 (MQ=255)                                                                                                                                                                                                        cTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttgtt                                                                                                                   >  1:327942/1‑158 (MQ=255)                                                                                                                                                                                                        cTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGAGTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttgtt                                                                                                                   <  2:327942/158‑1 (MQ=255)                                                                                                                                                                                                              gTACAGCGCCAGCCAACGCTACTGTTTAATGGTCACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGCTCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAtt                >  1:67980/1‑251 (MQ=255)                                                                                                                                                                                                                            aaCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAGCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCTACCt  <  2:135799/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  GCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGT‑ATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCT  >  NC_000913/2040188‑2040657

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