Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,344	T→G	P122P (CCT→CCG)	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,344	0	T	G	98.2%	180.5 / NA	55	P122P (CCT→CCG)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (20/34);  total (20/35)
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.89e-01

GCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATT  >  NC_000913/3796102‑3796583                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                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:75312/5‑250 (MQ=255)                                                                                             agctgAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAgc                                                                                                                                              >  1:174758/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)                                                                                                                                                                                                                    atggattgg‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACGCTGGCCg                      <  2:270425/245‑1 (MQ=255)                                                                                                                                                                                                                                tGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATTATAAACACTGGCCGGAAGACCACTg           >  2:417850/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaaccccaagcgAACCGTTAGCCAGCCTGTTTTACAACCGTAAGCATCATATTGCAACACCGCAGCACGCCGTAGAACGCACCCGCGAACTGCTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCCGACGAACCTGACTTCCGATGCTGGCGAATATGACGTATTTCTTCATGCGACGACCCGTGATGATAAACAATTGCCGGAAAAaccattgggg       >  2:434817/10‑243 (MQ=255)                                                                                                                                                                                                                                         gCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCTAGAatt  <  2:69861/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                GCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATT  >  NC_000913/3796102‑3796583

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