Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,100,529	C→T	intergenic (‑4/‑281)	rhaT ← / → sodA	rhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,100,529	0	C	T	97.9%	159.0 / NA	48	intergenic (‑4/‑281)	rhaT/sodA	rhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/29);  minor base G (0/1);  total (18/30)
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.96e-01

GCCGAACAGAAAAACAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCG  >  NC_000913/4100286‑4100770                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  gcCGAACAGAAAAACAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATAACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCt                                                                                                                                                                                                                                            <  2:469903/251‑1 (MQ=255)        aGAAAAACAGGCAGTCGCGTATAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCAtt                                                                                                                                                                                                                                     <  2:610995/251‑1 (MQ=255)           aaaaCAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCAttttt                                                                                                                                                                                                                                  <  1:164237/251‑1 (MQ=255)                    ggtagcgTAGAGACACGAAACGAGCCTTAATACGCCCAGAAATTAGGTAGTAACAGGGCGCTGATGGCCCAAGGCAGAATAACTCACGAAACAAGCCCACCGACTGACCAAATGGTTTCCCATGACCATTTTTTTACTGTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAAGATCCCCATCGTAATCGCGTTACTCATCTTTTGAGCCTCATCATTTTTCGTCGCGTc                                                                                                                                                                                                                         <  2:574736/247‑1 (MQ=255)                     gTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTc                                                                                                                                                                                                                         >  1:242361/1‑250 (MQ=255)                            agagagACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTcc                                                                                                                                                                                                                 <  1:346905/251‑1 (MQ=255)                                                        cAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGt                                                                                                                                                                                     <  1:527525/251‑1 (MQ=255)                                                         aGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTc                                                                                                                                                                                    <  2:414159/251‑1 (MQ=255)                                                          gAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCcgac                                                                                                                                                                                                              >  2:396616/1‑224 (MQ=255)                                                          gAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCcgac                                                                                                                                                                                                              <  1:396616/224‑1 (MQ=255)                                                                  gtagtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCt                                                                                                                                                                           >  2:438150/1‑251 (MQ=255)                                                                    agtaACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCttt                                                                                                                                                                         <  1:314617/251‑1 (MQ=255)                                                                       aaCAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGc                                                                                                                                                                      <  1:650063/251‑1 (MQ=255)                                                                         cAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCt                                                                                                                                                                    >  2:35319/1‑251 (MQ=255)                                                                                    tGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGt                                                                                                                                                         >  2:631138/1‑251 (MQ=255)                                                                                           cGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGAt                                                                                                                                                  >  1:530727/1‑251 (MQ=255)                                                                                           cGGCAGAATAATCCACGAAACAAGCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTATCCTTTCGCCTGTTGCCGCCGTTGTCGAt                                                                                                                                                  <  1:538907/251‑1 (MQ=255)                                                                                                         aCGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCAc                                                                                                                                    <  1:122215/251‑1 (MQ=255)                                                                                                                     cACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTgg                                                                                                                        >  1:607589/1‑251 (MQ=255)                                                                                                                        cGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTg                                                                                                                     >  2:46139/1‑251 (MQ=255)                                                                                                                          ccTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTACCCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGAt                                                                                                                   <  1:148282/250‑1 (MQ=255)                                                                                                                             gACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTg                                                                                                                >  1:60502/1‑251 (MQ=255)                                                                                                                                 acaTGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTc                                                                                                            <  1:438150/251‑1 (MQ=255)                                                                                                                                 acaTGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTAACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATAAAGGCATTAAGTGGGGGATTTGCTTc                                                                                                            >  1:13520/1‑251 (MQ=255)                                                                                                                                   aTGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTcac                                                                                                          >  2:522201/1‑251 (MQ=255)                                                                                                                                   aTGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTcac                                                                                                          >  2:252950/1‑251 (MQ=255)                                                                                                                                           ccANGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCGCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTTGGTGATTTGCTTTACATCTCggg                                                                                                  <  2:607589/251‑1 (MQ=255)                                                                                                                                               gACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCAt                                                                                               <  1:252950/250‑1 (MQ=255)                                                                                                                                                              tttttGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCa                                                                                <  1:464361/250‑1 (MQ=255)                                                                                                                                                                     aCGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCtt                                                                         >  2:273883/1‑250 (MQ=255)                                                                                                                                                                       ggAGCGTAAAAACAGGCTGCACTGTCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTGTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCATGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACg                                                                      <  2:530727/251‑1 (MQ=255)                                                                                                                                                                                     ggCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCAtt                                                        >  2:496548/1‑251 (MQ=255)                                                                                                                                                                                         gCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGAt                                                     <  1:522201/250‑1 (MQ=255)                                                                                                                                                                                          cacgggCCGCGCCGATCACATGCCAAACTATCCCCATCGTACTCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCCCATCTCCGACGAGATGAGGGTAAAAATCGTTCTTTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGTCATTTTCCTGCAACACCATACCCTTACGAAAAGTACGGCATTGATaa                                                   <  2:355815/247‑1 (MQ=255)                                                                                                                                                                                              ggCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATc                                                 <  2:60502/249‑1 (MQ=255)                                                                                                                                                                                                  gcgcCGATCAAATGCCAAAATGTCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGCTGAGTGTCAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCa                                                                                                                                     >  2:586525/1‑161 (MQ=255)                                                                                                                                                                                                  gcgcCGATCAAATGCCAAAATATCCCCAGCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTTTCGATTAGCCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCa                                                                                                                                     <  1:586525/161‑1 (MQ=255)                                                                                                                                                                                                       gATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCCTCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAatat                                      >  1:210554/1‑251 (MQ=255)                                                                                                                                                                                                            aaTGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCAttt                                 <  1:264160/251‑1 (MQ=255)                                                                                                                                                                                                               tccAAAATATCCCCATCGTAATCGCGTTGCTCATCTTTTTGTCCTCATCATTTTTCGGCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTTAGTGGGTGATTTGCGTCACATCTCGGGCATTTTCCTGTATAAACATGCCCTTACAAAAAGTACGGCAGTGATAATCATTTACAATACCATttaa                               <  2:263322/249‑1 (MQ=255)                                                                                                                                                                                                               gCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATttaa                               <  1:496548/250‑1 (MQ=255)                                                                                                                                                                                                                ccAAAATATCCCCATCGTAATCGCGTTACTCATCTTTTTATCCACATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCGTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATttaatt                             <  2:412422/251‑1 (MQ=255)                                                                                                                                                                                                                   aaaTATCACAATCGTAATCACGTTACTCATCTTTTTATCATCATCATTTTTAGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATAGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGAATTAAGTGGGTGATTTGCTTCACATCTCGAGCATTTTCATGCAAAACCATACCCTTACGAAAAGTACGGAATTGATAATCATTTTCAATATCATTTAATTAAc                          >  2:274326/1‑251 (MQ=255)                                                                                                                                                                                                                     atatCCCCTGCGTAATCGCGTTAATCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACAGCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAATCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATGTAAGTAACta                        <  2:342434/251‑1 (MQ=255)                                                                                                                                                                                                                               cGTAATCGCGTTAATCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACTAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCa              <  2:325117/251‑1 (MQ=255)                                                                                                                                                                                                                                     cgcgTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCt        <  1:37449/251‑1 (MQ=255)                                                                                                                                                                                                                                       cgTTACTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTa      <  1:426764/251‑1 (MQ=255)                                                                                                                                                                                                                                           aCTCATCTTTTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGGGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTAAGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTAcgcg  >  1:385126/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  GCCGAACAGAAAAACAGGCAGTCGCGTAGAGAGACTAAACGAGCTGTAATACGCCCAGAAATTCGGTAGTAACAGGGCGCTGATGGCCCACGGCAGAATAATCCACGAAACAATCCCACCGACTGACCACATGGTTTCCCATGACCATTTTTTTACTTTTTTGAACGGAGCGTAAAAACAGGCTGCACTGGCCGCGCCGATCAAATGCCAAAATATCCCCATCGTAATCGCGTTACTCATCTTCTTATCCTCATCATTTTTCGTCGCGTCACATCTCCGACGAGATGAGTGTAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCG  >  NC_000913/4100286‑4100770

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