Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,317	2 bp→AT	coding (339‑340/960 nt)	waaC →	ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,796,317	0	C	A	98.2%	188.1 / ‑4.5	56	coding (339/960 nt)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (19/36);  minor base . (0/1);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
*	NC_000913	3,796,318	0	A	T	96.4%	187.4 / ‑5.0	56	M114L (ATG→TTG)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (19/35);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACT  >  NC_000913/3796077‑3796558                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  tGACTGGGTGGTGTAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAATCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGTTAAAAAGCGCGGCGCGGGTGACGCGTCTGGCGCAGGGCGTAAAGCATTGATTGG‑ACTgg                                                                                                                                                                                                                                          <  2:434070/250‑1 (MQ=255)      gggTGGTGGAAGAAGGGTTCGCACAGATTCCGTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAAcc                                                                                                                                                                                                                                    <  2:190962/251‑1 (MQ=255)        gtggtgGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGc                                                                                                                                                                                                                                  >  1:128292/1‑251 (MQ=255)          ggtggAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTc                                                                                                                                                                                                                                <  2:326800/251‑1 (MQ=255)                    gggTGCGCACAGAGTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCTGCCCCCATAAAAGCTGAACGCATAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGTTGACGCGTCGGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAATCCGCTCGCGAAccgtt                                                                                                                                                                                                                      <  2:37955/251‑4 (MQ=255)                          gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGCTGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTagccac                                                                                                                                                                                                                >  2:330714/1‑250 (MQ=255)                          gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTagcc                                                                                                                                                                                                                  >  1:66211/1‑249 (MQ=255)                           cacaGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTagccagc                                                                                                                                                                                                               >  2:258296/1‑251 (MQ=255)                             caGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAATCCGCTCGCTAACCTTTagccagc                                                                                                                                                                                                               <  2:398355/249‑1 (MQ=255)                              aGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTg                                                                                                                                                                                                            <  1:207815/251‑1 (MQ=255)                                aTTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGtt                                                                                                                                                                                                          >  1:222317/1‑251 (MQ=255)                                 ttccttccTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGttt                                                                                                                                                                                                         >  1:93760/1‑251 (MQ=255)                                      tcctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACa                                                                                                                                                                                                     <  2:66211/250‑1 (MQ=255)                                       cctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACaa                                                                                                                                                                                                    <  1:258296/250‑1 (MQ=255)                                          ggCACGCTGCCGTTGATCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTTGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTGCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTTGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGTCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                                                                                                                                                <  1:173339/251‑1 (MQ=255)                                          ggCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                                                                                                                                                >  1:90213/1‑251 (MQ=255)                                                       ggAGCGAGTTATTCCTGTGGCAATACGTCGCTGGTGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAATCATCATATTg                                                                                                                                                                                   <  1:102405/250‑1 (MQ=255)                                                          gtgAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCaa                                                                                                                                                                                <  2:112926/249‑1 (MQ=255)                                                              gTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca                                                                                                                                                                            <  2:222317/251‑1 (MQ=255)                                                                             aTACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGa                                                                                                                                                             <  2:128292/251‑1 (MQ=255)                                                                               aCGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAAc                                                                                                                                                           <  2:22249/251‑1 (MQ=255)                                                                                    gCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAcc                                                                                                                                                      <  1:363267/251‑1 (MQ=255)                                                                                       ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTTG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAGTGCAAAACAGCAGCACGCCGTAGAGCGCAccc                                                                                                                                                     <  2:254669/249‑1 (MQ=255)                                                                                       ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCcg                                                                                                                                                    <  1:205151/250‑1 (MQ=255)                                                                                       ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACGATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAccc                                                                                                                                                     <  1:43721/249‑1 (MQ=255)                                                                                          gTAAAGCGTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCATTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGTCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCTCTCTCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGaa                                                                                                                                                <  1:267040/251‑1 (MQ=255)                                                                                                           gCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGtct                                                                                                                               >  1:415183/1‑249 (MQ=255)                                                                                                               ccATAAACGCTGAACGCAACGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGc                                                                                                                           >  1:318111/1‑251 (MQ=255)                                                                                                                      agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                                                                                     >  1:174758/5‑250 (MQ=255)                                                                                                                      agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                                                                                     >  1:75312/5‑250 (MQ=255)                                                                                                                            aCGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgcaaaccgc                                                                                                              >  2:347232/1‑251 (MQ=255)                                                                                                                               cAAAGCGAGTCGTAAAGCGCTACAAGCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAGAAAGCGGGGCGCTGGTGACGCGGCTGGCGCATGGCGTAACGCATGGATTGT‑ACTGGCAAACCGCACGGGAGCCGTGAGCCAGCCTGTGTGGTAATCGGAAGCATCATGTTGCAAATCAGGAGCACGCAGGAGAACGCACGCGCGAACGGTTGGCCAAAAGTGTGGGCTATTgcaatcggcaaa                                                                                                           <  2:23009/251‑1 (MQ=255)                                                                                                                                     gTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCACACCCAgg                                                                                                     <  1:401801/251‑1 (MQ=255)                                                                                                                                           ggAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGAtt                                                                                               <  2:157048/250‑1 (MQ=255)                                                                                                                                                 gcTACAAGCGGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAACAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCCAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCt                                                                                          >  2:181613/1‑250 (MQ=255)                                                                                                                                                  cTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAt                                                                                        >  1:335658/1‑251 (MQ=255)                                                                                                                                                    aCAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTAttg                                                                                      <  1:180347/251‑3 (MQ=255)                                                                                                                                                       aGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcac                                                                                   <  2:415183/251‑1 (MQ=255)                                                                                                                                                        gCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcaca                                                                                  <  2:376643/251‑1 (MQ=255)                                                                                                                                                            agaACTATGACGCAGTTATCGATGCACAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGGAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAGCCCAGGGCGATTATGCTATTGCACAGCAt                                                                              <  2:416528/251‑1 (MQ=255)                                                                                                                                                                  aTGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtttcct                                                                         <  2:287234/250‑4 (MQ=255)                                                                                                                                                                    gACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACTCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGa                                                                       <  1:187813/250‑1 (MQ=255)                                                                                                                                                                             aTCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGAAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGcc                                                             >  1:51684/1‑251 (MQ=255)                                                                                                                                                                                          ggCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCAGGGCGTAAAGCATGGATTGG‑ACTGGCAACGCGTTTGCGACCCGTTAGCCTGCCTGTTTTACAATCGTAAGGATCATATTGCAAAACAGCGGCACGCCGTAGGACGCTCCAGCGAACTGTTTTCCAAAAGTCTGGGCTATATCAATCCGTAGACCCAGTGCGATTGGGCTATTTCTCAGCAGTTCCTGTGGAACCTGCGTACGGATGCTGGc                                                <  2:193809/251‑1 (MQ=255)                                                                                                                                                                                           gCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGGTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATGATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCg                                               <  1:414378/251‑1 (MQ=255)                                                                                                                                                                                                    aaGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCg                                      <  1:330714/251‑1 (MQ=255)                                                                                                                                                                                                        gcggcgCTGATGACGCGTCTGGCGCATGGCGTACAGCATGGATTGG‑ACTGGCAAACCGCTCGCGCACCGTTAGCCAGCCTGTTTTACACTCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGACCTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCAGTAtt                                  >  2:359712/1‑251 (MQ=255)                                                                                                                                                                                                         cggcgCTGGTGACGCGTCTGGCGCATGGCTGAAAGAATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTAttt                                 <  2:196702/251‑1 (MQ=255)                                                                                                                                                                                                             gcTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATttcttc                             <  1:283845/251‑1 (MQ=255)                                                                                                                                                                                                                         tCTGGCGCATGGCGAAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGAAAGCCTGGTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACACGCGAACTGATTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGATATTGCACAGCATTTCCTGCCGAACCTGCCTACAGAAGCTGGCGAATCTGCCGAATTTCTTCATGCGAAGAccc                 >  1:334565/1‑251 (MQ=255)                                                                                                                                                                                                                            ggCGCATGGCGGAAAGCATGG‑ATGGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCATGTTTTACTATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACTGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATTCCGTATTTGGTTATGCGATGACCCg                <  2:366772/249‑1 (MQ=255)                                                                                                                                                                                                                               gcagggCGTAAAGCATGGAATGG‑ACGGGCAAACCTCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCATCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGtgatg           <  2:68302/247‑1 (MQ=255)                                                                                                                                                                                                                                      gTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAac     >  1:291820/1‑250 (MQ=255)                                                                                                                                                                                                                                       tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                             >  1:203011/1‑177 (MQ=255)                                                                                                                                                                                                                                       tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                             <  2:203011/177‑1 (MQ=255)                                                                                                                                                                                                                                        aaaGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACt  <  2:93760/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  TGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACT  >  NC_000913/3796077‑3796558

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