Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,566,734	G→A	P51S (CCG→TCG)	eutA ←	ethanolamine ammonia‑lyase reactivase EutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,566,734	0	G	A	100.0%	217.4 / NA	60	P51S (CCG→TCG)	eutA	ethanolamine ammonia‑lyase reactivase EutA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/29);  total (31/29)

GTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGGACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAA  >  NC_000913/2566440‑2567027                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      gTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACttt                                                                                                                                                                                                                                                                                                  >  1:547356/1‑300 (MQ=255)      cGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGcaag                                                                                                                                                                                                                                                                                             >  2:279031/1‑297 (MQ=255)             gCCGTTGTTCAGAAATGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAAttt                                                                                                                                                                                                                                                                                      <  2:125479/300‑1 (MQ=255)                            gggTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGATAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAAtt                                                                                                                                                                                                                                                                       <  2:480306/300‑1 (MQ=255)                            gggTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAAtt                                                                                                                                                                                                                                                                       >  1:438591/1‑300 (MQ=255)                             ggTTTGTGCCCAGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAACGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCAGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTc                                                                                                                                                                                                                                                                      >  2:169797/1‑300 (MQ=255)                                         ggCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCAc                                                                                                                                                                                                                                                          >  1:388135/1‑300 (MQ=255)                                                  aCCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGTTGATGATGATGGCACCAGAATCAACGCTTTCCNGCTCAATACCCGCAGCATGATATTGCNCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGacac                                                                                                                                                                                                                                                 <  2:561882/300‑1 (MQ=255)                                                      gCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTNTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGANAccgcc                                                                                                                                                                                                                                             >  2:707584/1‑300 (MQ=255)                                                          tCACGGATTCTAGGGGCGGCCCTGCGCTGGCAGCTACAAAATCGCCCAGCGACTGAGATAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTAGTGAGGATGGCACCAGAATCAACGCTTTCCTGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTCCCGCTTCTTTTAAACCACCCTGTTTATCGACAGAGGTAAAGAACACCGAACTTTGCCAACTAACTTCGCGCTTAATGAATTCTTAGCCCGGCACCTGCGACACCCCCGCAc                                                                                                                                                                                                                                         <  2:255098/300‑1 (MQ=255)                                                          tCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTNAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTNCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCacac                                                                                                                                                                                                                                         <  2:797204/300‑1 (MQ=255)                                                            aCGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGTTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACgg                                                                                                                                                                                                                                       <  1:421129/300‑1 (MQ=255)                                                                       gTGCGGCCCGGCGCTGGCAACGACAAAATCGACCAGCGAATGAGAGAGCGCCATCACAGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCGGAATCAACGCTTTCCGGCTCACTACCCGCCACATGATATTGATCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTCATTTCGCATTTAATGAATCCATAGCTCGGCACCTGTGACACCGACGCGGGGTTAACCAGCTc                                                                                                                                                                                                                            >  2:36618/1‑300 (MQ=255)                                                                            gCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGcc                                                                                                                                                                                                                       <  1:533387/300‑1 (MQ=255)                                                                            gCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTATTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGcc                                                                                                                                                                                                                       >  2:239342/1‑300 (MQ=255)                                                                                 gcgcTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGgag                                                                                                                                                                                                                  >  1:518846/1‑300 (MQ=255)                                                                                                 aaTCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCTGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAACCACCTGggtgg                                                                                                                                                                                                  <  2:437052/300‑1 (MQ=255)                                                                                                    cGCCCAGCGACTGAGATAGCGCCATCACCGCCGGGCGAGCATTACGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCTGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCTTAGCGCTGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGTCGGGAGAAAATCACCTGggtggtgg                                                                                                                                                                                               <  2:274956/300‑1 (MQ=255)                                                                                                    cGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGggtggtgg                                                                                                                                                                                               <  1:823890/300‑1 (MQ=255)                                                                                                     gCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCTGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGggtggtggt                                                                                                                                                                                              <  1:354687/300‑1 (MQ=255)                                                                                                      cccAGCGACTGAGATAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGTTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTGCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGggtggtggtg                                                                                                                                                                                             <  2:170640/300‑1 (MQ=255)                                                                                                           cGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGcc                                                                                                                                                                                        <  1:106145/300‑1 (MQ=255)                                                                                                            gACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCg                                                                                                                                                                                       >  1:302119/1‑300 (MQ=255)                                                                                                                gagagaGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGatat                                                                                                                                                                                   <  1:183004/300‑1 (MQ=255)                                                                                                                  gagagCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCCGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCg                                                                                                                                                                                 >  2:123382/1‑300 (MQ=255)                                                                                                                       cgcCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGAAACCGCCGCACGATTAACCAGCTCCAGCCAAGAGAAAATCAAATGGGTGGTGGTGGTGCCGATATCGCTAcc                                                                                                                                                                            >  2:353536/1‑300 (MQ=255)                                                                                                                       cgcCATCACCGCCGGGCGAGAATTGTGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCTGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGTGGTAAAGAACACCGAACTTTGCCAGCGAATTTCGCGTTTAATGAATGCGTAGCGCGGTACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCATCGGGGTGGTGGTGGTGCCGATATCGATAcc                                                                                                                                                                            <  1:313423/300‑1 (MQ=255)                                                                                                                                    gggCGAGCATTGTGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTAGTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGc                                                                                                                                                               <  2:109964/300‑1 (MQ=255)                                                                                                                                         aGCATTGCGGGTTTTCGCGCTTTCACCGTTGATGAGGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGTTTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCg                                                                                                                                                          <  2:512001/300‑1 (MQ=255)                                                                                                                                           cATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGNGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAg                                                                                                                                                        >  2:761229/1‑300 (MQ=255)                                                                                                                                                     tttCGCGCTTTCACCTGTGATGATGATGGCACCAGAATCAACGCTTTCCTGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCATTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAATTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGTAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGAt                                                                                                                                              <  2:353561/300‑1 (MQ=255)                                                                                                                                                     tttCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGAt                                                                                                                                              <  1:361948/300‑1 (MQ=255)                                                                                                                                                                          gatgatgGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTa                                                                                                                         >  1:94743/1‑300 (MQ=255)                                                                                                                                                                                        aaTCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGtc                                                                                                           >  1:14942/1‑300 (MQ=255)                                                                                                                                                                                           cAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGtcttc                                                                                                        <  1:437150/300‑1 (MQ=255)                                                                                                                                                                                                  ttCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAANTCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGc                                                                                                 >  1:605152/1‑300 (MQ=255)                                                                                                                                                                                                              aCCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCatcatc                                                                                     >  1:401024/1‑300 (MQ=255)                                                                                                                                                                                                                     gcaTGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCCCGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATcgccaagc                                                                              >  2:433601/1‑297 (MQ=255)                                                                                                                                                                                                                     gcaTGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATcgccacgc                                                                              >  2:805782/1‑300 (MQ=255)                                                                                                                                                                                                                            atTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGc                                                                       <  2:836688/300‑1 (MQ=255)                                                                                                                                                                                                                                  cGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTAc                                                                 <  2:508977/300‑1 (MQ=255)                                                                                                                                                                                                                                     gTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGCTACCGACGATCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTAcgcc                                                              >  2:458016/1‑300 (MQ=255)                                                                                                                                                                                                                                     gTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTAcgcc                                                              >  2:16892/1‑300 (MQ=255)                                                                                                                                                                                                                                         tAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAANAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACGAGCTCCAGCCGGTAGAAAATCACCTGGGTGGTGTTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGccgccg                                                          <  2:632868/300‑1 (MQ=255)                                                                                                                                                                                                                                               tttCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAAc                                                    <  1:557286/300‑1 (MQ=255)                                                                                                                                                                                                                                                          cttcttTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGTAGAAAATCACCTTGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGACTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAat                                         <  2:500081/300‑1 (MQ=255)                                                                                                                                                                                                                                                          cttcttTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAat                                         <  2:500227/300‑1 (MQ=255)                                                                                                                                                                                                                                                             cttTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAatcat                                      >  1:700573/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            gTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGc                       >  2:813850/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            gTTTATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGc                       <  1:335145/300‑1 (MQ=255)                                                                                                                                                                                                                                                                            gTTTATCGACAGGGGTAAAGAACACCGAACTTTACCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTACGCTTAGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCTTTACTCCGCCTCTCAACTTGCCGACACTCATCGGGAAGATCATGGc                       >  2:1303/1‑300 (MQ=255)                                                                                                                                                                                                                                                                             tttATCGACAGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCg                      <  1:385214/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                      aGGGGTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCNCCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTg             <  1:598780/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        gggTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGGAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCAGACAATCATCGGGAAGATCATGGCGTTCATGTTGGc           >  2:422355/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        gggTAAAGAACACCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGc           >  1:300885/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        gggTAAAGAACACCGAACTTTGCCAGCTAATNTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTNCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGc           >  2:614771/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 acacCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTACCCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGaa  >  1:361901/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 acacCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGaa  >  2:733992/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 acacCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGaa  >  1:361914/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 acacCGAACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTACGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGaa  >  1:221068/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      GTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTTTCCGGCTCAATACCCGCAGCATGATATTGCTCGAGTATTAAGGTTTTCAGTTCCGCTTCTTTTAAACCGCCCTGTTTATCGACAGGGGTAAAGAACACCGGACTTTGCCAGCTAATTTCGCGTTTAATGAATTCGTAGCGCGGCACCTGCGACACCGCCGCACGGTTAACCAGCTCCAGCCGGGAGAAAATCACCTGGGTGGTGGTGGTGCCGATATCGATACCGACGCTCAATAGCTGGCGAGTGTTCACGATTGTGCCTCCGCTTCGGTTTTAGTCGCGGTCGCGTCTTCTTTTGGCACCAGCATCATCGCCACGCCAATCGCCGTTACGCCGCCGATCAACTTGCCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAA  >  NC_000913/2566440‑2567027

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