Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	419,111	A→G	15.4%	E408G (GAG→GGG)	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	419,111	0	A	G	15.4%	156.8 / 3.0	75	E408G (GAG→GGG)	phoR	sensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base A (23/37);  new base G (10/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 8.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGT  >  NC_000913/418887‑419338                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    gCCGGAAGGCACGCATATCCACGGACGCTGGCAGCGAGTGCCGCAAGGTGCAGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGagag                                                                                                                                                                                                                                  <  2:466867/228‑1 (MQ=255)  ccGGAAGGCACGCATAGCACCGTACGCTGGCAGCGAGTGCCTCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGagag                                                                                                                                                                                                                                  <  2:100109/227‑1 (MQ=255)   cGGAAGGCACGCATATCACCGTACGGTGGCCGCGAGTGCCGAACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTa                                                                                                                                                                                                                                <  2:137932/228‑1 (MQ=255)    ggAAGGCCACCATACCCCCGTACGCTGGCAGGGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCACACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTa                                                                                                                                                                                                                                <  2:241623/227‑1 (MQ=255)                   cACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTCCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCCGTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGAAGACGATTAGGGTTGCCGAGCGGGAAACACGCTGT‑GAATTGTCACGAAATGCGCCTGATTATTAGGGGGAAAGAAGGAAAAAGaaa                                                                                                                                                                                                                >  1:149227/1‑227 (MQ=255)                        tACGCTGCCAGCGAGTGCCGCCCGGAGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGGTAAAGCGCGT‑TCCCGGCAACCCGACGGTAGCCGATTGGGGTTTACTTTCGGGCAACATGCTGT‑TAAACATCACCAAAGTCGCATGACAATTACGAGGACAGTAAGAAAAGTAACACGtt                                                                                                                                                                                                           >  1:133563/1‑228 (MQ=255)                                 aGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTg                                                                                                                                                                                                  <  2:204976/228‑1 (MQ=255)                                                    ctaatgtaggtgtGAAGATAACGGGCCGGTCATTGGACCGGAGCATATTCCGCGCCTGACCGGACGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAATCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTaa                                                                                                                                                                                <  2:588681/215‑1 (MQ=255)                                                    cGAATGTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTaa                                                                                                                                                                                <  2:618826/227‑1 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGCTAAAGCGCGT‑TCCCGGCCAACCGGCGGTACCGGATTAGGGTTAGCGATCGTGCAAAGTGCTGG‑CAATCGTCACCGAAGTCGCCTCAATATTTGGAGGACATTAGTAAAATGGACACGTTTCAGATTTGTTTTCCCgggaattttgatttcccaa                                                                                                                                                                        >  1:376754/1‑211 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGCTAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGGTAGAGATCGTGTAACATGCCGT‑GAATCATCACGGAAGTCGCCTGGATATTGAGGGTAAAGTAGGGAAAGGAACAAGTTTGAGTTTTGTTATCCCGGAAAGCTTAAATGACaaa                                                                                                                                                                        >  1:287303/1‑228 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGGTTAGGGTTAGCGATGGTGAAACATGCTGT‑GGATCCTCACGGAAGTCGCCTGAATATTGAGGGTACAATCGGGAAAGGAAACCGGTTGAGTTTTTTTATCCCGGGACGATTAAATGACaaa                                                                                                                                                                        >  1:93411/1‑228 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCGTCACGACAGTCGCCTGGATATTGATAGTACAGTAGGAAAATGAACCCGGTCAATTTTTGTTATCCCGGGACGTTTAATTGCCaaa                                                                                                                                                                        >  1:534788/1‑228 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGGTTAATTGCCaaa                                                                                                                                                                        >  2:223948/1‑228 (MQ=255)                                                           aGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGGGGAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGTGAGTACAGGAGGAAAAGGAACACGTTTCAGCTTTGTTATCCCGGAAAGTTTAATGGCCaaa                                                                                                                                                                        >  1:98706/1‑228 (MQ=255)                                                                  aaGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGGGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGGACAGTAGGGAAAGGAACACGTTTCAGTTTTTTTATCCCGGAACGTTTAATTGCCAAAAACAg                                                                                                                                                                   >  2:576167/1‑226 (MQ=255)                                                                       aaCGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTGGGGTTAGCGGTCGTGACAAATCCTAT‑GAACCATCACGGAAGTCCGCTGAATATTGGTAGGACAGAGGGGCAAGGAACACAATTCGAGTTTGTTATCCCCGCAAGATTATTTGCCaaaaaacagtctct                                                                                                                                                             >  1:94482/1‑218 (MQ=255)                                                                        aCGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTaa                                                                                                                                                                                <  2:94482/207‑1 (MQ=255)                                                                               gggCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAACGCGCGT‑TCCCGGCAAACCGGCGGGAGCGGGTTTAGGGTAGCGATCGTGAAACAGGCTGT‑GAATCATCACGGAAGTCGGCTGGATATTGTGAGTAAAGTAGGGAAAGGACAACGATGCAGATTTGTTATCCCGGAAGTGTTAATTGCCAAAAAAAGCGCATAATCCTCCtt                                                                                                                                                    >  1:641417/1‑228 (MQ=255)                                                                               gggCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTa                                                                                                                                                                                                                                <  2:95796/152‑1 (MQ=255)                                                                               gggCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTCTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGCAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAAATGACAAAAAAACCGATTAAACCACCtt                                                                                                                                                    >  1:107028/1‑228 (MQ=255)                                                                                 gCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCAACACGAAAGTCGCCTGAATATTGGGAGTAAAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACATGGATTAATCCGCCTTTg                                                                                                                                                  >  2:410647/1‑228 (MQ=255)                                                                                    ttGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTCTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCa                                                                                                                                               <  2:107028/228‑1 (MQ=255)                                                                                     tGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGGACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTTTCAt                                                                                                                                              >  2:25146/1‑228 (MQ=255)                                                                                           ggAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGATTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTGGTCATCtttt                                                                                                                                         >  2:473162/1‑227 (MQ=255)                                                                                           ggAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGGCGCCTGGATATTGGGGGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATACCGGAACGTTTAATTGCCCAAAACAGCGATTAATCCGCCTTTGGCATCtttt                                                                                                                                         >  2:49987/1‑227 (MQ=255)                                                                                            gacgattttCCGCGCCTGGCCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTa                                                                                                                                        <  2:4524/220‑1 (MQ=255)                                                                                            gAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGGGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGGATATTGGGGGTACAGTAGGAAAAGGAACACGGTTCAGGTTTTTTATCCCGGAACGTTTAATTGCAAAAAAAAGCGATTAATCCGGCTTTGTCATCTTTTAt                                                                                                                                       >  1:79665/1‑228 (MQ=255)                                                                                             gtcATATTCCGCGCCTGCCCGAGCTTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTAt                                                                                                                                       <  2:373930/225‑1 (MQ=255)                                                                                             aggtaatTCCGCGCATGACCGGGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGGGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACCCGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTAt                                                                                                                                       <  1:382533/222‑1 (MQ=255)                                                                                             agagtgttCCGGGCCTGAGCGAGGGTTTTGTGCGCGGTGATAAAGCGCGG‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAAGCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTAt                                                                                                                                       <  1:93443/221‑1 (MQ=255)                                                                                               cATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGGTTAGGGTTAGCGATCGTGAAACATGCTGT‑GAAACATAACGAAAGTCGCCTGGAAATTGGGAGCACAGGAGGAAAAGGAACACGTTTCCGTTTTTTTATCCCGGAAACTTCAATTGCCAAAAAAAACGATAAATCTGCCTTTCTCAAtttttttttc                                                                                                                                    >  1:631074/1‑226 (MQ=255)                                                                                                atatTCCGCGCCTTACCGAGCGTTTTTATCGGGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGcc                                                                                                                                   <  2:526219/228‑1 (MQ=255)                                                                                                   ttCCGCGCCTGGCCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCAt                                                                                                                                 <  2:532047/227‑1 (MQ=255)                                                                                                   ttCCGCGCCTGACCGAGCGTTTTTAACGCGTTGAGAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCAGCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGtt                                                                                                                                                                                                <  2:594196/164‑1 (MQ=255)                                                                                                    tCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGGTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATa                                                                                                                                >  2:147221/1‑227 (MQ=255)                                                                                                    cccccgcctgcgggggCGTTTTTATCGGGTTGAAAAAG‑GCGTGTCCCCGCAAACCGGCGGTAGCGGATTAGGGTTAGGGATGGTGAAACATGCGGT‑GAAGCATCACGAAAGTCGGCTGAATAGTGAGAGTCTAGTAGGAAAAGGAACACGTGTCAGTTTTTGTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCAtca                                                                                                                               <  1:667777/213‑3 (MQ=255)                                                                                                              aCCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTTGCGGTTTAGGGGTAGCGGTCGTGAAACATGCTGC‑GAATCATCACGGAAGTCGCCTGAATATTTGGAGTACAGCTGGCAAAGGAAACCCATTCCGGTTGTTTTTCCCGGAACGTTTAATTGACAAAAAAGACGATTAATTCGCCTTTTGCAACTttttttatcaaaaaacaaacacg                                                                                                                     >  1:310418/1‑211 (MQ=255)                                                                                                              aCCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTTGCGGATTAGGGGTAGTGACCGGGAAACATGCTGG‑GAAACATCCCGAAAGTAGCCTGAATATTGGGCATACAGTAGGAAAAGGAACCAGATTCCGTTTTGTTATACCGGAAAGTTTAATTGGCAAAAAAAGGGATTAATTCTCCTTTGTAATCTTTTATTTCCATAAACAAATCGac                                                                                                                     >  1:14221/1‑227 (MQ=255)                                                                                                              aCCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTTTTGCCATAAGCCAGCCgct                                                                                                                     >  2:442382/1‑226 (MQ=255)                                                                                                               ccGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGGGGTAGCGGATTAGGGTTGGCGATCGTGAAACATGCTGG‑GAATCATAGCGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATg                                                                                                                    <  2:470980/228‑1 (MQ=255)                                                                                                                 gggCGTTTTTATCGCGTTGAAAAAGCGCGT‑TCCCGGAAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTTAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCt                                                                                                                  <  2:631074/226‑1 (MQ=255)                                                                                                                   gggTTTTTTTCGGGTTGATAAAGCGCGT‑TCCCGGCAAACGGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTgg                                                                                                                <  1:644280/226‑1 (MQ=255)                                                                                                                                 ttgaaaaaggcggT‑TCCCCGCAATCGGGCGGTAGCGGATGTGGGTTAGGTGTCGTGAAACATGCTGT‑AAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTg                                                                                                   <  1:121607/215‑1 (MQ=255)                                                                                                                                 ttGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGTTTAGGTTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGACTATTGGGAGTACAGTAGGGAAAGGAAAACGTTTCAGATTTGTTATTCCGGAACGTTTAATTGCCAAAAACAGCGGTTAATCCGCCTTTCTCACTTTTTATTTCCATAAACCAATACATGCTTGTTTCTTTTCTTTgt                                                                                                  >  1:221594/1‑227 (MQ=255)                                                                                                                                 ttGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGACTATTGAGAGTAAAGTAGGGAAAGGAACACGTTTCAGTTTTGTTATCCCGGGACGTTTAATTGCCGAAAAAAGCGATTAATCCGCCTTTTCCAGCTTTTTTTGCCCTAAACCAGGCGGTGCTGGATTAtttttttttc                                                                                                  >  1:408279/1‑226 (MQ=255)                                                                                                                                    aTAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGt                                                                                               <  2:476866/228‑1 (MQ=255)                                                                                                                                            cgT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGGGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTTTCATCTTTTATTTGCATAAGCCAATCGATGCTGGCTttttttttttttattctaaaaaa                                                                                       >  1:422917/1‑215 (MQ=255)                                                                                                                                             gT‑TCCCGGCAAACCGGCGGTAGCGGATTAGTGTTAGCGATCGTGAAACATGCTGT‑GAATCATCCCGAAAGTCGCCTGAATATTGAGAGTACAGTAGAAAAAGTAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATAc                                                                                       <  2:98706/227‑1 (MQ=255)                                                                                                                                              t‑tCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAAAATTGAGAGTACAGTAGGGAAAGGAACACGTTTCAGTTTTGTTATGCCGGAACGTTTAATTGCCAAAAAAAGCGATTAATCTGCCTTTGTGATCCTTTATTTCCATAAACCAGCCGCTGCTGGCTCTTTTTTTTTCCATTACAAATaatg                                                                                     >  1:548866/1‑225 (MQ=255)                                                                                                                                                  ccgggaatacggcagTAGCGGGTTAGGGTTAGCGTTCGTGAAACATGGTTT‑GAATCATCACGAAAGTCGCCTGAATTTTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGACAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCg                                                                                  <  2:133563/219‑1 (MQ=255)                                                                                                                                                  ccGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGG‑TAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCg                                                                                  <  2:287303/228‑1 (MQ=255)                                                                                                                                                             tcgggtgtgtaTTGGGGTTAGGTATCGTTAAAATTGTTGT‑GAATTGTTACGAAAGGCGCTTGGATAGTGAGAGTACAGTAGGAAAAAGAACAAGTTTCAGTTTTGTTAGCCCGGAACTTTTAATTGCCAAAACCAGCGATTAATCCGCCTTTGTCATTTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTATGCAGTCAAAATACGGGCGTTAGATTTTAc                                                                       <  2:318511/218‑1 (MQ=255)                                                                                                                                                                           gggTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCATTAAGCCACTCGATGCTGGCTTATTTTCTTTCTATTCAAAATCCGGGCGTTAGATTTTACAACGTTTGGTGAtt                                                         >  1:100429/1‑228 (MQ=255)                                                                                                                                                                                 gCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTGGCCATTGTTTTTTTCCATAAGCCAGTCGATGTTGGCTTATTTTTTTTGAAGTTTAAATACGGGGCTTAGATTTTACAACGATGTGGGTtttttttt                                                   >  1:483050/1‑226 (MQ=255)                                                                                                                                                                                    tttgTGAAACATGCTGG‑AAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCTTCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTc                                                 <  1:25146/224‑1 (MQ=255)                                                                                                                                                                                    aTCGTGAAACATG‑TGTGGAATCTGCACGGAATTCGCCTGAATATTGAGAGGACATGAGGAAAAGGAACCCGTTTCAGTTTTGTTATCCCGAATGGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTTCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTc                                                 <  1:236567/227‑1 (MQ=255)                                                                                                                                                                                      cGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTgaga                                                                                                                                                                                                                                   >  1:179888/1‑47 (MQ=255)                                                                                                                                                                                       gTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGagag                                                                                                                                                                                                                                  <  2:179888/47‑1 (MQ=255)                                                                                                                                                                                          aaaCATGCTGT‑GAATCAGCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGAt                                          <  2:408279/228‑1 (MQ=255)                                                                                                                                                                                               tGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACACCGATTGGTGATTTTTTGTTCGCATAATTAGc                                      >  2:249207/1‑227 (MQ=255)                                                                                                                                                                                               tGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTTTTCGCATGATTAGcc                                     >  2:696709/1‑228 (MQ=255)                                                                                                                                                                                                       aaTCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGGGATTTTTTGTTCGCTTGATTAGCCATGTCtt                              >  2:426090/1‑228 (MQ=255)                                                                                                                                                                                                       aaTCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAAAACGTTTCAGTTTTGATATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCTCCTTTGTCATCTTTTATTGCCATAAGCAAGTCGTTGCTGGCTGATTTTATTTGCAGTTAAAATACGGGCGTTACATTTTACAACGATTGGTGCTTTTTTGTTGGCAtatatggacaatgctt                              >  1:399145/1‑213 (MQ=255)                                                                                                                                                                                                          cctcaCGAAAGTTGCCTGAATATTGAGAGTACAGGAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGTCTTTGTCATCTTTTATTGCCATAAGCCAGTCTATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCttttt                           <  2:380527/226‑1 (MQ=255)                                                                                                                                                                                                            tcccGAAATTCGCCTTAATATTGAGATTACAGTAGGAAAAGGAACCCGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCa                         <  1:49987/225‑1 (MQ=255)                                                                                                                                                                                                              aCGAAAGTCAGCTGAATATTGAGAGTTCAGTAGGAAAAGGAACACGTTTCCGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTAATGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCa                         <  2:243981/226‑1 (MQ=255)                                                                                                                                                                                                               cGAAAGTCGCCTGAATATTGAGAGTACAGTGGGAAAAGGAACACGTTTCAGCTTTGTTATCCCGGCACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCACCTTTTTTTGTCATAAGCCAGTCGATGCTGGCTTATTTTTTTTGCAGCCTAAATACGTGCGTTAGAGTTTACAACGATTGGTGTTTTTTTGTTTTCGTGTTTAGCCATGACTTTTTCACgg                      >  1:75622/1‑228 (MQ=255)                                                                                                                                                                                                                  aaGTCGCCTGAATTATGGGGGTTAAGGAGGAAAAGGAACACTTTTTAGTTTTGTTATCCCGGAACGTTTAATTGGCAAAAACAGCGATTAATCCGCCTTTTTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGaaa                   <  1:442382/228‑1 (MQ=255)                                                                                                                                                                                                                    ggcGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  2:34847/226‑1 (MQ=255)                                                                                                                                                                                                                    gTCTCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGGTTCAGTTTTGTTATCCCGGAACGTTTAATTTCCACAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCAATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCTTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  2:65546/228‑1 (MQ=255)                                                                                                                                                                                                                    gTCGCCTTAATTTGGAAAGTACAGTAGGAAAAGGAACCCGTTTCAGTTTTTTTATCCCGGAACGTTTAAATGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  1:473162/228‑1 (MQ=255)                                                                                                                                                                                                                    gTCGCCTTAATATTGGAATGACAGTAGGAAAAGGAACCCGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  2:550781/228‑1 (MQ=255)                                                                                                                                                                                                                    gTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACACCGATTGGTGATTTTTTGTTCGCATGGTTAGCCATGTCTTTTCCAGGGAAATa                 >  1:587012/1‑228 (MQ=255)                                                                                                                                                                                                                    gTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCCTGATTAGCCATGTCTTTTTCACGGaaa                   >  2:312232/1‑226 (MQ=255)                                                                                                                                                                                                                    gTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  2:420320/228‑1 (MQ=255)                                                                                                                                                                                                                    gTCGCCTGAATATTGAAAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATa                 <  2:442700/228‑1 (MQ=255)                                                                                                                                                                                                                        ccTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGCGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATatggt             >  2:157911/1‑224 (MQ=255)                                                                                                                                                                                                                              tatTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGGTATCCCGGAACGTTTAATTGCCAAAAGCAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTtata       <  2:14221/228‑1 (MQ=255)                                                                                                                                                                                                                                  gagaGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTTTTTTTCACGGAAATAGTGtttttttttgg   >  1:472967/1‑221 (MQ=255)                                                                                                                                                                                                                                   agagTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGCTGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTTCAACGATTGTTGTTTTTTTTTTGCAATGATTTGCATTTTCTTTTTCACGCAATTATGGTTTTTTACtggt  >  1:696938/1‑228 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    GCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTGAAGATAACGGACCGGGCATTGCACCGGAGCATATTCCGCGCCTGACCGAGCGTTTTTATCGCGTTGATAAAGCGCGT‑TCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGT‑GAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGT  >  NC_000913/418887‑419338

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