Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,708,813	G→A	P555S (CCG→TCG)	eptB ←	phosphoethanolamine transferase EptB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,708,813	0	G	A	100.0%	217.0 / NA	62	P555S (CCG→TCG)	eptB	phosphoethanolamine transferase EptB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/21);  total (41/21)

CGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCT  >  NC_000913/3708519‑3709108                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        ggcctacgtcgagtcgaacctccctgcccccggaggttctcatcctcatgagctgcagatgcaaaaaagcctgctcgttgagcaggcttttcgaatttggtcggtgatagaggattCGCACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGTATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTACTCGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCACTTAGTTAGCCGCTGCATCTTTTGCCTACGAGAtgt                                                                                                                                                                                                                                                                                                    >  1:22068‑M2/117‑300 (MQ=255) ggcctacgtcgagtcgaacctccctgcccccggaggttctcatcctcatgagctgcagatgcaaaaaagcctgctcgttgagcaggcttttcgaatttggtcggtgatagaggattCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCCTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTGGCCGCTGCCTCTTTTGCCTGCGAGAtgt                                                                                                                                                                                                                                                                                                    >  2:4235‑M2/117‑300 (MQ=255) ggcctacgtcgagtcgaacctccctgcccccggaggttctcatcctcatgagctgcagatgcaaaaaagcctgctcgttgagcaggcttttcgaatttggtcggtgatagaggattCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCCTAATAGCAAAATGCCTTTTGATCGGCGAGAAGGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGAtgt                                                                                                                                                                                                                                                                                                    >  2:78051‑M2/117‑300 (MQ=255) ggcctacgtcgagtcgaacctccctgcccccggaggttctcatcctcatgagctgcagatgcaaaaaagcctgctcgttgagcaggcttttcgaatttggtcggtgatagaggattCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTNACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGAtgt                                                                                                                                                                                                                                                                                                    >  1:507733‑M2/117‑300 (MQ=255) ggcctaagtcgcgtcgaacctccctgcccccggaggttctcatcctcatgagctgcagatgtaaaaaagcctgctcgttgagtaggcttttcgaatttggtcggtgatagaggattCGAACCTCCGACCCCTTCGTCCAGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTAATGCGCAGATACGCCCTCGTCAATCCCTTAATAGCCAAATGTCTTTTGATCGGNGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCNAGAtgt                                                                                                                                                                                                                                                                                                    >  2:574674‑M2/117‑300 (MQ=255)                          gCCCCCGGAGGTTCTCATCCTCATGAGCTGCCGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTAGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAAGGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTNAGCCGCTGCCTCTTTTGCCTGCGAGGTGTGACACCAGTTGTTGTTTTNGTTAAt                                                                                                                                                                                                                                                                           >  2:576195/1‑300 (MQ=255)                          gCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAAt                                                                                                                                                                                                                                                                           >  1:520935/1‑300 (MQ=255)                            ccccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTc                                                                                                                                                                                                                                                                         >  1:514056/1‑300 (MQ=255)                                ggAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTccacc                                                                                                                                                                                                                                                                     <  1:367824/300‑1 (MQ=255)                                 gAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTccacca                                                                                                                                                                                                                                                                    >  1:49245/1‑300 (MQ=255)                                     ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGAAACCAGTTGTTGTTTTCGTTAATTCCACCATCCg                                                                                                                                                                                                                                                                >  1:261216/1‑300 (MQ=255)                                             cTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGta                                                                                                                                                                                                                                                        >  1:376890/1‑300 (MQ=255)                                                       gCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGAc                                                                                                                                                                                                                                              <  2:371480/300‑1 (MQ=255)                                                                  aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCatgatg                                                                                                                                                                                                                                   <  1:271877/300‑1 (MQ=255)                                                                  aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCatgatg                                                                                                                                                                                                                                   <  2:114608/300‑1 (MQ=255)                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGNGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTAc                                                                                                                                                                                                                          <  1:508835/300‑1 (MQ=255)                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTAc                                                                                                                                                                                                                          >  1:436246/1‑300 (MQ=255)                                                                              ctgcgagatgtgtataagagacagcGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGc                                                                                                                                                                                                                       >  1:697295/25‑300 (MQ=255)                                                                                     ggCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCTAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGt                                                                                                                                                                                                                <  2:668304/300‑1 (MQ=255)                                                                                     ggCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGt                                                                                                                                                                                                                <  1:341537/300‑1 (MQ=255)                                                                                     ggCTTTGCGAATTTGGTCGTTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCTGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCTGGAAAGTCAGCAGGCTTCTTAGTTAGCCGCTGCCTCTTTTGCCTGCTAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGTTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGt                                                                                                                                                                                                                <  2:238579/300‑1 (MQ=255)                                                                                        ttttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTGAGCCACTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGAc                                                                                                                                                                                                             <  2:42884/300‑1 (MQ=255)                                                                                          ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACgg                                                                                                                                                                                                           <  1:211122/300‑1 (MQ=255)                                                                                                   ggTCGGTGGTAGAGGATTCGAACCTCCGACCCCTTCGTACCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTCTGNCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACt                                                                                                                                                                                                  >  1:467205/1‑300 (MQ=255)                                                                                                   ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACt                                                                                                                                                                                                  >  1:110515/1‑300 (MQ=255)                                                                                                       ggTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCa                                                                                                                                                                                              <  1:384436/300‑1 (MQ=255)                                                                                                        gTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCAt                                                                                                                                                                                             >  1:221950/1‑300 (MQ=255)                                                                                                        gTGATAGAGGATTCAAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGAAACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCATACAGCTCTACGTGACGGCGTGGCACTTTCAt                                                                                                                                                                                             >  1:225225/1‑300 (MQ=255)                                                                                                                ggATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGctt                                                                                                                                                                                     >  1:452784/1‑300 (MQ=255)                                                                                                                           tCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGc                                                                                                                                                                          >  1:205052/1‑300 (MQ=255)                                                                                                                           tCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGc                                                                                                                                                                          <  2:660567/300‑1 (MQ=255)                                                                                                                                 cccTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCa                                                                                                                                                                    >  1:518495/1‑300 (MQ=255)                                                                                                                                       ggcccGAACGAAGTGCGCTACCAGGCTACGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGAGACGCCCTCGGCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGcc                                                                                                                                                              <  2:361955/298‑1 (MQ=255)                                                                                                                                       gTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGCAGTATAGCCAAGACAACCCATGATGGTATCGTAAAGCTCTACGTGACGACGTGGCACTTTCATGTCGGCTTATTTTTTCAGCTGCGCaaccgac                                                                                                                                                              >  2:366021/1‑295 (MQ=255)                                                                                                                                       gTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGcc                                                                                                                                                              >  1:112363/1‑300 (MQ=255)                                                                                                                                            gaacgaaGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTgcgc                                                                                                                                                         <  1:581402/300‑1 (MQ=255)                                                                                                                                              acgaaGTGCGCTACCAGGCTGCGCCAATCACCGAATGCTGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAGAATGCCTTTTGATCGGCTAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCAt                                                                                                                                                       <  2:379222/300‑1 (MQ=255)                                                                                                                                                        cTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGAGACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATNGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGAtt                                                                                                                                             >  2:603515/1‑300 (MQ=255)                                                                                                                                                        cTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTNGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCATGCGCATTGGCCGGAtt                                                                                                                                             >  2:150605/1‑300 (MQ=255)                                                                                                                                                        cTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGAtt                                                                                                                                             >  1:129781/1‑300 (MQ=255)                                                                                                                                                                  gcgcCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATCCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCATGCGCATTGGCCGGATTTTCCAGatat                                                                                                                                   >  2:10712/1‑300 (MQ=255)                                                                                                                                                                       aaTCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCACACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATc                                                                                                                              >  2:521514/1‑300 (MQ=255)                                                                                                                                                                           aCCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGAc                                                                                                                          >  2:322615/1‑300 (MQ=255)                                                                                                                                                                                          cATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCatcatc                                                                                                           >  1:197348/1‑300 (MQ=255)                                                                                                                                                                                              ttACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTa                                                                                                       >  1:551175/1‑300 (MQ=255)                                                                                                                                                                                                    tgcaaaTACGCCCTCTGCAATCCCTTAATAGCAAAATGCCTTTTGGTCGGCTAGAAAGTCAGTAGGCCGCTTAGTGAGTCGTTGCCTCTTTTGCCTGCGAGATGTTACACCAGTTGTTGTTTTCGTTAATTCCACCGTCCGGTGAAGTATAGCCAAGACAACCCATGATGGGATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTGCAGCTGTGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGa                                                                                                 <  2:335099/295‑1 (MQ=255)                                                                                                                                                                                                        gATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTg                                                                                             >  1:203505/1‑300 (MQ=255)                                                                                                                                                                                                           aCGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTc                                                                                          <  1:96039/300‑1 (MQ=255)                                                                                                                                                                                                                cTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGGGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCcggcg                                                                                     >  2:589155/1‑300 (MQ=255)                                                                                                                                                                                                                cTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCcggcg                                                                                     >  2:378959/1‑300 (MQ=255)                                                                                                                                                                                                                  cGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGt                                                                                   >  2:625292/1‑300 (MQ=255)                                                                                                                                                                                                                               aTAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCAcgcg                                                                      <  1:239905/300‑1 (MQ=255)                                                                                                                                                                                                                                               tGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTg                                                      <  1:151333/300‑1 (MQ=255)                                                                                                                                                                                                                                                            aGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCAGTGCAGGTGCTCGCGTTCATTa                                         >  2:23066/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           tagttagCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCANATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTgg                          >  1:506490/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           tagttagCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTgg                          >  2:252278/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           tagttagCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTgg                          >  1:520540/1‑300 (MQ=255)                                                                                                                                                                                                                                                                             gttagCCGCTGCCTCTTTTGCCTGCTAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTATGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTc                        <  2:111772/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             gttagCCGCTGCCTCTTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTc                        >  1:631585/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          tctTTTGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCAGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAac           >  2:428732/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                               tGCCTGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAAt      <  1:422337/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   tGCGAGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCt  >  1:604227/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        CGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATAGCCAAGACAACCCATGATGGTATCGTACAGCTCTACGTGACGGCGTGGCACTTTCATGTCGGCTTCTTTTTTCAGCTGCGCAAACGCCTGCGCATTGGCCGGATTTTCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCT  >  NC_000913/3708519‑3709108

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