Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,043	G→A	G169G (GGC→GGT)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,043	0	G	A	98.3%	184.8 / NA	59	G169G (GGC→GGT)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (36/21);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATT  >  NC_000913/3815801‑3816289                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCACCGTATTCCAGATCGCAAATTGCTTCACCGttta                                                                                                                                                                                                                                                 <  1:366916/251‑3 (MQ=255)       aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTAGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtac                                                                                                                                                                                                                                           >  1:235132/1‑249 (MQ=255)           cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                                                                                                                                                                       >  1:306623/1‑251 (MQ=255)              cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                                                                                                                                                                       >  1:438059/1‑248 (MQ=255)                            tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAcc                                                                                                                                                                                                                       <  1:243955/250‑1 (MQ=255)                                         gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCCGCGCCAACAAGGTCAGTCGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTCCCATCCCTTTCAttgg                                                                                                                                                                                                         >  1:288126/1‑248 (MQ=255)                                              cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCattggattgg                                                                                                                                                                                                    <  2:82711/251‑1 (MQ=255)                                                 cgccTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTAGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGttt                                                                                                                                                                                                 >  1:261527/1‑251 (MQ=255)                                                        tGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCCTTGGCTTGGTTTTCAGCtt                                                                                                                                                                                          >  2:360591/1‑251 (MQ=255)                                                             cGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGt                                                                                                                                                                                     >  2:86489/1‑251 (MQ=255)                                                             cGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGt                                                                                                                                                                                     >  2:184819/1‑251 (MQ=255)                                                                tACAATGGATTCGATTCCACC‑TCGGGCCAGCGCCACCAAGGCGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCACCGTATTCCAGATCGCAAATTCCTTCACCGTTTACAATACCGACAGAAACCGCGGCTCCCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGtttt                                                                                                                                                                                  >  2:67286/1‑251 (MQ=255)                                                                          cgagtccc‑cc‑TCGGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCGTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCGCCGTTTACAATTCCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTTCACCAGCtt                                                                                                                                                                        <  2:227980/244‑1 (MQ=255)                                                                             ttccc‑cc‑TCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt                                                                                                                                                                    >  2:396126/4‑251 (MQ=255)                                                                                    c‑TCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc                                                                                                                                                               >  1:396543/1‑250 (MQ=255)                                                                                         ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg                                                                                                                                                          <  1:122627/251‑2 (MQ=255)                                                                                            aagcgcaaacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAACTTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCCGCTTCTGTAGCGCATCtgca                                                                                                                                                        >  2:69555/8‑247 (MQ=255)                                                                                             agcgCCAACAAGGTGAGTAGCTCATCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCGGAGTCGTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACAGACAGAAACCGCGGCTACCATCCCTTTAATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTACTGTAgccgctctgcc                                                                                                                                                        >  2:372273/4‑240 (MQ=255)                                                                                               cgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCCGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGAATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcc                                                                                                                                                    >  1:83033/2‑250 (MQ=255)                                                                                                  caacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGcc                                                                                                                                                  >  1:254615/1‑250 (MQ=255)                                                                                                    acaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc                                                                                                                                                >  2:356388/1‑250 (MQ=255)                                                                                                          tGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAcc                                                                                                                                                                                                                       >  2:351291/1‑173 (MQ=255)                                                                                                          tGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAcc                                                                                                                                                                                                                       <  1:351291/173‑1 (MQ=255)                                                                                                           gAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAc                                                                                                                                        <  1:262726/251‑1 (MQ=255)                                                                                                             gTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                      <  2:306623/251‑1 (MQ=255)                                                                                                             gTAGCCCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCACGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                      >  2:222902/1‑251 (MQ=255)                                                                                                               aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGt                                                                                                                                    <  1:211593/251‑1 (MQ=255)                                                                                                                  tctTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                 >  1:136764/1‑251 (MQ=255)                                                                                                                   ctTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTaa                                                                                                                                  <  2:358316/249‑1 (MQ=255)                                                                                                                    ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCGACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCg                                                                                                                               >  1:164877/1‑251 (MQ=255)                                                                                                                    ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                 <  1:50134/249‑1 (MQ=255)                                                                                                                    ttCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCCCACAGGCACCCggcat                                                                                                                                 >  2:50134/1‑245 (MQ=255)                                                                                                                       aTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaag                                                                                                                            <  1:213632/251‑3 (MQ=255)                                                                                                                          ggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCgg                                                                                                                         >  2:354385/1‑251 (MQ=255)                                                                                                                                       ccTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACCATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaagc                                                                                                                           >  1:41527/1‑233 (MQ=255)                                                                                                                                       ccTTCTGCCGTCCCCGTCACTGCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCGGTAGCGCATCTGCCATCGCCACACAGGCACCCGTatgcgaagc                                                                                                                           <  2:41527/236‑9 (MQ=255)                                                                                                                                                          cTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCAAAGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCg                                                                                         >  2:207745/1‑251 (MQ=255)                                                                                                                                                                 gatgCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCAGCAGCCTGAAGCACGTCGCAGTCCa                                                                                  <  2:438059/251‑1 (MQ=255)                                                                                                                                                                   tgCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCTTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAg                                                                                 <  1:29982/250‑1 (MQ=255)                                                                                                                                                                      gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                              >  2:90046/1‑250 (MQ=255)                                                                                                                                                                      gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTa                                                                             <  1:180941/251‑1 (MQ=255)                                                                                                                                                                            cTTCGGTCATCCCTAAGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAACATCCCTTTCATTGGATTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCCCCATCAGCCCGACGCACGTCCCCGTCCAGCGTAATGGTg                                                                       >  1:359794/1‑251 (MQ=255)                                                                                                                                                                               cGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAc                                                                    <  1:32142/251‑1 (MQ=255)                                                                                                                                                                                     tcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAATAACGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                              >  1:431867/1‑251 (MQ=255)                                                                                                                                                                                       aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                              >  1:98349/1‑249 (MQ=255)                                                                                                                                                                                            gTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCTCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCtt                                                       <  2:39345/251‑1 (MQ=255)                                                                                                                                                                                                    cGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTTCCACCATCAGCCTGAAGCACGTCGCAGGCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATc                                                <  2:261527/250‑1 (MQ=255)                                                                                                                                                                                                            gCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACGCACCCAGCGCTTTCAAATCTACTGCcg                                        <  2:43264/250‑1 (MQ=255)                                                                                                                                                                                                                 tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTCATGGTGAACTCACCCAGCGCTTTCAAATCTCCTGCCGCGCGa                                   >  2:75994/2‑250 (MQ=255)                                                                                                                                                                                                                 tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                              <  2:41135/174‑1 (MQ=255)                                                                                                                                                                                                                 tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                              >  1:41135/2‑175 (MQ=255)                                                                                                                                                                                                                 tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAACGCTTTCAAATCTACTGCCGCGCGa                                   >  1:261153/2‑250 (MQ=255)                                                                                                                                                                                                                    gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAagtg                               >  2:85361/1‑249 (MQ=255)                                                                                                                                                                                                                           aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGCAGTGCACGGGc                        >  2:200923/1‑251 (MQ=255)                                                                                                                                                                                                                              tATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                     >  1:243691/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCCGCGGCTACCATCCCTTTCATTGGATTGGTTGTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGatg               <  2:302192/251‑3 (MQ=255)                                                                                                                                                                                                                                      aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt             <  1:396126/251‑1 (MQ=255)                                                                                                                                                                                                                                                gCTTCACCGTTTACAATACAGACAGAAACCCCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCAGTTTTCCAACAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAt   >  1:159969/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAtt  >  1:151756/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATT  >  NC_000913/3815801‑3816289

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