Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,149,252	G→T	5.6%	T338N (ACC→AAC)	eptC ←	LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,149,252	0	G	T	5.6%	197.5 / 2.7	61	T338N (ACC→AAC)	eptC	LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base T (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTCCAGCTCTTCCGCGTTTAATCCCGGCGGAACATGATCGGTATTGCCATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCG  >  NC_000913/4148993‑4149520                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                             cTCCAGCTCTTCCGCGTTTAATCCCGGCGGAACATGATCGGTATTGCCATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGTCTTCAGCACGTTGGTGTCGTATTCACCCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGACA‑GCATGGTATTGCGGGCGGTCa                                                                                                                                                                                                                                                               >  1:452635/1‑275 (MQ=255)    cAGCTCTTCCGCGTTTAATCCCGGCGGAACATGATCGGTATTGCCATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTTTCGGTCTGGCGCCAAAAT‑CCCTCACGCATGGTATTGCGGGCGGTCATCg                                                                                                                                                                                                                                                            >  2:232795/1‑275 (MQ=255)                                   ttaTCGGTATGGCCAGCAAACTTGCCCTGATTTTCCTGGTAGCGGTATTTGTATGTGATATGCGTACCCAGCAGATGAACAATGTTCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGGCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCATAAGGTTTtata                                                                                                                                                                                                                             <  1:232796/273‑1 (MQ=255)                                                    aaaCTTGCCCTGATTTTCCGGGTAGCGGAATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGAACGTTGGTGTCGTATTCACGAGCACTCTGCGTACGTTGCTGGTTCATGTAGTCCTGCTTGTCGGTCTGGCGCGAACATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTTTGCTGGTTGGTGATACCGACGGTTTTATAAACAGCCTGTTtcatta                                                                                                                                                                                                            >  1:457517/1‑273 (MQ=255)                                                         tGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCa                                                                                                                                                                                                                                         >  1:89629/1‑241 (MQ=255)                                                          gccggtaTTTGCCGGGTAGCGGTAGTGGTATTGGAGATGCGTAGCCGGGAGAGGACAAATGATGAGTGTCTGCGGCGGAGGGTCATTCAGCACTTCCTGGAAGAGCTTCAGCACTTGGGGGTCGTATGAACGGGCACTCGGCGGACGTTGCTGGTTCATGTAGTACGGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTGTCATCATGTtca                                                                                                                                                                                                      <  2:49572/269‑1 (MQ=255)                                                             cTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAg                                                                                                                                                                                                                                     <  2:89628/241‑1 (MQ=255)                                                                    tCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGTTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGAc                                                                                                                                                                                              >  1:496856/1‑273 (MQ=255)                                                                      cGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCt                                                                                                                                                                                          <  2:199569/275‑1 (MQ=255)                                                                        ggTGGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGTTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCt                                                                                                                                                                                          <  2:496855/273‑1 (MQ=255)                                                                                       atttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTg                                                                                                                                                                                         >  1:72831/1‑259 (MQ=255)                                                                                          tgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTctgc                                                                                                                                                                                                                                                       >  2:167092/1‑194 (MQ=255)                                                                                           gATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTc                                                                                                                                                                                     <  2:72830/259‑1 (MQ=255)                                                                                              atGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCtggt                                                                                                                                                                                                                                                   <  1:167093/194‑1 (MQ=255)                                                                                                        agcagATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAACTACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGCAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCAtttg                                                                                                                                                         >  2:472914/1‑272 (MQ=255)                                                                                                          cagATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTggtac                                                                                                                                                      >  1:405208/1‑272 (MQ=255)                                                                                                          cagATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCAGCAGCGACGGCTGCGTTAGATACAGATCCGGGTTCTTTTCATTGGCaa                                                                                                                                                      >  1:279832/1‑275 (MQ=255)                                                                                                             aTGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCCTGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATATGCCAAgg                                                                                                                                                   >  2:141633/1‑275 (MQ=255)                                                                                                               gAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCCTGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTc                                                                                                                                                 >  1:190200/1‑275 (MQ=255)                                                                                                                       aTCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTTCTGGTTCATGTAGTACTGCTCGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTCTTGCGGGCGGGCATCGTCTGCTGGTTGGTGACCCAGATGCTTTTCTAACCCGCCTGTTTCATAATGTTCCTCAACGACGCCCGCGTCCGATACAGATACTGGCTCTGTTCACTGTCAAAGGTACGCGCCTg                                                                                                                                         >  2:444773/1‑275 (MQ=255)                                                                                                                       aTCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTg                                                                                                                                         >  1:253167/1‑275 (MQ=255)                                                                                                                                 tcGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGAtt                                                                                                                               >  1:151991/1‑275 (MQ=255)                                                                                                                                     cgcAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGCTTTCaa                                                                                                                           >  1:406386/1‑275 (MQ=255)                                                                                                                                      gcAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCaa                                                                                                                           <  2:253166/274‑1 (MQ=255)                                                                                                                                      gcAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCaa                                                                                                                           <  2:151990/274‑1 (MQ=255)                                                                                                                                          ggTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGgtg                                                                                                                      >  1:102289/1‑275 (MQ=255)                                                                                                                                                tCAGCCCGTCCGGGAACGGGTTCAGGACTTGGGGGTCGTATTCACGCGCACTCTGCGTACGTTGCGTGTTCATGTTGTACGGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTGTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGa                                                                                                                <  2:452634/275‑1 (MQ=255)                                                                                                                                                gcatcACTCCCTGGAACGGCTTCAGCACGTTGGTGTGGTGTTAACGCGCGCTCTGCGTACGTTGCTGGTGCATGTAGTACTGCTTTTGGGTCTGGCGCGAAAATACCGTCA‑GCATGTTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAGCCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGa                                                                                                                <  2:339564/271‑1 (MQ=255)                                                                                                                                                    cACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGCTCCAGAAGGTTTCATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCAGCGTCAGACACATCTCAGGGTGTGTTTCTTTACGAAACGTCCGCGCCGGGTGCACGACTTTAAGCGTGTCCgagcctg                                                                                                            >  2:126243/1‑269 (MQ=255)                                                                                                                                                      cTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCGGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGGCAGGTAGAAACCCGGGTTCTTTTCATTGCAACAGGTCAGCGCCCGTTGCAGGATTTCAATGGTGTATGGGCGAGAc                                                                                                          >  2:362995/1‑275 (MQ=255)                                                                                                                                                         ccTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCCGTTGCAGGATTTCAATGGTGTACGGACGAACCGTa                                                                                                       >  2:494558/1‑275 (MQ=255)                                                                                                                                                             gAACGGGTTCAGCTCGTTGGTTGCGTATTCACGCGCACTCGGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGGTGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTa                                                                                                   <  2:406385/275‑1 (MQ=255)                                                                                                                                                               aCGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACg                                                                                                 <  1:494559/275‑1 (MQ=255)                                                                                                                                                                 ggCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCCAAGGTCAGCCGCTGCTGCAGGATTTCACTGGTGTACGGATGGGCCGTACCTACGtt                                                                                               >  2:208179/1‑275 (MQ=255)                                                                                                                                                                             ttGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATACAGCAGGTTTTATAACCCGCCTGATTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATACGGGTTCTTTTCATTGGCAACGGTCTGCGCCTGTTGCAGGCTTTCAATGGTGTCCGGACGAGACGTa                                                                                                       >  1:255212/1‑255 (MQ=255)                                                                                                                                                                             tgtgttgtgtgttcccgcgcctgctGCGTCCGTTGGTGGTTGAGGTAGTACTGCGTGTCGGTCTGGGGCGAAAATACCTTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGTTCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGt                                                                                   <  2:279831/252‑1 (MQ=255)                                                                                                                                                                                 gttCGTATTCCCGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTa                                                                                                   <  2:255211/253‑1 (MQ=255)                                                                                                                                                                                          ccgcggcACTCCGCGTACGTTGTTGGTTCATGTAGTACTGCTTGTCTGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCCGTCATCGTCTGCTGGTTGGTGAGCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCg                                                                      <  1:260126/270‑1 (MQ=255)                                                                                                                                                                                                cACTCTGCGTACGTTGCTGGTTTATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGtttt                                                                 >  2:474513/1‑274 (MQ=255)                                                                                                                                                                                                cACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGtt                                                                                                                                       >  1:325041/1‑204 (MQ=255)                                                                                                                                                                                                  ctctGCGTACGTTGCTGATTCATTTAGTACTGCTTGTCGGTCTAGCGCGAAAATACCGTCA‑GCATGTTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCATAAGGTATTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGAtt                                                                                                                               >  1:219953/1‑210 (MQ=255)                                                                                                                                                                                                    ctGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAg                                                                                                                                   <  2:325040/204‑1 (MQ=255)                                                                                                                                                                                                      gcttACGTTGCGGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCaa                                                                                                                           <  2:219952/207‑1 (MQ=255)                                                                                                                                                                                                               gCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGCCGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTTCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGCAAGATTCGGATCGGTTTTATGCAGGGCCTCCAGc                                                 >  2:300000/1‑275 (MQ=255)                                                                                                                                                                                                               gCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTACGATTCGGATCGGTTTTATGCAGCGCATCAAGc                                                 >  1:515497/1‑275 (MQ=255)                                                                                                                                                                                                                             taCTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGCCGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAAGTACGTTATTAAACACGGCCAGTTTCGGATCGGTTTTATGAAGCGCATCCAGCCCCGGCGTGGTTTc                                   >  2:177563/1‑275 (MQ=255)                                                                                                                                                                                                                                  cTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTTTAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCAcgcg                              >  2:454195/1‑275 (MQ=255)                                                                                                                                                                                                                                  cTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGAtt                                                                                                                               >  1:488638/1‑178 (MQ=255)                                                                                                                                                                                                                                  cTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCAcgc                               >  1:476741/1‑274 (MQ=255)                                                                                                                                                                                                                                  cTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGAAGCGCAACCAGCCCCGGCTTGGTTTCAcgcg                              >  2:512967/1‑275 (MQ=255)                                                                                                                                                                                                                                   ttGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCAcgcg                              <  1:474514/274‑1 (MQ=255)                                                                                                                                                                                                                                     gTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACAGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCaag                                >  1:247789/1‑268 (MQ=255)                                                                                                                                                                                                                                     gTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACAGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGa                            <  2:247788/274‑1 (MQ=255)                                                                                                                                                                                                                                      tgggTCGGGCGTGAAAATCCCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATGACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATa                          <  2:190199/273‑1 (MQ=255)                                                                                                                                                                                                                                      tCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCaa                                                                                                                           <  2:488637/178‑1 (MQ=255)                                                                                                                                                                                                                                       ccgTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGGGGTCCTCGGCTGCTGGTTGGTTATCCAGGAGATTTTATACCCCGCCTGTTGCATCATGTTCATCAGCGACCGCTGCTTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGGACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCTGCGCATCCAGCTCCCGCGTGGTTTCACGCGGATaa                         <  1:362996/273‑1 (MQ=255)                                                                                                                                                                                                                                       cGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGAGTGGTTTCACGCGGATaa                         >  1:285340/1‑275 (MQ=255)                                                                                                                                                                                                                                       cGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGAAGTAACTACGTTATTAAACACGGTCAGATTCGGCTCGTTTTTATGCAACGCATCCAGCTCCGGCCTGGTTTCACGCGGATaa                         >  1:17774/1‑275 (MQ=255)                                                                                                                                                                                                                                            cggctgcaaaaTACCGTCA‑GCATGGTATTGCGGGCGGTCAGCGTCTGCTGGTTGGTGATCCAGAATGTTTTATAACCCGCCTGTTTCATCATGTTCCTCAGCGACGGCTGCGTCAGATACAGATCCTGGGTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTa                    <  2:125002/268‑1 (MQ=255)                                                                                                                                                                                                                                                              a‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTTTAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCcgcg  <  1:454196/275‑1 (MQ=255)                                                                                                                                                                                                                                                               cgCATGGTATTGCGGGCGGTCATCGTCGGCTGGTTGGTGAGCCAGAAGGTTTTAGAACCCGCCGGTTTCATCATGTTCATCGGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCcgcg  <  1:208180/274‑1 (MQ=255)                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                             CTCCAGCTCTTCCGCGTTTAATCCCGGCGGAACATGATCGGTATTGCCATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCA‑GCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCAGCGACGGCTGCGTCAGATACAGATCCGGGTTCTTTTCATTGGCAAAGGTCAGCGCCTGTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCG  >  NC_000913/4148993‑4149520

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