Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	15,226	C→T	G353G (GGC→GGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,226	0	C	T	94.6%	100.1 / ‑4.0	39	G353G (GGC→GGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (13/22);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAA  >  NC_000913/14986‑15464                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                ggcggcGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCacaaca                                                                                                                                                                                                                                       <  1:229422/251‑1 (MQ=255)            gTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCa                                                                                                                                                                                                                            <  1:321470/251‑1 (MQ=255)            gTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCa                                                                                                                                                                                                                            <  1:182578/251‑1 (MQ=255)                   aaTCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGct                                                                                                                                                                                                                    <  1:315207/251‑1 (MQ=255)                      cGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGGCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGcttc                                                                                                                                                                                                                  <  1:221192/250‑1 (MQ=255)                      cGAAGTACCGAACCTTGATGGTCGCGGCAAACTGACAGTGCCTGGCGAAACCCAGACCGGTACGCTTTTCCGTATGCGCGGTACAGGCGTCAAGTCTGTCCGCGGTGCCCCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGCAGCGGAAACCGCTGCTGCAAGAGCCGCACGAAAGCTTCGGGGGCCCAACCGGAGAGCACCACACCCCGCGCTCC‑Aagacatact                                                                                                                                                                                                                 >  2:350922/1‑245 (MQ=255)                            aCCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGAt                                                                                                                                                                                                            <  1:239708/250‑1 (MQ=255)                             ccggcgcTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTACGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGTCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGATGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCAAAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATgg                                                                                                                                                                                                          <  1:98661/245‑1 (MQ=255)                                   cTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCAGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGgtgaga                                                                                                                                                                                                     >  1:239972/1‑247 (MQ=255)                                           tCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGGGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGTTCA‑AAAAACTTATTTGATGGGGGgaagaacgatc                                                                                                                                                                                            >  2:91296/1‑246 (MQ=255)                                           tCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTgaagaa                                                                                                                                                                                                 <  1:91296/246‑1 (MQ=255)                                            cgcttCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCCCAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAACCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGGGGCCCAACCGGTGAGCACAAACGCCCGCGCTCA‑CAGAGCTTCTTTGATGGTGTGGATAAGttttt                                                                                                                                                                                           >  2:121477/5‑251 (MQ=255)                                            cgcgTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGtttt                                                                                                                                                                                            <  1:251252/250‑1 (MQ=255)                                                     cTGAAAGTGCCTGGAGAAACCCAGACCGGTAAGCTTTTCCGGATGCGCGGCAAAGGCGCCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTGGCAGGGCCGCGTTGACGTCGGAACACCGGTCGACGTGAACGAGAACAAGAGAGAGCTGCTGCAAGAGATGCAAGAACGCTTCGGGGGGCCAACAGGAGAGCAAACCCGCACGCGC‑CACAAGACCATCTTTGATGGGGAGAAGACGATTTTTGACGACCt                                                                                                                                                                                  >  2:239963/1‑251 (MQ=255)                                                                 ggCGAAACCCAGACCGGTAAGCTGTTCCGTATGCGCGGTAAAGGCGGCAAGTCTTACCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGGCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAAcc                                                                                                                                                                      <  1:259046/251‑1 (MQ=255)                                                                             aCCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGcc                                                                                                                                                          <  1:218698/251‑1 (MQ=255)                                                                                     gcttttCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTgaa                                                                                                                                                                                                    >  1:110665/5‑201 (MQ=255)                                                                                               aTGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGcttg                                                                                                                                        >  1:11322/1‑248 (MQ=255)                                                                                                 gcgcGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTggg                                                                                                                                      <  1:56903/251‑1 (MQ=255)                                                                                                           ggcggcaggTCTGTCCGGGGTGGCGCTCAGGGTGATTTGCTGTGCCGGGGTGTCGTCGAAACTCCGGTAGGCCCGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGAAGACCTGACCCGCTAACCTCCCTAAAAGCCGGCCCGCGGGCAGGCCTGGGTAAAAATAgg                                                                                                                            <  2:217816/243‑1 (MQ=255)                                                                                                                 aaGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGCGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAAAAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGAAAAATTAGGGTgagt                                                                                                                      >  1:226369/1‑248 (MQ=255)                                                                                                                     cTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGa                                                                                                                   >  1:314710/1‑250 (MQ=255)                                                                                                                       atCCGCGGAGGCGCACAGGGAGATTTGATGAGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGCGACGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCACCCGGTGAGCACAACAACCAGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAACAGATATTTGACGCCCAGACCCGCTAACCTACCCAAAAGCCGGCCCGTGGGAAGGCTTGGGTaaccttaaagcgcattgaggt                                                                                                                >  2:263484/2‑232 (MQ=255)                                                                                                                               gggCGCACGGGGTGATTTGGTAGGACGCGTTGTCGACAAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGTCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTTCGTTGAAGATATGCGAg                                                                                                        <  2:166403/250‑1 (MQ=255)                                                                                                                                  cgcACAGGGTGATTTGCTGTTCCGCGTTGTCGTCGAAACACCGCTTGGCCCTAACGAGAAGCAGAAACAGCTGCTTCAAGCGCTGCAAGAAAGCTTCGGTGGCCCAATCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGACGGGGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCATTGAAGATATGCGAGCAc                                                                                                     <  2:126426/251‑1 (MQ=255)                                                                                                                                     acaGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTg                                                                                                  <  1:32788/251‑1 (MQ=255)                                                                                                                                                      gcctcgTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACTCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCa                                                                                 <  1:141880/247‑1 (MQ=255)                                                                                                                                                        cgagTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGCAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCTTGAGCAGGCTTGGGAAAAATTAGGGCACGTTGAAGATATTCGAGCACATGAAAAGGGGCggggtgacac                                                                               >  2:21931/4‑245 (MQ=255)                                                                                                                                                                cgtcgAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggagacacccctaa                                                                        >  1:382962/1‑237 (MQ=255)                                                                                                                                                                cgtcgAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCCGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCTCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGATCACTCCCATAAg                                                                       <  1:190544/251‑1 (MQ=255)                                                                                                                                                                         acCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggacaatcccatcaacgcaaactt                                                              >  1:351540/1‑228 (MQ=255)                                                                                                                                                                         acCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACtt                                                              <  1:315028/251‑1 (MQ=255)                                                                                                                                                                           cGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCAAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGaatacacccagacaaacgacatta                                                             >  1:202854/1‑227 (MQ=255)                                                                                                                                                                             gTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCACGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGACGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGATGCGTTGAAGATATGCGAGCAACTGAAAAGTGGCgggggatactactacaagagcaacaaga                                                             >  1:335145/1‑224 (MQ=255)                                                                                                                                                                                          agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                             <  1:316109/247‑1 (MQ=255)                                                                                                                                                                                                aGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                       <  1:148710/251‑1 (MQ=255)                                                                                                                                                                                                       gctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                                <  1:350922/251‑1 (MQ=255)                                                                                                                                                                                                          gctgcaagagctgcaagaAAGCTTCGGTGGCCCACCCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                             <  1:69142/251‑1 (MQ=255)                                                                                                                                                                                                                                     gggCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTACCGATGAATTACTCTCACGATaa  <  1:121477/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                GGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAA  >  NC_000913/14986‑15464

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