Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,159,530	G→A	5.4%	G539G (GGC→GGT)	ygiQ ←	Radical SAM superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,159,530	0	G	A	5.4%	240.4 / 4.0	74	G539G (GGC→GGT)	ygiQ	Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (41/28);  new base A (2/2);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

TGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCA  >  NC_000913/3159244‑3159817                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                tGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGcc                                                                                                                                                                                                                                                                                               <  2:52266/289‑1 (MQ=255) gtgtgCCCTTCGGTCTGAAGACGTCTTCGCTCTTTTAACCAATCTTCGCCAGCGGGTTTTTCCCGGGTGAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGcc                                                                                                                                                                                                                                                                                               <  2:469580/285‑1 (MQ=255)         ttCGGTACGAAGACGTCTTCGCTCTTATAACCGATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtca                                                                                                                                                                                                                                                                                       <  1:3/289‑1 (MQ=255)         ttCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtca                                                                                                                                                                                                                                                                                       >  2:6036/1‑289 (MQ=255)                 gAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCtt                                                                                                                                                                                                                                                                                >  1:317946/1‑288 (MQ=255)                       gTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTACAGCGGTCATAGCTGCCCATTCCCGGCTTCATTATTTTCGATAAc                                                                                                                                                                                                                                                                         >  2:106053/1‑289 (MQ=255)                         cTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACgg                                                                                                                                                                                                                                                                       <  1:223647/289‑1 (MQ=255)                           tCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGcc                                                                                                                                                                                                                                                                     <  2:289286/289‑1 (MQ=255)                                   ataACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATACCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCg                                                                                                                                                                                                                                                             >  1:125291/1‑289 (MQ=255)                                    taACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCgg                                                                                                                                                                                                                                                            >  1:37281/1‑289 (MQ=255)                                           tCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATaa                                                                                                                                                                                                                                                                          >  2:35329/1‑268 (MQ=255)                                           tCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATaa                                                                                                                                                                                                                                                                          <  1:35329/268‑1 (MQ=255)                                           tCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGtt                                                                                                                                                                                                                                                      >  2:479446/1‑288 (MQ=255)                                                gCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgggg                                                                                                                                                                                                                                                <  1:6036/289‑1 (MQ=255)                                                      gggTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAAt                                                                                                                                                                                                                                            <  1:479446/287‑1 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATCGCTGCCCATGCCCGGCTTCAACATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCCATTTTCGAATAAcaccccaccgtttgg                                                                                                                                                                                                                   >  2:17826/1‑275 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  2:29067/1‑289 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  2:494947/1‑289 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  2:125676/1‑289 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  1:486834/1‑289 (MQ=255)                                                                             ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  2:39799/1‑289 (MQ=255)                                                                               tggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATggg                                                                                                                                                                                                                  >  2:220740/1‑288 (MQ=255)                                                                               tggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATggg                                                                                                                                                                                                                  >  1:56/1‑288 (MQ=255)                                                                               tggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATggg                                                                                                                                                                                                                  >  1:523698/1‑288 (MQ=255)                                                                                      gTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAg                                                                                                                                                                                                           <  1:39799/288‑1 (MQ=255)                                                                                      gTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAg                                                                                                                                                                                                           <  2:37281/288‑1 (MQ=255)                                                                                                ggcgACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGata                                                                                                                                                                                                 <  1:220740/288‑1 (MQ=255)                                                                                                                 cTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACg                                                                                                                                                                               >  2:174763/1‑289 (MQ=255)                                                                                                                         ggTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATACCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTc                                                                                                                                                                       <  2:125291/289‑1 (MQ=255)                                                                                                                               ggCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCCGCTCTTTGATAGTGCGAGGATCTTCTACGGCTATGTcctcac                                                                                                                                                                  >  1:285597/1‑283 (MQ=255)                                                                                                                                gCGGAAGCGGTGCTTTTGCAGCCACAGCGCCAGATTCACCATATCTTCATCACGGGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTNATATAGCGCGGATCTTCTACGGCTATGTCATAACg                                                                                                                                                                 <  1:17826/288‑1 (MQ=255)                                                                                                                                          tGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCagag                                                                                                                                                       <  1:125676/288‑1 (MQ=255)                                                                                                                                          ggCTTGTTCAGCCACAGCGCCAGATTCACCATATCTTCAGCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCagag                                                                                                                                                       <  2:317946/287‑1 (MQ=255)                                                                                                                                           gCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAgg                                                                                                                                                      <  1:174763/288‑1 (MQ=255)                                                                                                                                                        caGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCttt                                                                                                                                        >  1:269576/1‑289 (MQ=255)                                                                                                                                                         aGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCtttt                                                                                                                                       >  1:375029/1‑289 (MQ=255)                                                                                                                                                         aGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCCATGTCATAACGTACACCAGAGGCCATCAGGA‑Cttttt                                                                                                                                      >  2:245963/1‑289 (MQ=255)                                                                                                                                                            gcCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTa                                                                                                                                     <  1:494947/288‑1 (MQ=255)                                                                                                                                                                 aTTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAccgcc                                                                                                                                                                                                                             >  1:20137/1‑195 (MQ=255)                                                                                                                                                                 aTTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAccgcc                                                                                                                                                                                                                             <  2:20137/195‑1 (MQ=255)                                                                                                                                                                     aCCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTACGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGTCTATGGCATAACGCACACCAGAGGCAATCAGGATCTTGTTAATGCCtttta                                                                                                                           >  1:89769/1‑287 (MQ=255)                                                                                                                                                                       cATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATACCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGa                                                                                                                         >  2:49650/1‑289 (MQ=255)                                                                                                                                                                           tCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCa                                                                                                                      >  2:501043/1‑288 (MQ=255)                                                                                                                                                                               catcaCGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACgcgc                                                                                                                 >  1:436311/1‑289 (MQ=255)                                                                                                                                                                                atccCGCTTACCGGGGTGGGGGGGAATGAAATAGGGGGTCAGATTCTGCTCTTTACTGGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCa                                                                                                                <  2:285597/289‑1 (MQ=255)                                                                                                                                                                                 tcaCGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCAc                                                                                                               >  1:96564/1‑289 (MQ=255)                                                                                                                                                                                 tcaCGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCAc                                                                                                               <  1:106053/289‑1 (MQ=255)                                                                                                                                                                                    cgcgTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGc                                                                                                            >  1:468885/1‑289 (MQ=255)                                                                                                                                                                                       gTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGa                                                                                                          <  2:523698/288‑1 (MQ=255)                                                                                                                                                                                       gTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATACCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGAt                                                                                                         <  1:49650/289‑1 (MQ=255)                                                                                                                                                                                         aCCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATag                                                                                                       >  2:340522/1‑289 (MQ=255)                                                                                                                                                                                                 gcgcGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGa                                                                                                >  2:327980/1‑288 (MQ=255)                                                                                                                                                                                                  cgcgGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGt                                                                                                                                                                                                                                                           <  1:311433/132‑1 (MQ=255)                                                                                                                                                                                                  cgcgGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGt                                                                                                                                                                                                                                                           >  2:311433/1‑132 (MQ=255)                                                                                                                                                                                                  cgcgGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATc                                                                                              >  1:379315/1‑289 (MQ=255)                                                                                                                                                                                                       agaTGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtgg                                                                                         >  1:12270/1‑289 (MQ=255)                                                                                                                                                                                                                   ggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCAcgcg                                                                                                                  >  2:98868/1‑252 (MQ=255)                                                                                                                                                                                                                   ggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCAcgcg                                                                                                                  <  1:98868/252‑1 (MQ=255)                                                                                                                                                                                                                             cTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGc                                                                   <  2:379315/289‑1 (MQ=255)                                                                                                                                                                                                                             cTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGc                                                                   <  2:486834/289‑1 (MQ=255)                                                                                                                                                                                                                             cTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCCCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGCCGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATCCCTTTCAGATCCCGCGCACGGCGATAGGGGTTTATCGTTGGTTCGTGGTTCGTTTCCAtgtt                                                                    >  2:85302/1‑287 (MQ=255)                                                                                                                                                                                                                                ctctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCgg                                                                 >  1:40856/1‑288 (MQ=255)                                                                                                                                                                                                                                ctctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCgg                                                                 >  2:163035/1‑288 (MQ=255)                                                                                                                                                                                                                                   tttACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACaa                                                             >  2:107935/1‑289 (MQ=255)                                                                                                                                                                                                                                     tACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACgg                                                                                                                                                                              >  1:250242/1‑174 (MQ=255)                                                                                                                                                                                                                                     tACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACgg                                                                                                                                                                              <  2:250242/174‑1 (MQ=255)                                                                                                                                                                                                                                          gcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCAcgcg                                                                                                                  <  1:49999/229‑1 (MQ=255)                                                                                                                                                                                                                                          gcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCAcgcg                                                                                                                  >  2:49999/1‑229 (MQ=255)                                                                                                                                                                                                                                                                    ctctTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAACTATCCGGATAAACGCACGACAAACGGCGACaa                            >  1:346964/1‑289 (MQ=255)                                                                                                                                                                                                                                                                    ctctTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACaa                            >  1:304483/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                 tCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   >  1:357256/1‑93 (MQ=255)                                                                                                                                                                                                                                                                                 tCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                                                                                                                                                                                                                   <  2:357256/93‑1 (MQ=255)                                                                                                                                                                                                                                                                                 tCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGTTTCGTGTCCATGTGCGGACAAATATCCGGCAAAACGCACGACAAACGGCGAAAAGTTTGTTCAgcgc               >  2:38342/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                       cTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCg         <  2:375029/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCa  >  2:246146/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                               gCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCa  >  2:79840/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                TGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCA  >  NC_000913/3159244‑3159817

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