Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,682,459	A→G	K101E (AAA→GAA)	rstA →	DNA‑binding transcriptional regulator RstA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,682,459	0	A	G	100.0%	223.4 / NA	62	K101E (AAA→GAA)	rstA	DNA‑binding transcriptional regulator RstA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/27);  total (35/27)

GTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTACCCATGCCGGGCAAATCATGGA  >  NC_000913/1682215‑1682703                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCAGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAacg                                                                                                                                                                                                                                                 >  2:516126/1‑250 (MQ=255)        gCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTAGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGAACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGCCGccccc                                                                                                                                                                                                                                         >  2:699729/1‑251 (MQ=255)         cAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCct                                                                                                                                                                                                                                        <  1:494247/251‑1 (MQ=255)          aaaaCATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCctg                                                                                                                                                                                                                                       >  2:499648/1‑251 (MQ=255)               aTGATATGCAGGTTACCGGAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGtt                                                                                                                                                                                                                                  <  2:213919/251‑1 (MQ=255)                  atatGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATGTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTg                                                                                                                                                                                                                               <  2:212338/251‑1 (MQ=255)                    atGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCt                                                                                                                                                                                                                             >  2:804195/1‑251 (MQ=255)                        aGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAgcg                                                                                                                                                                                                                         <  1:807441/251‑1 (MQ=255)                             aCCGTAGAGCCGCGCAGCGACCAGGCCGAATAAATCATTTTGCGAGAAAATCCGGATTTGGTGTTACGCGACATCATGCTACCAGGCACGGACGGCATGACCATTTGTCGTGATTTACGCGCAACGTGGTCTGGACCGATTGTTCCTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGACATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGtt                                                                                                                                                                                                                     <  1:481438/250‑1 (MQ=255)                              ccGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGttt                                                                                                                                                                                                                    >  1:157957/1‑250 (MQ=255)                               cGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGCCTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTAc                                                                                                                                                                                                                  >  2:741563/1‑251 (MQ=255)                                 tAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGGTTTGGGGTTACTCGACATCATGCTACCAGGCAAGGACGGCAGGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTTTAACCTCTCTCGATAGCGATATGAACCACGTCCTGGCAATGGAAATGAGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGAGTTTACGt                                                                                                                                                                                                                >  2:432930/1‑251 (MQ=255)                                    agCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGCCTATATTCTCGAAACGACGCACCCTGCTGTTTTGCTAGCGCGTTTACGTTTg                                                                                                                                                                                                             >  2:598577/1‑251 (MQ=255)                                       cgcgcgGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCAt                                                                                                                                                                                                          >  1:666082/1‑251 (MQ=255)                                                  cAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGGTTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGaa                                                                                                                                                                                               <  2:299905/251‑1 (MQ=255)                                                   aGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAAt                                                                                                                                                                                              >  2:198937/1‑251 (MQ=255)                                                         aagaaACCATTGTGCGAGAAAAACCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTTGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCaa                                                                                                                                                                                        <  2:10365/251‑1 (MQ=255)                                                               ccATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATATCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCcac                                                                                                                                                                                   >  1:179914/1‑250 (MQ=255)                                                                cATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCacac                                                                                                                                                                                 <  2:292117/251‑1 (MQ=255)                                                                 aTTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGTCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTCCGTCAGAATGAGCAAGCCCCACt                                                                                                                                                                                >  2:547644/1‑251 (MQ=255)                                                                 aTTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACt                                                                                                                                                                                >  1:582451/1‑251 (MQ=255)                                                                               cccGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTCCGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTc                                                                                                                                                                  <  2:163676/250‑1 (MQ=255)                                                                                ccGGATTTGGTGTTACTCGAAATCATGCTACCAGGCAAGGACGGCAGGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTGCTTCTAACCTCTCTCGATAGCGATAGGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCa                                                                                                                                                                 <  2:488814/251‑1 (MQ=255)                                                                                  ggATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAgg                                                                                                                                                               <  1:88654/251‑1 (MQ=255)                                                                                                 cGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAGGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTc                                                                                                                                                >  2:225177/1‑251 (MQ=255)                                                                                                                        aCGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGc                                                                                                                         <  1:678433/251‑1 (MQ=255)                                                                                                                          ggCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCAc                                                                                                                       >  2:547950/1‑251 (MQ=255)                                                                                                                          gacATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCAc                                                                                                                       >  2:603206/3‑251 (MQ=255)                                                                                                                             aTGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGtt                                                                                                                    >  1:529443/1‑251 (MQ=255)                                                                                                                               gACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCca                                                                                                                                                                                    >  1:238376/1‑185 (MQ=255)                                                                                                                               gACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCca                                                                                                                                                                                    <  2:238376/185‑1 (MQ=255)                                                                                                                                 ccATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGCTTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGAcc                                                                                                                >  1:82981/1‑251 (MQ=255)                                                                                                                                   aTTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCAt                                                                                                              >  1:32022/1‑251 (MQ=255)                                                                                                                                      tGTCGTGATTTACGCGTAAAGTGGTCGGGACCGTTTGTATTGTTAACCTCTCACGGTAGCGATATTAACCACATCCTTGCTCCGGAAATGGGTTCCTGCTACTGTATTCTCGAAACTACGCCCCTTGATGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCATAATGAGCAAGCTACACTGACCAAAGGTATTCAGGAAACTTCTCGGACTTCCTACAAAGCCCTGCATTTCGGCTCGTTGTCCaccg                                                                                                            <  2:546287/250‑4 (MQ=255)                                                                                                                                        tcgagagttACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGGTAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACAATATTCTAGAAACGACGCCCCCTGCTGTTTTGCGAGCGCGTTTACGTTTGTATTTGCGTCAGCATGAGCAAGCCACACTGACCAAAGGCCTTAAGGAACCGTCTCTCACTCCCTCCAAAGCCATGCATTTCGGCACGTTGACCagcggtc                                                                                                         >  2:636387/8‑245 (MQ=255)                                                                                                                                          gTGATTTACGCGCAAAGTGGTCGGGACCGATTGTTCTTCTAACCTCTCGCGATAGCGATACGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCGGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCGGCATTTCGGCACGTTGACCATCGATccc                                                                                                       <  2:14896/251‑1 (MQ=255)                                                                                                                                                tACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCaa                                                                                                  <  1:76942/250‑1 (MQ=255)                                                                                                                                                  cgcgCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACcg                                                                                               >  1:164287/1‑251 (MQ=255)                                                                                                                                                          gTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTc                                                                                        >  1:409058/1‑250 (MQ=255)                                                                                                                                                                 ggACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTg                                                                                  <  1:499648/249‑1 (MQ=255)                                                                                                                                                                   aCCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCt                                                                               >  1:155435/1‑250 (MQ=255)                                                                                                                                                                    ccGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGCCCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTa                                                                              >  2:33223/1‑250 (MQ=255)                                                                                                                                                                      gATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCAAATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAac                                                                            >  1:599918/1‑250 (MQ=255)                                                                                                                                                                            tcttctAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGaa                                                                      <  1:124018/250‑1 (MQ=255)                                                                                                                                                                                 tAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCg                                                                >  2:349790/1‑251 (MQ=255)                                                                                                                                                                                     ctctctCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGACTGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCACTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGctct                                                            >  2:396547/1‑251 (MQ=255)                                                                                                                                                                                       cgcgcgatatcgatatGTACCACATCCTGGCACTGGAATTGGGTGCCTGCGATTATATTCTCGAAACGACGCCCCCTGCTTTTTTGCTAGCGCTTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCTTTCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCg                                                          <  1:232270/247‑1 (MQ=255)                                                                                                                                                                                        tctcGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCg                                                          >  1:626797/1‑250 (MQ=255)                                                                                                                                                                                               aGCGATATGAACCACATCCTGGCACTGGCAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGATCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACCGCTg                                                   >  1:92632/1‑250 (MQ=255)                                                                                                                                                                                                     attaaCCCCATCCTGGCACTGGAAAGGGGTGCTGGCGACTATATTCTCGAAACGACGCCCCCGTCTTTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGa                                            <  2:719852/248‑1 (MQ=255)                                                                                                                                                                                                                ccTGGCACTGGAAAGGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATGATTGTggg                                 <  1:470909/251‑1 (MQ=255)                                                                                                                                                                                                                    gCACTGGAAATGGGTGCCTGCGACTATATTCCCGAAACGACGCCCCCGTCTGTTTTGCTAGCGCGTTTACGTTTGCATTGGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTAGTGTGGGAAtt                             <  2:773489/251‑1 (MQ=255)                                                                                                                                                                                                                     cACTGGAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAAtt                             >  1:86150/1‑250 (MQ=255)                                                                                                                                                                                                                         ggAAATGGGTGCCTGCGATTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGc                          <  2:157957/249‑1 (MQ=255)                                                                                                                                                                                                                         ggAAATGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCTTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCTCGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTTACTCTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCt                         <  1:516126/250‑1 (MQ=255)                                                                                                                                                                                                                         ggAAATGGGGGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGGATTAGCt                         <  2:280268/250‑1 (MQ=255)                                                                                                                                                                                                                            aaTGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACACAGCCCTGCATTTCGGCACGTTGCCCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTAccc                     >  2:212104/1‑251 (MQ=255)                                                                                                                                                                                                                             atggggtcccgCGACTATATTTTCGAAACGGCGCCCCGTGGTGTTTTGCTAGCGGGTTTACGTTTGCATTGGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTACCCa                    <  2:179914/241‑1 (MQ=255)                                                                                                                                                                                                                             aTGGGTGCCTGCGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTAccc                     >  1:197163/1‑250 (MQ=255)                                                                                                                                                                                                                                        cGACTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTCCCCATGCCGGGCaa          >  2:397676/1‑250 (MQ=255)                                                                                                                                                                                                                                           cTATATTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGCATTATTGTGGGAATTAGCTACCCATGCCGGGCAAATc       >  1:231504/1‑250 (MQ=255)                                                                                                                                                                                                                                               atTCTCGAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGAATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTACCCATGCCGGGCAAATCATGGa  <  1:603206/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGATCCCATCAACCGCGTAGTCACCCTGGCTAACACTGAAATCTCGCTCTCGACAGCTGATTTCGAATTATTGTGGGAATTAGCTACCCATGCCGGGCAAATCATGGA  >  NC_000913/1682215‑1682703

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