Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	92.3%	241.5 / ‑4.6	78	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base A (45/27);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 3.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558660‑4559152                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       gctcccGTTGCCGTCGGAGCGGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                                                                                                                                                                                                     <  2:282170/247‑1 (MQ=255) gctACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATTCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCGGCAGGTCATAAATAGGCTGTAACGCCattt                                                                                                                                                                                                                                                     <  1:333399/251‑1 (MQ=255)         tGCCGGCGGAGCTGAGGGGGACGCGTGCCAGCGTAATACCCGCCTGAACGGCGCGGGCAATACCTGCTGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑c                                                                                                                                                                                                                                             <  2:78909/250‑1 (MQ=255)            cGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑GGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGccc                                                                                                                                                                                                                                         >  2:421907/1‑251 (MQ=255)                ggAGCTGAGGGTGACGCGTGCCTGCTGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGCCCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCAtgt                                                                                                                                                                                                                                     <  1:8978/251‑1 (MQ=255)                  aGCTGAGGGTGCCGCGTGCCAGCGGAATACCCGCCTGAACTGCGCGGGCAATACCTTCGGCAGGGGCGACCGGATCGTCAATGCTGCCGGTGATATCGATGGTGCCGCCTTTGTGAGCGAACTCCA‑GCGCCTGCTCAACCAACGG‑GACGTTGCGGAT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGTCCAGAAATAGGCTTTAACGCCTTTTTCCTGT‑TGCCCACGTgg                                                                                                                                                                                                                                   >  2:285823/1‑251 (MQ=255)                     tGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAac                                                                                                                                                                                                                                >  1:80677/1‑251 (MQ=255)                          gtgacggttGCCAGCGGAAGACCCGCCTGAACGGCGCGTGCAATAGCTTCGGCAGGGGCGACCGGTGCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCATCAGCTTGCTGATCGGCACATCGCAGTTTGCCAGCAGGGCATAAATAGGCTGTAACGCCATTTTACTTT‑CGCCCATGTGGAACACGGt                                                                                                                                                                                                                           <  2:140208/243‑1 (MQ=255)                             aCGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                        <  2:366092/251‑1 (MQ=255)                                cgTGCAAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCACATGTGGAACACGGTGACGcc                                                                                                                                                                                                                     >  2:436255/1‑251 (MQ=255)                                        gggAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTAc                                                                                                                                                                                                             <  1:89664/250‑1 (MQ=255)                                         ggAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTAcc                                                                                                                                                                                                            >  1:158440/1‑251 (MQ=255)                                                      gAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGAGACGCCAGGTTTACCGCCGAGCAAACCg                                                                                                                                                                                               >  1:35015/1‑251 (MQ=255)                                                             gcgGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc                                                                                                                                                                                        >  1:367386/1‑251 (MQ=255)                                                               gggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGAGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg                                                                                                                                                                                      >  2:32915/1‑251 (MQ=255)                                                                            ggCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCa                                                                                                                                                                         <  1:40733/251‑1 (MQ=255)                                                                               gggggCGCCCGGTGCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCGTTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCatag                                                                                                                                                                      <  2:273634/250‑2 (MQ=255)                                                                                  ggCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAACCAACGG‑TACGTTGCGGTT‑ACC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATCTTgg                                                                                                                                                                   >  1:377923/1‑251 (MQ=255)                                                                                       ccGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAgg                                                                                                                                                              >  1:208919/1‑251 (MQ=255)                                                                                              gTCAATGCTGCTGGAGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGAGATaa                                                                                                                                                        >  2:346829/1‑250 (MQ=255)                                                                                                        gggtggATATCGATGCTGCCGCCGTGCCGCGGGAACTCCA‑GCTCGTGCCAAAACAACGG‑TCCGTTTCGGTT‑ACC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACTTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATGTTAC‑GTGCGCCCAGGTGGACCTCGGTGCCGCCGGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTTATAAACGTCCGGTGc                                                                                                                                             <  2:39887/246‑1 (MQ=255)                                                                                                        ctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGc                                                                                                                                             >  1:72266/1‑251 (MQ=255)                                                                                                             gATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                          >  2:48498/1‑249 (MQ=255)                                                                                                                atCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                     >  2:314090/1‑251 (MQ=255)                                                                                                                atCGATGGGGCGGCCTTGGCGCGCGAACTCCG‑GCGCCTGCCAAAACAGCGG‑GACGT‑GGGGTT‑GAC‑GGTGGGTGGGCAGAAGTTTGCGGAGCGGCACATCGAAGTTTTCCAGCAGGTAATAAATAGGCTGTAACGCCTTTTTACGGT‑CGCCCAGGTGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGGGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                     <  2:140112/251‑1 (MQ=255)                                                                                                                      ggTGCCGCCTTTGCGCGCGACCTCCA‑GCGCCTGCTCAACCAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                               <  1:432593/251‑1 (MQ=255)                                                                                                                      ggTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                               >  2:121164/1‑251 (MQ=255)                                                                                                                      ggTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCGGCTCAAACAAGGG‑TACTTTGCGGTT‑ATC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATGGGCTGTAACGCATTTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGa                                                                                                                                <  2:131968/250‑1 (MQ=255)                                                                                                                      ggTGCCGCCTGTCCGCGCGAAC‑CCACGCGCCGGCTCAAA‑AACGGGTGC‑TTGGCGTT‑AACCG‑TGGGTCGGCAGCAGCTGGCTGATCGGCACATCGCAGTTTGCCAGCAGGAAATAAAAAGGCTTTAACGCCATTTACCTGT‑CGCCCTTGTGGAACACGGTGACGCCGGTTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCGGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                               <  2:156056/251‑1 (MQ=255)                                                                                                                       gTGCCCCCTTTGCGCGCTAACTCCA‑GCGCCTGCTAAACCACCGG‑TCCGTTGCGGTT‑AAC‑G‑TGGGTCGGCATCAGCTTGCTGATCGGCACATCCCAGTTTTCCAGCGGGTCATAAAGAGGCTGTAACCCCATTTTACTGT‑CGCCCCTGGGGACCACGGGGCCGCCAGGTTTACCGCCGAGAAAACCGCCAACGCGGGATTCGGCCGCCATATTGGCCAGATGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                              <  2:289604/251‑1 (MQ=255)                                                                                                                              cTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATg                                                                                                                       >  1:84403/1‑251 (MQ=255)                                                                                                                                  ggcgcgggaaccccatggcCTGCTCAAACAACGG‑GACGTTGCGGTT‑AAC‑G‑GGGGTCGGCAGAAGCTTGCTGATCGGCACATCGAAGTTGTCAGGCAGGTCATAAATAGGCTGTAACGCCTTTTTACGGT‑CGCCCATGTGGACCACGGTGACCCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCAGAACGGTGATCAGAGATGgc                                                                                                                     <  2:159441/233‑1 (MQ=255)                                                                                                                                    gcgcgAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTG‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAt                                                                                                                 <  1:26602/251‑1 (MQ=255)                                                                                                                                    gcgcgAACTCCA‑GCGACTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑GGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAt                                                                                                                 <  1:41455/251‑1 (MQ=255)                                                                                                                                         aaccccA‑GCGCCTGCTCAAACAACGG‑TACGTGGCGGTT‑AAG‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCCCAGTTTTCCAGCAGGTCATAAATAGGCGGTAACGCCTTTTTACTGT‑CCCCCAGGTGGAACACGGTGACGCCAGGTTACCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCAc                                                                                                            <  2:399955/247‑1 (MQ=255)                                                                                                                                               a‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAAt                                                                                                      >  1:185984/1‑251 (MQ=255)                                                                                                                                                   gcCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATcac                                                                                                   >  1:61438/1‑251 (MQ=255)                                                                                                                                                          cAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatc                                                                                             <  1:61469/250‑1 (MQ=255)                                                                                                                                                          aaaaCAACGG‑TACGTTGCGGTT‑ACC‑G‑GGGGTCGGCAGCAGCTGGCTGATCGGCACAGCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcg                                                                                            <  2:138410/250‑1 (MQ=255)                                                                                                                                                           aaaCAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                                                                                                                                                                                                    <  1:315799/146‑1 (MQ=255)                                                                                                                                                           aaaCAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                                                                                                                                                                                                    >  2:315799/1‑146 (MQ=255)                                                                                                                                                                acGG‑TACGTTGCGGTT‑AAC‑G‑TGGGGCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGGGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATc                                                                                      <  2:282149/251‑1 (MQ=255)                                                                                                                                                                      aCGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCAc                                                                                 <  1:327480/251‑1 (MQ=255)                                                                                                                                                                       aggtGCGGTTTAAC‑G‑TGGGTGGGCAGAAGCGTGTTGATCGACACATCGCAGTTTTCCAGCGGGTAATAAAGAGGCTGTAACGCCATTTTACTGT‑CGCCCAGGTGGAACCCGGTGACGCCAGGTTTACCGCCGAGAAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCa                                                                                  <  2:110977/247‑1 (MQ=255)                                                                                                                                                                          tGCGGTT‑AAG‑T‑TGGGTCGGCAGCCGCTTGCTGTTCGGCCCATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCTGGTGGAACACGGGGACGCCGGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTc                                                                              <  2:253140/250‑1 (MQ=255)                                                                                                                                                                           gCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCAGAACGGTGCTCAGAGATGGCGCATTGCACGCCAATCACACGATCGATAATCGCCACGTCtt                                                                            >  2:44365/1‑251 (MQ=255)                                                                                                                                                                             ggTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                          >  1:286380/1‑251 (MQ=255)                                                                                                                                                                                t‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTc                                                                        <  1:364562/250‑1 (MQ=255)                                                                                                                                                                                  aaC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGcc                                                                                                         >  1:170903/1‑216 (MQ=255)                                                                                                                                                                                  aaC‑G‑TGGGCCGGCAGAAGCTTGCGGATCGGCACATCGCAGTTTTCCAGCAGGTCAAAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCGGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGcc                                                                                                         <  2:170903/216‑1 (MQ=255)                                                                                                                                                                                        tGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGCGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCaagg                                                                   >  2:487/1‑247 (MQ=255)                                                                                                                                                                                        tGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCAAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGCGATAAACGTCCGGTGCGGCAGAACGGTGATAAGAGCTGGCGCATTTCACGCACATCACACGATCGATAATAGCCACGCATTTTTCCACgg                                                                   >  2:127543/1‑250 (MQ=255)                                                                                                                                                                                            tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGCGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                              >  1:336011/1‑251 (MQ=255)                                                                                                                                                                                              ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTCTTCCACGGAACCCGt                                                            >  2:262262/1‑251 (MQ=255)                                                                                                                                                                                               gcagcgGCTTGCTGATCGCCACATCGCAGTTTTCCAGCAGGTCCTAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGACAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGTTCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                            >  2:329455/1‑250 (MQ=255)                                                                                                                                                                                                cagcagCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTa                                                           >  1:203650/1‑250 (MQ=255)                                                                                                                                                                                                 agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGCTTACCGCCGAGCAAACCGCCAACGCCGGATTCTGCCGCCATATTGGCCAGGTGATACACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCCCGCCCATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTaa                                                          >  1:123003/1‑250 (MQ=255)                                                                                                                                                                                                  gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATg                                                        <  1:329488/251‑1 (MQ=255)                                                                                                                                                                                                      cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGc                                                    >  1:415039/1‑251 (MQ=255)                                                                                                                                                                                                             aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAgg                                             >  2:43738/1‑251 (MQ=255)                                                                                                                                                                                                                ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGAc                                          >  1:12113/1‑251 (MQ=255)                                                                                                                                                                                                                  cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTTCTGT‑CGCCCATGTGGAACACGTTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                        <  1:118654/251‑1 (MQ=255)                                                                                                                                                                                                                     aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGAt                                     >  1:49/1‑251 (MQ=255)                                                                                                                                                                                                                          aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcg                                >  1:375885/1‑251 (MQ=255)                                                                                                                                                                                                                          aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcg                                >  2:413248/1‑251 (MQ=255)                                                                                                                                                                                                                                caagcagGTAATAAATGGGCTGTACCTCCTTTTTGCTGT‑CGCCCATGTGTAACACGGTGACGCCTGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGGCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCTCCACGTCTTTGTCCACGGAACCCGTAATGGTGCTGGAAGGGACATGATAAGCGCCGGTc                          <  1:44911/249‑1 (MQ=255)                                                                                                                                                                                                                                 cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGt                           <  1:57963/249‑1 (MQ=255)                                                                                                                                                                                                                                 cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGt                           >  2:57963/1‑249 (MQ=255)                                                                                                                                                                                                                                 cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGACAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                          >  1:135215/1‑250 (MQ=255)                                                                                                                                                                                                                                  agcagGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                          >  1:56577/1‑249 (MQ=255)                                                                                                                                                                                                                                     agGTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCAt                     >  1:26892/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTCATAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATcc                   >  1:172031/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                >  2:64384/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTAAATAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                >  1:132310/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaaaTAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc               <  1:237998/250‑1 (MQ=255)                                                                                                                                                                                                                                            aaaTAGGCTGTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgcg              >  2:96023/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTAGGCTGTAACGCCATTTTACTGT‑CGCCCAGGGGGAAACCGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt            <  1:149570/251‑1 (MQ=255)                                                                                                                                                                                                                                                     gTAACGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCct     >  1:57067/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aCGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttct  >  1:351230/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aCGCCATTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttct  >  2:339290/1‑251 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       GCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCA‑GCGCCTGCTCAAACAACGG‑TACGTTGCGGTT‑AAC‑G‑TGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGT‑CGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558660‑4559152

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