Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	418,428	Δ1 bp	50.0%	coding (540/1296 nt)	phoR →	sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	418,428	0	G	.	50.0%	72.7 / 82.0	42	coding (540/1296 nt)	phoR	sensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base . (10/11);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCT  >  NC_000913/418145‑418709                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                          aTAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      <  1:61320/289‑1 (MQ=255)                        ggCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGAtgg                                                                                                                                                                                                                                                               >  1:270614/1‑289 (MQ=255)                             tCTGATTAAACGCTTTCGTAGCGGCGCGTAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTGGCGCGCAACCAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTCCGTTACCCGGAGTTTACGCAATATCTGCAAATGTGTGATTTTTCTCGCCCGCCCAATCTGGTTCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACACACAAAACTTTGCTGATGGTGgc                                                                                                                                                                                                                                                           >  2:173766/1‑289 (MQ=255)                                   tAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATg                                                                                                                                                                                                                                                   <  1:20505/289‑1 (MQ=255)                                           ttCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCa                                                                                                                                                                                                                                            >  1:305813/1‑289 (MQ=255)                                                                                    aCCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCa                                                                                                                                                                                                   >  2:235915/1‑289 (MQ=255)                                                                                                ggCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATg                                                                                                                                                                                       >  1:99069/1‑289 (MQ=255)                                                                                                ggCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAGACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGCGCCATGa                                                                                                                                                                                      >  1:318333/1‑289 (MQ=255)                                                                                                                 tAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCa                                                                                                                                                                                                                                                <  2:278710/214‑1 (MQ=255)                                                                                                                 tAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCa                                                                                                                                                                                                                                                >  1:278710/1‑214 (MQ=255)                                                                                                                  aaCGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGa                                                                                                                                                                     <  2:270614/289‑1 (MQ=255)                                                                                                                  aaCGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGAc                                                                                                                                                                    <  1:110587/289‑1 (MQ=255)                                                                                                                    cGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGAcc                                                                                                                                                                   <  1:235915/289‑1 (MQ=255)                                                                                                                        cTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGgc                                                                                                                                                                                                                       <  2:67510/233‑1 (MQ=255)                                                                                                                        cTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGgc                                                                                                                                                                                                                       >  1:67510/1‑233 (MQ=255)                                                                                                                         tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      >  2:213592/1‑169 (MQ=255)                                                                                                                         tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      <  1:213592/169‑1 (MQ=255)                                                                                                                         tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTtata                                                                                                                                                                                                                                                                                    <  1:110777/170‑1 (MQ=255)                                                                                                                         tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTtata                                                                                                                                                                                                                                                                                    >  2:110777/1‑170 (MQ=255)                                                                                                                                 aacaaATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAggg                                                                                                                                                        <  2:58311/288‑1 (MQ=255)                                                                                                                                  acaaATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGtt                                                                                                                                                      >  2:120496/1‑289 (MQ=255)                                                                                                                                        tCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTg                                                                                                                                                 <  2:99069/289‑1 (MQ=255)                                                                                                                                                   gcTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGatgaat                                                                                                                                     <  1:120496/289‑1 (MQ=255)                                                                                                                                                                 aaCGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAACATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGa                                                                                                                       >  1:303360/1‑289 (MQ=255)                                                                                                                                                                      gCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGACGCTGGAAGgc                                                                                                                   >  1:39776/1‑288 (MQ=255)                                                                                                                                                                           aCATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGt                                                                                                              >  1:101059/1‑289 (MQ=255)                                                                                                                                                                                                   cgaggTTTACCCACTATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAACCTGGTGCTCAACCCCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAGACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATg                                                                                     <  2:318333/285‑1 (MQ=255)                                                                                                                                                                                                      aGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAg                                                                                                   <  1:69998/273‑1 (MQ=255)                                                                                                                                                                                                      aGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAg                                                                                                   >  2:69998/1‑273 (MQ=255)                                                                                                                                                                                                                atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGAcc                                                                         <  2:305813/289‑1 (MQ=255)                                                                                                                                                                                                                atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGa                                                                           <  2:304756/286‑1 (MQ=255)                                                                                                                                                                                                                atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGa                                                                           >  1:304756/1‑286 (MQ=255)                                                                                                                                                                                                                           cgcgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    <  2:231654/250‑1 (MQ=255)                                                                                                                                                                                                                           cgcgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    >  1:231654/1‑250 (MQ=255)                                                                                                                                                                                                                             cgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGTAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCgg                                                                   >  1:219023/1‑282 (MQ=255)                                                                                                                                                                                                                             cgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCgg                                                                   <  2:219023/282‑1 (MQ=255)                                                                                                                                                                                                                                       ctcGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAg                                                 <  2:303360/289‑1 (MQ=255)                                                                                                                                                                                                                                                   aTCTGGTGCTCAGCAACGGGCGGCAACTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTGGCTGGGGGTGGCGCGGGATGGCACGGAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACg                                     <  2:296585/289‑1 (MQ=255)                                                                                                                                                                                                                                                          gCTCAACACCGGGCGGCATCTGGAAATTCGCGTCAT‑CCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCAAGAGGATGGAAGGACTGGTGAAGCCATTGCTGACGCTGTCGa                              >  2:108864/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                    ttCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    <  1:236614/194‑1 (MQ=255)                                                                                                                                                                                                                                                                                    ttCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    >  2:236614/1‑194 (MQ=255)                                                                                                                                                                                                                                                                                      cgcgTCA‑GCCTGAATCCCCCAAGCAGTGGCGGAGGGTGGAGCGTGATGTGACGAGAATGCATAAACGGGAAGGGACGCGGCGTAACGTATTTGCCAACGTGAGCCATGAGTTACGTACGCCAGTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAAGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCt  <  1:173766/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                          ATAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCT  >  NC_000913/418145‑418709

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