Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,638,571	G→A	H33Y (CAT→TAT)	cspI ←	cold shock‑like protein CspI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,638,571	0	G	A	100.0%	226.6 / NA	66	H33Y (CAT→TAT)	cspI	cold shock‑like protein CspI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/32);  total (34/32)

TTGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAC‑TGCTTCAGATAAATTTGTAACAAACC  >  NC_000913/1638327‑1638810                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 caataccctggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAANCGTTGCTCTGAATTGCTGAGAAATAGACaa                                                                                                                                                                                                                                             >  2:330349‑M2/141‑250 (MQ=255)          ggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACCTCttt                                                                                                                                                                                                                                   >  2:541471‑M2/132‑251 (MQ=255)             aaggatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGGTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACTCATCTTtgct                                                                                                                                                                                                                                <  2:121684‑M2/123‑1 (MQ=255)               tgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctgc                                                                                                                                                                                                                              >  2:444243‑M2/127‑251 (MQ=255)                    ccggttacccaaagctaggccggactgaatgnccaccgtcagtaatttttatacccggcgtaactgccgggatattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATTAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTATGCTGCCATCt                                                                                                                                                                                                                         <  2:353976‑M2/130‑1 (MQ=255)                      ggttacccaaagatatgacggactgaatgaccaccgtcagtaatttttatacccggcgtaactaccgggttatagcttgtcacaaaaaagtggaagactcatgtagttaacccactcacTACATGAACAGCGGAAGGCCCTTTAGGTCAGTTCCCAATACCAAATTCAACTTCCTGATTCTCCGTTTATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAATTAGACAAAAACATCTTTGCTGCCATCttt                                                                                                                                                                                                                       >  2:687541‑M2/120‑251 (MQ=255)                       gttacccaaagctatggcggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTc                                                                                                                                                                                                                      <  2:114613‑M2/133‑1 (MQ=255)                       ggtacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCCCCATTAATGTTTTGAAAACGTTGCTCTGAATTGCTGAGAAATAGACAAACACAACTTTGCTGCCATCTTTc                                                                                                                                                                                                                      >  2:136393‑M2/119‑251 (MQ=255)                         tacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCg                                                                                                                                                                                                                     <  2:538284‑M2/134‑1 (MQ=255)                             caaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                                                                                                                                                                <  1:145485‑M2/139‑1 (MQ=255)                                agctaggccggaatgaatgtccaccttcagtaaattttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATTAACAGCGGCAGGACCTTTAGGGCCGTTTTCAATACCAAATTCAGCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTTAGAAATATATAAATACATCTTTGCTGCCAGCTTTCGGCatgatg                                                                                                                                                                                                             <  2:673109‑M2/142‑1 (MQ=255)                                  ctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGGCCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                                                                                                                                                            >  2:249255‑M2/108‑250 (MQ=255)                                      gccgggctgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                                                                                                                                                                        <  2:527865‑M2/147‑1 (MQ=255)                                         ggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaaa                                                                                                                                                                                                    <  1:205259‑M2/151‑1 (MQ=255)                                           actgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                                                                                                                                                                  >  2:707512‑M2/99‑251 (MQ=255)                                              gaatgtccgccaccagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCAATCTTTCGGCGTGATGAAACCAAAACCttt                                                                                                                                                                                               <  2:620971‑M2/156‑1 (MQ=255)                                              gaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                                                                                                                                                               <  2:190992‑M2/156‑1 (MQ=255)                                               aatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtttt                                                                                                                                                                                              <  1:202861‑M2/157‑1 (MQ=255)                                                         gtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGTTGAAACCAAAACCTTTTTCTGGGTTa                                                                                                                                                                                     >  2:496934‑M2/85‑250 (MQ=255)                                                           cagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                                                                                                                                                   <  1:485286‑M2/168‑1 (MQ=255)                                                            agtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaaggggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTTTAAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGATAAATAGACAAACTCATCTTTGCTGCCATCTTTCGGCGTGATTAAACCAAAACCTTTTTCAGGGTTAAAcc                                                                                                                                                                                 <  2:590228‑M2/170‑1 (MQ=255)                                                             gtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTTCTCTGAATTGCTGAGAAATAGACAACCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAACCCTTTTTCAGGGTTAAaccc                                                                                                                                                                                >  2:661034‑M2/81‑250 (MQ=255)                                                               aatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                                                                                                                                              >  1:350983‑M2/79‑251 (MQ=255)                                                                    ttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                                                                                                                                                          <  1:492313‑M2/177‑1 (MQ=255)                                                                    ttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                                                                                                                                                          <  2:128487‑M2/177‑1 (MQ=255)                                                                          ccggcgtaactgccgggttattgctggtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAAGACCAAATTCAACGTCCTGATGCTCAGTTAATGTTTTGAAATCGTTGCTCTGTATTGCTGAGAAATAGACAAACACATCGTTGCTGCCAGCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAGCCACGTCACTAAACCa                                                                                                                                                                   <  2:236762‑M2/184‑1 (MQ=255)                                                                            ggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGTGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAACAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                                                                                                                                                                 >  2:686669‑M2/66‑251 (MQ=255)                                                                            ggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCAGCTTTCGGCGtgat                                                                                                                                                                                                              <  2:84410‑M2/141‑1 (MQ=255)                                                                            ggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGCTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                                                                                                                                                                              >  1:84410‑M2/66‑206 (MQ=255)                                                                              cgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCa                                                                                                                                                               <  1:262846‑M2/188‑1 (MQ=255)                                                                                        gggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAACAGTCATTTTGtt                                                                                                                                                        >  1:521631‑M2/54‑248 (MQ=255)                                                                                                    gtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACTAGTCATTTTGTTAGACATAATTATTAc                                                                                                                                         <  1:640878‑M2/210‑1 (MQ=255)                                                                                                            aaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTACCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTTCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaa                                                                                                                                 <  2:413249‑M2/218‑1 (MQ=255)                                                                                                                   tagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt                                                                                                                           >  1:256429‑M2/27‑250 (MQ=255)                                                                                                                    agactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCGTTTGAAGAAATTAg                                                                                                                         <  2:88007‑M2/226‑1 (MQ=255)                                                                                                                     gactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                                                                                                        <  1:218998‑M2/227‑1 (MQ=255)                                                                                                                          ctgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGATCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGGCATTTTGTTAGACATAATTCTTACCTTTTGCAGAAATTAGCCCTTg                                                                                                                   >  2:232642‑M2/20‑251 (MQ=255)                                                                                                                               gttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCt                                                                                                                                                                                                                         >  2:273984‑M2/15‑144 (MQ=255)                                                                                                                               gttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCt                                                                                                                                                                                                                         <  1:273984‑M2/130‑1 (MQ=255)                                                                                                                                     ccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGt                                                                                                        <  1:444243‑M2/243‑1 (MQ=255)                                                                                                                                       actcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa                                                                                                                          <  1:474661‑M2/225‑1 (MQ=255)                                                                                                                                       actcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa                                                                                                                          >  2:474661‑M2/7‑231 (MQ=255)                                                                                                                                           aCTACATGAACAGCGGCAGGACCTTTAAGTCCGTTCTCAATAACAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATATTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaa                                                                                                   >  1:284966/1‑250 (MQ=255)                                                                                                                                              aCATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaaa                                                                                                 <  1:323068/249‑1 (MQ=255)                                                                                                                                                  gAACAGCGGCAGGGCCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGCAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAat                                                                                            >  2:517874/1‑250 (MQ=255)                                                                                                                                                   aaCAGCGGCAGGGCCTTTAGGTCGGTTCTCCATACCACATTCAACATCCTGATTCTCCACTAATGTTTTGAAATCGTTGTTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCACTCTTTCGGCGTGATGAAATCAAAAACTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAATAAATTAGCCCTTGGGCAGAATGGTCCGaaaaaaaaatt                                                                                          >  2:631461/1‑248 (MQ=255)                                                                                                                                                     cAGCGGCAGGACCTTTAGAACCGTTCTCAATCCCAAATTCAACTTCCTACTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGGTGAAACCAAAACCTTTTTCATGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGCGCAGAATGGTCCGAAAAAAAATATc                                                                                         >  2:612322/1‑250 (MQ=255)                                                                                                                                                                         ccGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAGTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCGTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAgg                                                                     <  2:162231/250‑1 (MQ=255)                                                                                                                                                                           gTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGaa                                                                   >  1:128883/1‑250 (MQ=255)                                                                                                                                                                                   gacaaaaTTGAAGTTCCTGTTTCTCAGTTAATGTATTGCAATAGCTGCTCTATATTGCTGAGATATAGACAACCACATCTTTGCTGCCATCTTTCGACCTGACGAAACCAAACACTTTTTCAGGGTTAACCCCTTTCACTAAACCAGTCATTTTGTTAGACATACTTATTACCTCTTGAAGAAATTAGCCCTTTGGAAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAAACAACAAGGAAAtaacacg                                                           >  2:77712/5‑244 (MQ=255)                                                                                                                                                                                    aCCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAagag                                                         <  1:301270/251‑1 (MQ=255)                                                                                                                                                                                     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAgg                                                        >  1:325279/1‑251 (MQ=255)                                                                                                                                                                                            cAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                                                 >  2:98897/1‑251 (MQ=255)                                                                                                                                                                                             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaata                                                <  2:256429/251‑1 (MQ=255)                                                                                                                                                                                                  ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAAtaataa                                            >  1:128610/1‑250 (MQ=255)                                                                                                                                                                                                           cAGTTAATGTTTTGAGATCGTTGCTCTAAATTGATGAGAAAAAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTCGACCTACTTATTACCTTTTGAAGAAATTAGCCCTTAGGCAGAATAGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTAAAGAGGAACAACAAATAAATTTATAAc                                   >  2:664433/1‑250 (MQ=255)                                                                                                                                                                                                           cAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGCCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGATAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAAc                                   <  1:548983/250‑1 (MQ=255)                                                                                                                                                                                                           cAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATaa                                    <  1:350667/249‑1 (MQ=255)                                                                                                                                                                                                           cAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAAc                                   >  2:418584/1‑250 (MQ=255)                                                                                                                                                                                                               tAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACactg                               >  2:232032/1‑247 (MQ=255)                                                                                                                                                                                                               tAATGTTTTGAAATCGTTGCACTGAATTGCTGAGAAATAGACAAACACAGCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCATGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAAGAATAAAATTATAACAATg                               <  2:284966/250‑1 (MQ=255)                                                                                                                                                                                                                 aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCCTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAC‑t                           >  2:666237/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACAAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCACTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGATAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATAATAACAAT‑ACATGCTTCAGATAAATTTGTAACa      >  2:312215/1‑250 (MQ=255)                                                                                                                                                                                                                                         ttGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAC‑TGCTTCAGATAAATTTGTAACAAAc   >  1:451576/1‑251 (MQ=255)                                                                                                                                                                                                                                         ttGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGGACAAATAATCAAATTATAACAATGAC‑TGCTTCAGATAAATTTGTAACAAAc   >  2:155077/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAC‑TGCTTCAGATAAATTTGTAACAAAcc  >  1:321770/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 TTGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAC‑TGCTTCAGATAAATTTGTAACAAACC  >  NC_000913/1638327‑1638810

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