Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,810,758	0	G	A	31.2%	56.7 / 37.9	48	intergenic (+79/‑12)	ygaH/mprA	L‑valine exporter, YgaH component/DNA‑binding transcriptional regulator MprA
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base A (11/4);  total (32/16)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTACCGCACCAGAAGTGATGCACGATACACGCCGTTTCGTGCCCACGCTGGTCGGCTTCGCGGTACTGGGTGCCAGTTTCTATAAAACACGCAGCATTATCATCCCAACACTGCTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTG  >  NC_000913/2810514‑2811007                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          ggcaatttgcccgcatcaatattcatgccacatttgccatcaggggttgcctcagatgctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatatTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATa                                                                                                                                                                                                                                                     <  2:33150‑M2/135‑1 (MQ=255)        tgctcgcatcaatattcatgccacatttgccatcaggggttgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAAt                                                                                                                                                                                                                                              >  2:69056‑M2/110‑251 (MQ=255)                tcaatattcatgccacatttgccatcaggggttgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAACGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGtt                                                                                                                                                                                                                                      >  1:692961‑M2/102‑251 (MQ=255)                        catgccacatttgccatcaggggttgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                                                                                                                                                                               <  1:424080‑M2/157‑1 (MQ=255)                                     ccatcaggagttgcctcagattatcattatgttagggtagaaaaaagcgactatttccattgggtaatctatcgacatagTTAGTACGCTGGCCTATGGGCTCGACTGGAAAGTGATGGCGATTATATAACACACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATAACCAAATGGATAGTTCGTTTATGCCCATTGAACaa                                                                                                                                                                                                                  >  2:320743‑M2/81‑250 (MQ=255)                                          aggggttgcctcagattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGGTGGCCTATGGGCTCGCCTGGAAAGTGATGACGATTATATAACCCACAAGGATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCAAACCCAAATGGTTAGTTCGTTTACGCCCCTTGAACAACTGCt                                                                                                                                                                                                             >  2:165537‑M2/76‑250 (MQ=255)                                            gggttgcctccgattctcagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAGCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaa                                                                                                                                                                                                           >  2:764011‑M2/74‑250 (MQ=255)                                                             cagtatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCAACAAGAATCATTTTTCTAAAAAAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACTCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTATCGCCAgcagcc                                                                                                                                                                                         >  2:439212‑M2/57‑251 (MQ=255)                                                                tatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAACACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACg                                                                                                                                                                                      >  2:734248‑M2/54‑251 (MQ=255)                                                                tatgttagggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACg                                                                                                                                                                                      >  1:71244‑M2/54‑251 (MQ=255)                                                                      agggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                                                                                                                                                                                >  1:625453‑M2/48‑251 (MQ=255)                                                                      agggtagaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAAAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                                                                                                                                                                                >  2:505538‑M2/48‑251 (MQ=255)                                                                       gggtagaaaaaattgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                                                                                                                                                                               <  2:810743‑M2/205‑1 (MQ=255)                                                                            gaaaaaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                                                                                                                                                                          >  1:617707‑M2/42‑251 (MQ=255)                                                                                aaagtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                                                                                                                                                                       >  1:475762‑M2/38‑250 (MQ=255)                                                                                   gtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAGCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga                                                                                                                                                                   <  1:764011‑M2/217‑1 (MQ=255)                                                                                   gtgactatttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg                                                                                                                                                                    >  1:875568‑M2/35‑250 (MQ=255)                                                                                         atttccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt                                                                                                                                                             <  2:82412‑M2/223‑1 (MQ=255)                                                                                             ccattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg                                                                                                                                                         <  2:381887‑M2/227‑1 (MQ=255)                                                                                              cattgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGa                                                                                                                                                        <  1:379592‑M2/228‑1 (MQ=255)                                                                                                 tgggtaatatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTc                                                                                                                                                     >  1:171069‑M2/21‑251 (MQ=255)                                                                                                       atatatcgacatagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCttt                                                                                                                                               <  1:770666‑M2/237‑1 (MQ=255)                                                                                                                 ctagTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcctgcacat                                                                                                                                     >  2:102057‑M2/5‑251 (MQ=255)                                                                                                                    gTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgc                                                                                                                                   >  2:4632‑M2/2‑250 (MQ=255)                                                                                                                                gCCTATGGGCTCGCCGGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCGGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTGATCAGGAGATCCTGCTGACTCGTCTTTGCATGCACATGCAAAGCAAgctgct                                                                                                                      <  2:623355/251‑1 (MQ=255)                                                                                                                                                    aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                                                                  >  1:302697/1‑251 (MQ=255)                                                                                                                                                          ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAgg                                                                                            >  2:130633/1‑251 (MQ=255)                                                                                                                                                              aTTATATAACCCTCAAGAATCATTTTTCTAAAACAATTCATTTACTTTAGTTGTCACTGTCGTTACTATATCGGGTGAAATAAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTGATCAGGAGAGCCTTCTGACACGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCTCAATAAAATGCTGAAGGCTCa                                                                                        <  2:692961/251‑1 (MQ=255)                                                                                                                                                               ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTCTATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTc                                                                                         >  2:617802/1‑249 (MQ=255)                                                                                                                                                               ttATATAAACCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTc                                                                                         >  1:632659/1‑249 (MQ=255)                                                                                                                                                                          acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACg                                                                             >  1:784372/1‑250 (MQ=255)                                                                                                                                                                                   aTTTTTCTAAAACAATACATTTACTTTATTTGTCATTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                                   <  1:170786/251‑1 (MQ=255)                                                                                                                                                                                        tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATg                                                               <  2:617707/250‑1 (MQ=255)                                                                                                                                                                                         cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGTCTCAGGGGATTAACGAGACGTTGTTTATGGc                                                             >  2:706870/1‑251 (MQ=255)                                                                                                                                                                                            aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                                          >  2:42855/1‑251 (MQ=255)                                                                                                                                                                                             aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                                                         >  2:139104/1‑251 (MQ=255)                                                                                                                                                                                              ccAATACCTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAATTCATACCCAAATGGATAGTTCGTTTACGCCCATTTAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGCTTTTCCTTATCAGGAAATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAAAGCTGAAGGCTCAGGGGATTAACGCGCCGTTGTTTATGTCGTTg                                                         >  2:564499/2‑250 (MQ=255)                                                                                                                                                                                                 aTACATTTACTTTATTTGTCACTGACGTACCTATATCGGCTGAAATTAATGAAGTCATACCCAAATGTATAGCTCTATTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATTCAAAGCAAGCTGCTGGAGAACCGCAATCAAATGCTGAACGCTCAGGCGATTAACGATGCGGTGTTTATGGCGTTGAttc                                                     >  2:626436/1‑250 (MQ=255)                                                                                                                                                                                                     aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt                                                 <  2:625453/251‑1 (MQ=255)                                                                                                                                                                                                       ttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                                               >  1:65040/1‑251 (MQ=255)                                                                                                                                                                                                            ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtct                                          >  2:182050/1‑251 (MQ=255)                                                                                                                                                                                                             tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtct                                          <  1:69056/250‑1 (MQ=255)                                                                                                                                                                                                                  gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTGCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACGCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                                    <  1:648845/251‑1 (MQ=255)                                                                                                                                                                                                                       tGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCCGGAAAAcc                                >  2:628714/1‑250 (MQ=255)                                                                                                                                                                                                                          cGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACcaca                             <  1:439212/250‑1 (MQ=255)                                                                                                                                                                                                                            ttACTATATCGGCTGAAATTCATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGTCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGt                           >  2:703481/1‑250 (MQ=255)                                                                                                                                                                                                                                      ggCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTc                >  1:832201/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgtg  >  2:827615/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          CTCTACCGCACCAGAAGTGATGCACGATACACGCCGTTTCGTGCCCACGCTGGTCGGCTTCGCGGTACTGGGTGCCAGTTTCTATAAAACACGCAGCATTATCATCCCAACACTGCTTAGTGCGCTGGCCTATGGGCTCGCCTGGAAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTG  >  NC_000913/2810514‑2811007

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