Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,946	T→C	Q464Q (CAA→CAG)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,946	0	T	C	95.9%	267.8 / 0.8	99	Q464Q (CAA→CAG)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (48/45);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGAT  >  NC_000913/66700‑67191                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGgc                                                                                                                                                                                                                                                      <  1:485226/251‑1 (MQ=25)         gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCg                                                                                                                                                                                                                                              <  1:203741/251‑1 (MQ=25)                    gaaTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTAGCAAT‑CACCGGGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAg                                                                                                                                                                                                                                   <  2:295393/249‑1 (MQ=25)                     cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgc                                                                                                                                                                                                                                  <  1:512529/251‑2 (MQ=37)                         aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCCTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTTGGGTTCAgcgcct                                                                                                                                                                                                                              >  1:173148/1‑246 (MQ=37)                          ggCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGGGATTTCAATGTCGGG‑CATCTCGGCGAACAGGCGCATATCGTTGAGGGTCAGCGCATg                                                                                                                                                                                                                             >  1:637094/1‑251 (MQ=25)                          ggCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCTTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATg                                                                                                                                                                                                                             <  1:699048/251‑1 (MQ=255)                              ttGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCG‑GAGATCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑AATCTTAGCGAACCGGCGCataacggttaggatcagcgccaggctg                                                                                                                                                                                                                         >  1:380731/1‑227 (MQ=21)                                      cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCCGCGCATCGC‑T‑GAAGACggg                                                                                                                                                                                                               >  1:428156/1‑250 (MQ=37)                                       ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACCGGCGCATATCGTTGGGGGTCCGCGCATGGC‑T‑GAAGACggg                                                                                                                                                                                                               >  1:71588/1‑249 (MQ=25)                                             gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTgc                                                                                                                                                                                                        <  1:31632/251‑1 (MQ=38)                                                   cTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCCATGTCGTG‑CATTTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGccc                                                                                                                                                                                                  >  1:234604/1‑250 (MQ=37)                                                      aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGGTTTTAAT‑CAACGTGTTTTCAATTTC‑TGTCATCTGGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCCTGGC‑TTGCAG‑CGGTCTGGTGCGCGccaccc                                                                                                                                                                                               >  2:298902/1‑250 (MQ=19)                                                         cGGCCTACAAAATCGCAACGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATTAATACACGGT‑ATTTCAATGTCGTG‑CATCTTGGCGTACTCGCGCATATCGGTGGGGGTCAGCGCATGGC‑T‑GGAGACGGGATTGGGCGGGCCACCcgcg                                                                                                                                                                                            >  2:347379/1‑251 (MQ=255)                                                           gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACAAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAG‑TCGGTATCAAA‑TAACGAGATTTCATTATCGTT‑CATTTCGGCGAACTGGCCCCTATAGGTGAGGGTTACGGCATGGC‑T‑GAAGACgggagaggtggtgccacccgca                                                                                                                                                                                            >  1:370894/1‑229 (MQ=255)                                                               aCAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCAGCGCGTCTTTAAACGCTGGAAGGCGGGTA‑TCGTTATCCCA‑CACCGTGTTTTCCATATAGTG‑CATATCGGCGAATTGGCGCATATCGTTTGGGGTCAGCGCCTGTCAT‑GA‑GACGGTATGGGGGGGGCCACCCGCGAGGATc                                                                                                                                                                                      >  2:534223/1‑251 (MQ=255)                                                                  aaagcGCAGCGTGTAGGCCAGATAAGACGCGGCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCAc                                                                                                                                                                                   <  2:301142/247‑1 (MQ=255)                                                                  aaaacGCAGCGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCAc                                                                                                                                                                                   <  2:328889/247‑1 (MQ=255)                                                                  aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCAc                                                                                                                                                                                   >  1:484964/1‑251 (MQ=255)                                                                        cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg                                                                                                                                                                             <  2:355268/251‑1 (MQ=255)                                                                         aGCGTGTAGGCCTGACAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCTGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCgg                                                                                                                                                                            >  1:432032/1‑251 (MQ=255)                                                                            gtgtAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGGTTCGGa                                                                                                                                                                           >  1:478823/1‑249 (MQ=255)                                                                            gtgtAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGAATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGa                                                                                                                                                                           >  1:66449/1‑249 (MQ=255)                                                                                aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGCTTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGTGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGGGCGCCCCCAGCGGGGAGCCCCGCTTCGGAAGCAgg                                                                                                                                                                     >  2:381247/1‑250 (MQ=255)                                                                                 ggccttaTAAGACGCGGCAGCGTCGCATCAGTCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGGCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGtt                                                                                                                                                                    <  2:238413/245‑1 (MQ=255)                                                                                 ggCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                                                                                                                                                                                                 <  1:586858/192‑3 (MQ=255)                                                                                     tgATAAAACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAA‑TCGTTATCAAT‑CACCGTGGTTTCAATGCAGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCAGGA‑T‑TAAGACGGTAAGTTTCGCGCCCCCATCGAGGGTCCACGCTTCGGGAGCAGTTGGc                                                                                                                                                                 >  1:127941/2‑250 (MQ=255)                                                                                        tAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACTGGCGCATAACGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGGAGCAGTTGGCAGa                                                                                                                                                              >  1:296606/1‑250 (MQ=255)                                                                                         aaGACGCGTCAGCGTTGTATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                                                                                                                                                             <  2:663465/250‑1 (MQ=255)                                                                                         aaGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                                                                                                                                                             <  1:381247/250‑1 (MQ=255)                                                                                          aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc                                                                                                                                                           <  1:671835/251‑1 (MQ=255)                                                                                          aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc                                                                                                                                                           <  2:194814/251‑1 (MQ=255)                                                                                            aCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCgg                                                                                                                                                         <  1:264377/251‑1 (MQ=255)                                                                                                  cAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTgcg                                                                                                                                                   <  1:290081/250‑1 (MQ=255)                                                                                                     cGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAACGCATGGC‑T‑GAAGACGGGATGGGGCGCGCCACCAGGGCAGATCCACGCCTTGGAAGcgggttggagatacgggttcgcttt                                                                                                                                                >  2:454739/1‑227 (MQ=255)                                                                                                      gTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTc                                                                                                                                              <  1:181078/251‑1 (MQ=255)                                                                                                                        cacctgATGCGACTACTTAGCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTT‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCatc                                                                                                                            <  2:285378/246‑3 (MQ=255)                                                                                                                         acctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGCG‑CATCTCGGCGGACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGGCCGTATTGGGCCCGCCACCAGAGAGGATCCACGCCTCCGGAGCCGTTGGGAGGTCCGGGTTG‑GCTTTccccaccgccatcgcccccc                                                                                                                           >  2:625289/5‑233 (MQ=255)                                                                                                                          cctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCgg                                                                                                                          >  1:318767/4‑251 (MQ=255)                                                                                                                          cctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCgg                                                                                                                          <  1:347379/248‑1 (MQ=255)                                                                                                                             gATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCa                                                                                                                        <  2:637094/250‑1 (MQ=255)                                                                                                                                 cgactactgatctacgaAATCCGTAATACACATCGTTCCAGCGCTGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCGCTGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                   <  2:622112/238‑1 (MQ=255)                                                                                                                                 cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                   <  1:514594/248‑1 (MQ=255)                                                                                                                                  gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGACACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCCCATTCGCCACCGGCAGTTTCg                                                                                                                  >  1:129301/3‑251 (MQ=255)                                                                                                                                  gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCATTTGGCAGATCCGG‑TTgc                                                                                                                                                    <  1:559675/215‑1 (MQ=255)                                                                                                                                   acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAA‑CACCGTGATTTCAATGTCGTG‑CATCCCGGCGAACTTGCGCCTATCGTTGAGGTTCAGCGCATGGC‑T‑GCAGACGGTATGGGGCGCGCCACCCGCGAGGCACGCCGCTCCGGAAGCAATTGGCCGGTCCGG‑TTGCGATTTCCACCACGCCTTCGCCACCGGCAggctctg                                                                                                                 >  2:22628/2‑245 (MQ=255)                                                                                                                                             gacgaAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGc                                                                                                       <  1:702729/251‑1 (MQ=255)                                                                                                                                                 aaatccGTAATACACTTCGATCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGTTA‑TCGTTATCAAT‑CACCGTGGTTTCAATGTGGTG‑CATCTCGGCGGACTGGCCCATATTGCTGAGGTTTCGCGCACTGC‑T‑GGAGACGGTATGGTGCGCGCCACCCGCGAGGCTCCCCGCCTCGGGAGCAGTTGGCTGGTCCGG‑TGGCGCGTTCCCCAACGCATTCCCCACCCGCaagttagcccgggggggcgggg                                                                                                   >  2:467533/5‑230 (MQ=255)                                                                                                                                                      cGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCCAA‑CACCGTGCTTTCAATGTCGTG‑CCTCTCGGCGAACTGGAGCATATCGTTGAGGTTCAGCGCATGGC‑T‑TAAGACGGTATGGTGCGAGCCACCCGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACCGACCATTCGCCCCCGGCAGCGTCTGCAGGGGGTGCGGGGtttat                                                                                              >  2:521982/1‑249 (MQ=255)                                                                                                                                                      cGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTg                                                                                                                                                     <  1:222797/196‑1 (MQ=255)                                                                                                                                                      cGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGGTGAGGTTCAGCGCATGGC‑C‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGGTGGCAGATCCGG‑TTg                                                                                                                                                     >  2:222797/1‑196 (MQ=255)                                                                                                                                                      cGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGCTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGCACTGGCGCATATCGTTGAGGTTCAGCGCCTGGC‑T‑GAAGACGGTATGGGGAGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCCCATTCGCCACCGGCAGTTTCGGCCGGGAGTGCgggggtttt                                                                                              >  1:46218/1‑244 (MQ=255)                                                                                                                                                      cGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCCATGTCGTG‑CATCTCGGCGAACTGGCGCCAATCGTTGAGGGTCAGCGCCTGGC‑T‑GAAGACGGGATGGGTCGCGCCCCCCGCGGGGATCCCCGCTTCGGAAGCAGTTGGCAGAACCGG‑TTGGGCTTTCCACTGGGCCTTGCCCACCGGGAGTTTCGGGAgggggagtggggttttt                                                                                              >  2:282663/1‑237 (MQ=255)                                                                                                                                                       gTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAA‑TTGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGGGGAACCACGCTTCGGAAGCAGATGGGAGATACGG‑TTGCGCTTTCCACAACGCCTTCGCCCCCGGCAGGTTCGGCAgggggggcgggggttttc                                                                                             >  2:497930/1‑236 (MQ=255)                                                                                                                                                         aaTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCAc                                                                                            <  1:534223/250‑1 (MQ=255)                                                                                                                                                           tACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CAACGTGATTTCCATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑GTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTCGGCAGGGAggtagggggtttcccacg                                                                                         >  1:395886/1‑234 (MQ=255)                                                                                                                                                                  cGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGGATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGGAGGGAGTgggggggtttccacgtgtcggt                                                                                   >  1:550948/1‑229 (MQ=255)                                                                                                                                                                   gTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAA‑AACCGTGATTTCAATGTCGTG‑CATCTCGGCGAAATGGAGCATATCGTTGAGGTTCAGCGCATGTC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCCCGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCACCACCGGCAGTTTATGCCGGGAGAGCGGTGgtttacccgtggcactac                                                                                 >  1:407548/1‑237 (MQ=255)                                                                                                                                                                    ttCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGAGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCGGTTTCGGCAGGGGGGGCGGGGGTTTCACCGAGTCGCTGc                                                                                 >  2:204136/1‑250 (MQ=255)                                                                                                                                                                        agcgcagcgcGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGGTTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCCGCGCATGGC‑T‑GAAGACGGGATGGGGCGCGCCACCAGCGCGGATCCCCGCTTCGGAAGCAGTTGGCAGATCCGG‑GTTCGCTTTCCACAGCGCATTTCCCACCGGCAGGTTCCGCCGGGGGGGCGGGGTTTTCACCTTGTCgcgccagat                                                                             >  2:663195/1‑242 (MQ=255)                                                                                                                                                                             agcgcGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGGGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATAGTTGTGGGTCAGAGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCCGTTGGCAGAACCGG‑TTGCGCTTTCCACAGAGCATTAGCACCCGGCAGTTTCGGCAGGGGGGGCGGTGGTTTTACCGTGTCGAAACAGTTAAcc                                                                         >  1:657410/1‑249 (MQ=255)                                                                                                                                                                               cgcgTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCcac                                                                                                                                           <  1:189007/181‑1 (MQ=255)                                                                                                                                                                               cgcgTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAgca                                                                     <  1:292579/251‑3 (MQ=255)                                                                                                                                                                                gcgTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGGTGAGGGTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGGAGGGGGGGCGGGGGTTTCACCggggcgataaaggtaaacaga                                                                      >  1:336907/1‑229 (MQ=255)                                                                                                                                                                                  gTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGAc                                                                  <  1:565098/251‑1 (MQ=255)                                                                                                                                                                                  gTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGGTCCACGCTTCGGGAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCCGCAGGGAGTGCGGTGGTTTTACCGTGTCCGTACAGTTAACCAGCAGa                                                                   >  1:104578/1‑250 (MQ=255)                                                                                                                                                                                  gTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CAACGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGGTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTCGGCAGGGGGAGGGGTGTTTTTCCCGTGTCGATACAGTTAACCAGCAGAc                                                                  >  1:227312/1‑251 (MQ=255)                                                                                                                                                                                    cTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑GTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGGTTTCACCGTGTCGATACAGTTAACCAGCAGACgg                                                                >  1:166480/1‑251 (MQ=255)                                                                                                                                                                                          aCGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTTAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTTGGCAGGGAGTGCGGGGTTTTTACCGTTTCGGTACAGGTAACCAACAGACGGTaaagg                                                          >  1:653375/1‑248 (MQ=255)                                                                                                                                                                                            gCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                        <  1:202882/251‑1 (MQ=255)                                                                                                                                                                                                        ggtatCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGCTTTCACCGTGTCGATACAGTTAACCTGCAGACGGTAACGATCGCCGAGatca                                              >  1:70243/5‑250 (MQ=255)                                                                                                                                                                                                         gtatCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGAGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGTCAGTTTCGGCAGGGAGTGCGGGGTTTTCACCGTGTCGATACAGTTCACCCGCAGACGGGAAAGATCTCCGGGatcaa                                             >  2:21207/4‑250 (MQ=255)                                                                                                                                                                                                           atCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGgtcaatc                                           >  1:385415/2‑250 (MQ=255)                                                                                                                                                                                                             cgtaagcaat‑caCCGTGAGTTCAATGTCGTG‑CATCGCGGCGAACTGGCGCATATCGTGGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcag                                         <  2:318767/244‑2 (MQ=255)                                                                                                                                                                                                             cGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcagg                                        <  1:369377/251‑3 (MQ=255)                                                                                                                                                                                                                tATCAAT‑CACCGTGATTTCAATGTCGGG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTTCGCTTTCCACAGCGCATTCGCCACCGGCAGTTCCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCt                                      <  1:341586/250‑1 (MQ=255)                                                                                                                                                                                                                 atcaat‑caCCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTgg                                    >  1:205240/1‑251 (MQ=255)                                                                                                                                                                                                                   caat‑caCCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCg                                  <  2:280445/251‑1 (MQ=255)                                                                                                                                                                                                                    aat‑caCCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCCGGGAGGGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTCAAGATCTCCGAGATCAATCCGGCTGGCg                                  >  2:486977/1‑250 (MQ=255)                                                                                                                                                                                                                     at‑caCCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGAc                                >  1:527212/1‑251 (MQ=255)                                                                                                                                                                                                                          ccGTGATTTCAATGTCGGG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATc                            <  2:162451/251‑1 (MQ=255)                                                                                                                                                                                                                            gTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCCTTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAAAGATCGCCGCGTTCCATCAGGATGGCGCCAATCGc                          >  2:119405/1‑251 (MQ=255)                                                                                                                                                                                                                               gattgcaagggcgttcATCTCGGCGAACTGGCGCATATCGTTGGTGTTCATCTCATGGC‑T‑GAAGACGGTATGTTGCGCTCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGGGCTGTGTTTTCTCCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGTGACAATCGCt                         <  2:671791/235‑1 (MQ=255)                                                                                                                                                                                                                                 ttCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                   >  1:297667/1‑157 (MQ=255)                                                                                                                                                                                                                                 ttCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                   <  2:297667/157‑1 (MQ=255)                                                                                                                                                                                                                                  tCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATTGCCGAGATCAATCAGGCTGGCGACAATCGCTggac                     <  1:511851/250‑3 (MQ=255)                                                                                                                                                                                                                                    aaTGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGGTTTCACCGTTTCGATACAGTTTACCAGCAGACGGGCACGATCGCCGGGATCCATCCGGCTGGCGACCATCGCTggacagg                  >  2:166779/1‑246 (MQ=255)                                                                                                                                                                                                                                       gTCGTG‑CATCTCGGCGACCTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCCCCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGCGTGCGGGGTTTTCACCGAGTCGGTCAAGTTAACAAGCAGACGGGAACGGTCGCCGGGATCAAACAGGCTGGCGACAATAGCTGGGCCGGttt               >  2:677786/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGttt               <  1:490953/251‑1 (MQ=255)                                                                                                                                                                                                                                       gTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGtt                <  1:169781/250‑1 (MQ=255)                                                                                                                                                                                                                                        gcGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTg              <  1:625289/250‑1 (MQ=255)                                                                                                                                                                                                                                         cGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTgg             >  1:269347/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ctcGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTAGGCAGGGCGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGAACGCCGGGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTgg      >  2:308995/1‑250 (MQ=255)                                                                                                                                                                                                                                                   cGGCGAACTGGCGCATTTCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTagt        <  2:478823/247‑3 (MQ=255)                                                                                                                                                                                                                                                   cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTAtt        <  1:55003/247‑1 (MQ=255)                                                                                                                                                                                                                                                   cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGGTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa      <  1:163051/249‑1 (MQ=255)                                                                                                                                                                                                                                                    ggCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa      >  1:263108/1‑248 (MQ=255)                                                                                                                                                                                                                                                     gCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGAt  <  1:512616/251‑1 (MQ=255)                                                                                                                                                                                                                                                      cGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa      >  1:568041/1‑246 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         GATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGAT  >  NC_000913/66700‑67191

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