Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,317,944	C→T	intergenic (‑49/+62)	phnL ← / ← phnK	methylphosphonate degradation complex subunit PhnL/carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,317,944	0	C	T	100.0%	214.7 / NA	61	intergenic (‑49/+62)	phnL/phnK	methylphosphonate degradation complex subunit PhnL/carbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/33);  total (28/33)

TGCGCGCTGGCGCGGTGACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCC  >  NC_000913/4317649‑4318237                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      tGCGCGTTGGCGCGGTGACCAGGTCTACCCACTCGTCACCGGGTTTGATCTGGATTGGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTGGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGTGCAGGCGTACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGAc                                                                                                                                                                                                                                                                                                   <  2:303229/300‑1 (MQ=255) tGCGCGCTGGCGCGGTGACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGAc                                                                                                                                                                                                                                                                                                   <  1:33932/300‑1 (MQ=255)  gcgcgcTGGCGCGGTGACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCAGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGAcc                                                                                                                                                                                                                                                                                                  <  1:812776/300‑1 (MQ=255)        tGGCGCGGTGACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACggg                                                                                                                                                                                                                                                                                            <  1:685592/300‑1 (MQ=255)                 gACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCAtt                                                                                                                                                                                                                                                                                   >  2:536805/1‑300 (MQ=255)                 gACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCAtt                                                                                                                                                                                                                                                                                   >  1:489960/1‑300 (MQ=255)                 gACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGNTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGANCGCAtt                                                                                                                                                                                                                                                                                   >  2:673285/1‑300 (MQ=255)                      ggTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAg                                                                                                                                                                                                                                                                              >  1:680147/1‑300 (MQ=255)                         cTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCATAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGcc                                                                                                                                                                                                                                                                           <  1:118967/300‑1 (MQ=255)                               aCTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACGGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATaa                                                                                                                                                                                                                                                                     >  2:520623/1‑300 (MQ=255)                                cTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGCCGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATagg                                                                                                                                                                                                                                                                    >  2:257397/1‑298 (MQ=255)                                  cggcACCGTGTTTGATCTGGATTTGACCTTCGTCAGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCGCCACGCATTCGCCCACGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGATCTACGGGGAGCGCATTTGTAGGCCGGATAAgcc                                                                                                                                                                                                                                                                  <  2:418690/297‑3 (MQ=255)                                   gTCACCGTGTTTGATCTGGATTTGACCTTCGTCGTGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCNCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCNGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCg                                                                                                                                                                                                                                                                 <  2:647581/300‑1 (MQ=255)                                                aTCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                                                                                                                    <  1:265798/300‑1 (MQ=255)                                                    ggATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGNTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGGGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                                                                                                                                >  2:182852/1‑300 (MQ=255)                                                      aTTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGATCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAc                                                                                                                                                                                                                                              <  1:137068/300‑1 (MQ=255)                                                                 gtcgGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTc                                                                                                                                                                                                                                   <  1:257873/300‑1 (MQ=255)                                                                     ggCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCa                                                                                                                                                                                                                               >  1:750078/1‑300 (MQ=255)                                                                         gctcgTTGGCGTACAGCGAGCGTCGCAGAGTTGATTTGCCGCTGCCGGAATGGCAGTAGAGCACCACGCATTCGCCCGCGTTGCCGGTGAGCGAGGCGCGATTGAGGACGGGCAGTCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTATTATTCGACCTACGGGGCGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAAc                                                                                                                                                                                                                           >  2:362191/5‑300 (MQ=255)                                                                             gTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCACTTCTGCAAAACCGAt                                                                                                                                                                                                                       >  2:43450/1‑300 (MQ=255)                                                                                ggCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGACCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGac                                                                                                                                                                                                                    >  2:391720/1‑300 (MQ=255)                                                                                         cgagcgTAGCAGATTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCGCGCATTCGCCCGCGGTGGCGGTGAGCTAGGCGCGATTGAGGACGGGCAGGCGTACGCCGTGTTGCTGGTGCAGGATGAAGGTTTTATTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTAGGGTGTCGGATGCAACGCTGACGCGTCTTATTCGACCTACGGGGAGCGTATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTTCAAAACCGATGACACCagcagc                                                                                                                                                                                                           <  2:52084/300‑1 (MQ=255)                                                                                          gagcgTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCt                                                                                                                                                                                                          >  1:314927/1‑300 (MQ=255)                                                                                            gcgTAGCAGAGTTGATTTGCCGCTGCCGGAATGACCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGCTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTAGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCCAAACCGATGACACCAGCAGCtgt                                                                                                                                                                                                        >  2:738245/1‑300 (MQ=255)                                                                                               tcgCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGTGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGCGta                                                                                                                                                                                                     >  2:150051/3‑300 (MQ=255)                                                                                                aGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGtat                                                                                                                                                                                                    <  1:166306/300‑1 (MQ=255)                                                                                                   agagTTGATTTGCCGCTGCCGGAATGGCCTTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACg                                                                                                                                                                                                 <  2:336298/300‑1 (MQ=255)                                                                                                   agagTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACg                                                                                                                                                                                                 <  1:165500/300‑1 (MQ=255)                                                                                                            ttgtcgttgctggAATTGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGGGAGAGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTGTGGCTGGTGCAGGAGGAAGGTTTTACTGACGTTGTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCTACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACAGATGACACCAGCAGCTGTGTATACGGATGATGCg                                                                                                                                                                                        <  2:151895/289‑1 (MQ=255)                                                                                                                      cGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCCGCAGCTGTGTATACGGATGATGCGGGTCGTCGAg                                                                                                                                                                              >  1:189990/1‑300 (MQ=255)                                                                                                                      cGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAg                                                                                                                                                                              >  1:190008/1‑300 (MQ=255)                                                                                                                                gTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTCTTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCCCGtggtcg                                                                                                                                                                    >  2:307375/1‑300 (MQ=255)                                                                                                                                gTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtcg                                                                                                                                                                    >  1:220782/1‑300 (MQ=255)                                                                                                                                   ggggggcggggatttgccagcgTTGAGGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACATCGCTTTGCTGGTGCAGGATGAAGGTCGTACTGACGTCTTGTACGTTCATCATTTTTCGCCTCGTGGTTCGTGTATCGGATGTGATGCTGAAGCGTCTTATTTGACCTATGGGGAGCGTCTTTGTAGGCCGAATGAGGCGTTGACGCCGCATCCGGCACCGGATCGAGTCTGCACAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGATAGTCGAGCCAGcggtcgggg                                                                                                                                                                 >  2:93143/20‑298 (MQ=255)                                                                                                                                   gAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGtt                                                                                                                                                                 >  2:92467/1‑300 (MQ=255)                                                                                                                                        ccacGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATNATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAccc                                                                                                                                                            <  1:566397/300‑1 (MQ=255)                                                                                                                                                 cGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTccac                                                                                                                                                   <  1:176391/300‑1 (MQ=255)                                                                                                                                                    ccGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCcaccac                                                                                                                                                <  1:61664/300‑1 (MQ=255)                                                                                                                                                    ccGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTACAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCcaccac                                                                                                                                                <  1:390275/300‑1 (MQ=255)                                                                                                                                                       cgTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGAGGAGCGCATTTGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCACCGGCTCAATTCTGCGAAACCGATGACACCAGCAGCTGTGTATNCGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTg                                                                                                                                             >  2:603636/1‑300 (MQ=255)                                                                                                                                                       cgTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGCCTTGCTGATGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGACTCCGGCACCGGCTCCACTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCTAGCACGCGGCCGGCTAACCCACCCTCCACCACTTg                                                                                                                                             >  2:19470/1‑300 (MQ=255)                                                                                                                                                           gACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGAACTAAGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCAATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCAACACTTGcccc                                                                                                                                         >  2:264454/1‑300 (MQ=255)                                                                                                                                                                     gAGGCGCAATTGAAGACGTGCAGACGCACTCCCTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCTTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcatca                                                                                                                               <  2:66423/300‑1 (MQ=255)                                                                                                                                                                             aTTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaaa                                                                                                                       >  1:708029/1‑300 (MQ=255)                                                                                                                                                                                       ggCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCca                                                                                                             >  1:186684/1‑300 (MQ=255)                                                                                                                                                                                                     gTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAAACGATGACACCAGCAGCTGTGTATACGGATGCTGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACg                                                                                               >  2:736500/1‑300 (MQ=255)                                                                                                                                                                                                           cTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTaaa                                                                                         <  1:396653/300‑1 (MQ=255)                                                                                                                                                                                                            tGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAAt                                                                                        <  1:677373/300‑1 (MQ=255)                                                                                                                                                                                                                   ggATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTAGTGATTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCCCCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGAGGGCGAAGCCTAAATCATGGGt                                                                                 >  2:440991/1‑300 (MQ=255)                                                                                                                                                                                                                     aTGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGa                                                                               >  1:251941/1‑300 (MQ=255)                                                                                                                                                                                                                                    actttttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGTATCCGGCACCGGCTCNATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAACACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGa                                                                <  2:208741/297‑1 (MQ=255)                                                                                                                                                                                                                                    actttttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGa                                                                <  2:631373/297‑1 (MQ=255)                                                                                                                                                                                                                                              cGTTAAGCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCACCACCGAGAGGTTCAGCTc                                                      <  2:308601/300‑1 (MQ=255)                                                                                                                                                                                                                                                     cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTccaccacc                                               <  1:715927/300‑1 (MQ=255)                                                                                                                                                                                                                                                      aTTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTccaccacca                                              <  1:347865/300‑1 (MQ=255)                                                                                                                                                                                                                                                         tttCGCCTCGTGGTTCTGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGGTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGc                                           <  1:469049/300‑1 (MQ=255)                                                                                                                                                                                                                                                         tttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGc                                           <  1:227239/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        cGCTGACGCGTCTTATTCGACCTACGGGGAGCGTATTTGTAGTCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCAc            <  2:307047/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        cGCTGACGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGNTTAACCCACTCTCCACCACTTGCCCCNGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCAc            <  2:555900/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             aCGCGTCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGaca       <  2:623866/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  tCTTATTCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATcc  <  1:540632/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                      TGCGCGCTGGCGCGGTGACCAGGTCTACCCACTCGTCACCGTGTTTGATCTGGATTTGACCTTCGTCGGGCAGATAGTTGGCGTACAGCGAGCGTAGCAGAGTTGATTTGCCGCTGCCGGAATGGCCGTGGAGCACCACGCATTCGCCCGCGTTGACGGTGAGCGAGGCGCGATTGAGGACGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCC  >  NC_000913/4317649‑4318237

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: ‑