Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,764	T→A	A138A (GCT→GCA)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,764	0	T	A	96.5%	241.6 / 1.7	87	A138A (GCT→GCA)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (39/43);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

GCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGG  >  NC_000913/71521‑72007                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      gCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttca                                                                                                                                                                                                                                                <  2:876822/251‑3 (MQ=255)   aCGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATtcatt                                                                                                                                                                                                                                              >  1:799652/1‑251 (MQ=255)     gCCTGGCTGGCAGGGCTTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGATCTGCCGGGGGCAAAATTCATTAc                                                                                                                                                                                                                                            >  1:735934/1‑251 (MQ=255)       ctggctggCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTGCCGGGGGCAAAATTCATTACGc                                                                                                                                                                                                                                          >  2:663762/1‑251 (MQ=255)         ggctggcAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTAGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGACCGACGCGTCCGCTGGTGCCAATGGCGGCGGGAATGCTGGGTCTGCCGGTGGCAAAAttcaatccgccc                                                                                                                                                                                                                                        >  1:919264/1‑241 (MQ=255)            tggcagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCACAATTCATTACGCCGAat                                                                                                                                                                                                                                     >  1:830481/25‑251 (MQ=255)             ggcagggattgtcggttgttagctGGGCGACTGGATCTCCTTATGGCTGGGATGGCGGTTTAAAAAGCCGTTGCATCGCTGGTCATGTCTGAAGAAAAACAAATCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGGCGGGCGCTTTGTGGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAata                                                                                                                                                                                                                                    <  2:216042/228‑1 (MQ=255)             ggcagggattgtcgggtgtctgctGGGCGACTGGATCTCCTTCGGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAATCATTACGCGATAAAACTGAACAGGCGTTGCATCAACACAGCATGTGCGCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAata                                                                                                                                                                                                                                    <  2:47003/228‑1 (MQ=255)               cagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTTGATCTGCCGGTGGCAAAATTCATTACGCCGAAtatt                                                                                                                                                                                                                                  <  1:542979/230‑1 (MQ=255)                       gtcggttgtctgctGGGCGACTGGATCTCCTTCTGGATGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGct                                                                                                                                                                                                                           <  1:276296/237‑1 (MQ=255)                         cggttgtctgctggtcGACTGGATCTCCTTCGGGCGTGGCTGGCGTTTTAAAA‑GCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGGTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGctgcct                                                                                                                                                                                                                       <  2:479180/236‑1 (MQ=255)                                  gctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCtgt                                                                                                                                                                                                                  <  1:898311/246‑1 (MQ=255)                                    ggggcgatgggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATGTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTGGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcc                                                                                                                                                                                                              <  2:315832/233‑1 (MQ=255)                                     gggcgaccggtgttccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGTGGGTCATTTCGGAAGAAAACCAAAGCATTACTCGATAAAACTGAACATTCGTTGCAGCAACACAGCATGTTCACCTTTCTGGTCGGTCTCTTTGTTGGCCCGACGCGTCCGCTGGTTCAAATGGTGGCGGGAATGCTGGTACTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcc                                                                                                                                                                                                              <  2:883500/233‑1 (MQ=255)                                        cGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCCACCGt                                                                                                                                                                                                         >  1:592820/1‑251 (MQ=255)                                         gattGGATCTCCTTCTGGCTGGGCTGGTTGTTTAAAATGCCGTTGCATCGCTGGGCATTTCTTAAGAAAAACAAAGCATTCCCCGATAAAACTGAACATGCTTTGCATCAACAAAGCATGTTCACCATTCTGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAGGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcaaccgtt                                                                                                                                                                                                        <  2:657920/248‑8 (MQ=255)                                             ggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTAc                                                                                                                                                                                                    >  1:285485/9‑251 (MQ=255)                                               atctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACtt                                                                                                                                                                                                  >  1:937636/7‑251 (MQ=255)                                                   ccttctgtgtggGGTGGCTTTTTAAAAAGCCGTTGCAACGCTGGTAATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTGGGCCCGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTGTACTTCCTg                                                                                                                                                                                              <  2:552719/242‑1 (MQ=255)                                                   ccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTg                                                                                                                                                                                              >  1:677797/3‑251 (MQ=255)                                                     ttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGGATGCTGGGTCTGCCGGGGGCAAAATTCCTTACGCCGATTATTTTCGGCTGCCTTATGGGGCCCCCGTTTTACTTTCTGcc                                                                                                                                                                                            >  2:73009/1‑251 (MQ=255)                                                      tctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGc                                                                                                                                                                                             <  1:210389/249‑1 (MQ=255)                                                              ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAATCAATACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGTCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                   <  1:398158/248‑1 (MQ=255)                                                              ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                   <  1:971861/248‑1 (MQ=255)                                                               gTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTg                                                                                                                                                                                  <  1:308489/251‑1 (MQ=255)                                                                gcaaccatttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCACCATTCATTACACCGAATATTATCGGCTGCCCGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTgg                                                                                                                                                                                 >  2:505415/8‑251 (MQ=255)                                                                      ttttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggg‑cg                                                                                                                                                                          >  1:336791/1‑251 (MQ=255)                                                                       tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggg‑cg                                                                                                                                                                          <  1:691650/250‑1 (MQ=255)                                                                       tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGg‑cgc                                                                                                                                                                         >  1:715484/1‑251 (MQ=255)                                                                            aaaGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGAGTCTGGATCTGCCGGTGGCAAAATTCATTACCCCGAATATTATCGGCCGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGagtcgtgcgggggcgggcg                                                                                                                                                                     >  2:356411/1‑233 (MQ=255)                                                                             aaGCCGTTGCATCGCTGGTCATTTCTGTAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGcgat                                                                                                                                                                   <  2:677797/251‑1 (MQ=255)                                                                               ggcGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCTGG‑CGCGGcgatcg                                                                                                                                                                 <  1:585811/249‑1 (MQ=255)                                                                                       atcccTGGTCATTTTTTAAGAAAAACAAAGCATTACCCGAGAAAACTGAACATGCGTTGCACTACCACAGCATGTTCACCATCCGGGTCGGTCGCTGTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCTCCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑Tcca                                                                                                                                                        <  2:592820/247‑2 (MQ=255)                                                                                         cGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGGACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGACCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGCTGGCAAAATTCCTTACGCCGAATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGGGGG‑CGCGGGGATCGATAt‑ttccgc                                                                                                                                                      >  2:132311/1‑246 (MQ=255)                                                                                          acTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGACGCTTTGTTGGCCCGACGCGTCCGCTGGTGCAAATGGTGGCGGGAATGCTGGATATGACGGTAGCAAAATTCATTACGCCGCATATTATAGGCTGCCTTCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggggctgggggtcgatcattcagcc                                                                                                                                                       >  2:508318/2‑229 (MQ=255)                                                                                             gggcATTTCGGAAGAAAAAAAAAGCATTACTCGATAAAACGGAACAGGCGTTGCATCAACACAGCATGGTCACCATTCGGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCg                                                                                                                                                    <  2:285485/246‑1 (MQ=255)                                                                                             ggTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCCGCCggg                                                                                                                                                  >  1:114978/1‑250 (MQ=255)                                                                                                cATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGTCGATAT‑TCCAGCCGGAATg                                                                                                                                               >  1:232588/1‑251 (MQ=255)                                                                                                 aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATGTTGGGTCTGCCGGTGGCCACATTCGTTACGCCGGATATTATCGGCGTCCTGCTGTGGCCACCGTTTTAACTCCTTCCTGGGCATCTGGCGGG‑CGGGGcgagcgggttttcagccgggaagg                                                                                                                                               >  2:16574/1‑230 (MQ=255)                                                                                                 aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGGCGATAT‑TCCAACCGGAATgg                                                                                                                                              >  1:648349/1‑250 (MQ=255)                                                                                                   ttCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTTCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCAGCCGGGATGc                                                                                                                                              >  1:635148/1‑249 (MQ=255)                                                                                                    cccgccgcccttccgatctttactCGATAAAACTGAACATGCGTTGCATCAACCCAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCaga                                                                                                                                           <  2:503087/231‑1 (MQ=255)                                                                                                        aagaaAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCgg                                                                                                                                       <  1:572563/251‑1 (MQ=255)                                                                                                          gaaAAACAAAGCATTACGCGATAAAACTTAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGGCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCggcg                                                                                                                                     <  1:132311/251‑3 (MQ=255)                                                                                                           aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGG‑CGCGGCGCTCGCTAT‑TTCAGCCGGAATGCAGAGCggcgg                                                                                                                                    >  1:674943/1‑248 (MQ=255)                                                                                                           aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCAGCCGGAATGCAGAGCggggg                                                                                                                                    >  1:891673/1‑248 (MQ=255)                                                                                                            aaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGCCGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGAGTCTGGCGGG‑CGCGGCGGTCGAATT‑TTCAGCCGGAAATCAGAGCggcggg                                                                                                                                   >  1:97749/1‑247 (MQ=255)                                                                                                                cAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCggcgcgttt                                                                                                                                >  1:983021/1‑243 (MQ=255)                                                                                                                     cattactCGATAAAACTGAACATGCGATTCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTACGCTGGTGCCAATGGTGGCGGGAATGCCGGAACTGCCGGTGGCAAAATTCATTACGCCGAATAATATCGGATGCCCTCTATGGCCACCGTTTTAATTCCTACCTGGGATTCTGGCGGG‑CGCGGAGATCGA‑ATATCCAGCCGGAATGAAGAGCgggaggattaaatgg                                                                                                                          >  1:917214/4‑238 (MQ=255)                                                                                                                      attactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGCATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑GGGGGGGGTCGATAT‑TCCCGCCCGAATgacgaggggcgaggtttaatggt                                                                                                                         >  2:307788/3‑229 (MQ=255)                                                                                                                       ttactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCaga                                                                                                                                           >  1:503087/2‑232 (MQ=255)                                                                                                                        tacccgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGtt                                                                                                                        <  1:920128/245‑1 (MQ=255)                                                                                                                         actcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACACCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTtgc                                                                                                                      <  1:505415/247‑1 (MQ=255)                                                                                                                           tcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCCGAGCGGCGGGTTTAAATgggttctg                                                                                                                    >  1:677428/3‑245 (MQ=255)                                                                                                                              aTAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                 <  2:336791/251‑1 (MQ=255)                                                                                                                              aTAAAACTGAACATGCGTTGCATCAACACAGCATTTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                 <  1:191263/251‑1 (MQ=255)                                                                                                                               aaaaactaaaccggcgctGCATCCACACTGCAGGTGCACCATTCGGGGCGGGCGCTTTGTTGGCCCGACGCTGCCGCGGGTGCCAATGTGGGCGGGAAGGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                                <  2:388327/234‑1 (MQ=255)                                                                                                                                  aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCAGAGCGGGGGGTTTAAATGGTTGCTGCTGGcgg                                                                                                             >  1:708062/1‑249 (MQ=255)                                                                                                                                      gAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATg                                                                                                                                               <  1:630070/213‑1 (MQ=255)                                                                                                                                       aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGGGGGCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGAAGAGCGGCGGGGTTAAATGGTTGGTGCTGGcggaggcg                                                                                                        >  1:142919/1‑244 (MQ=255)                                                                                                                                          aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGATGGGTATGCCGGTGGCAAAATTAATTAAGCCGAATATTATCGGCTGACTGGTGTGGACACCGTTTTAATTCCCGCCTGGGGTTCTGGGGGGGCGCGG‑GGTCGATAT‑TCCAGGCGGCATCCAGAGCGGGGGGTTTTAATGGTtgctggcggggacggaggc                                                                                                      >  1:415293/1‑236 (MQ=255)                                                                                                                                            gctttGCATCAACCCAGCAGTTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtat                                                                                                    <  2:114978/247‑3 (MQ=255)                                                                                                                                            gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGAGAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtatt                                                                                                   <  1:637119/251‑4 (MQ=255)                                                                                                                                                tGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGCATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATTT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTTAATGGTTGCTGCTggggacggcggtatttgtg                                                                                               >  1:81655/1‑234 (MQ=255)                                                                                                                                                           aGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGGATGCAAAGCGGCGAGTTTAAATGGTTTCTTACGGCGACGGCGGTATTTGTGTGGGTttgtgg                                                                                    >  1:611649/1‑251 (MQ=255)                                                                                                                                                                 ttCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGGTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTTGTggctggct                                                                              >  1:269377/1‑251 (MQ=255)                                                                                                                                                                         tCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑Tccag                                                                                                                                                       <  1:651795/170‑3 (MQ=255)                                                                                                                                                                             gacggtcgCTTTGTGGGCCCGCCGCGTCCGCGGGTGCCAATGGTGGCGGGAATGGTGGATCTGCCGGTGGCAAAATTCATTGCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCGGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGtt                                                                  <  2:203310/249‑1 (MQ=255)                                                                                                                                                                               cggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGTCGATAT‑TCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGGCGGCGGGATTTGTGTGGGGTGGTGGCTGGCTGTGCTGGCgggttt                                                                >  1:603635/1‑247 (MQ=255)                                                                                                                                                                                ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTGGGGGGTGGCTGTGCTGGCgggtgtt                                                               >  1:230587/1‑246 (MQ=255)                                                                                                                                                                                      tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGGATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTTGGTGGGGTGGGGGCTGGCTGtggtgggggggtttgggcgg                                                         >  1:117252/1‑237 (MQ=255)                                                                                                                                                                                      tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGCGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTTGTGGCTGGCTGTGCTGGCgggtttggggggg                                                         >  1:92557/1‑240 (MQ=255)                                                                                                                                                                                        tggtgtcccggcgggtcctcgggTGCAAATGGTTGCGGGATTGTTGGATCTGCCGGTGGCAAAATGCATTACGCCGAATATTATCGTCTGCCGGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCTGCGATTGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGTTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                        <  2:64203/229‑1 (MQ=255)                                                                                                                                                                                        ggtttgcccaccgcgTCCGCTGGTGCCAATGGTGGCGGGAAGGCTGGCTCTGCCGGTGTCCAAATTCACTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGGTTGACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGGTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                        <  2:394733/239‑1 (MQ=255)                                                                                                                                                                                             gCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TTCCGCCGGGATGCAGAGCGGCGAGGTTAAATGGGTGCTGGTGGCGGCGGCGGTTTTTTTGTGGGGTGGGGGCTGGCGGTGCTTGCgggttttgtggtccgggg                                                    >  2:307855/1‑233 (MQ=255)                                                                                                                                                                                                    gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                           <  1:473692/251‑6 (MQ=255)                                                                                                                                                                                                    gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                           <  1:149110/251‑6 (MQ=255)                                                                                                                                                                                                       cccGCTGGTGCCACTGGTGGCGGGAATGCTGGAGCTGCCGGTGGCAAAAGTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccgatc                                        <  2:427747/250‑9 (MQ=255)                                                                                                                                                                                                                      ggtggCGGGAATGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTATTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCattat                          <  1:993768/250‑1 (MQ=255)                                                                                                                                                                                                                            gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                       <  2:118778/167‑1 (MQ=255)                                                                                                                                                                                                                            gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                       >  1:118778/1‑167 (MQ=255)                                                                                                                                                                                                                               aaTGCTGGAGCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                <  1:865995/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGCTGGATCTGCAGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCgg               <  1:783426/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTTGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGTCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttg        <  1:242311/247‑1 (MQ=255)                                                                                                                                                                                                                                           gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttgt       <  1:137733/248‑1 (MQ=255)                                                                                                                                                                                                                                            ccGGTGGCAAAATTCATTACGCCGAATATTTTCGGCTGCCTGCTGTGGCCACCGTTTTTCTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTtgtg   <  2:648349/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTgg  <  1:4234/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      GCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGG  >  NC_000913/71521‑72007

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