Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	635,136	C→T	11.3%	E81K (GAA→AAA)	ybdM ←	Spo0J family protein, ParB‑like nuclease domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	635,136	0	C	T	11.3%	209.3 / 2.6	72	E81K (GAA→AAA)	ybdM	Spo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (26/37);  new base T (6/3);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCC  >  NC_000913/634915‑635374                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   ggTGATTTGAAGGCGTCCACGAGCGCGATTATGGCGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATATCCTTTTAAGCGTAATTTTAATGATGATGATCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTgcagc                                                                                                                                                                                                                                   <  1:600512/238‑1 (MQ=255)    gATTTGGTGGCCTCCCCGCACGCGATTATGACGAATAGTCCGTGTGATGCTTTGATTACGCGTACCTTCCAGACAGGTCCCCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGt                                                                                                                                                                                                                                <  1:474909/238‑1 (MQ=255)         tgtgggttCCCCGAGCGCGTTTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGttttttt                                                                                                                                                                                                                          <  2:317734/232‑1 (MQ=255)               gTCCACGAGCGCGATTATGACGAAAAGGCGGTGCGATTCGTTTAATACGCGTACCTTCCAGACAGGTCCCCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCggg                                                                                                                                                                                                                     <  1:442921/238‑1 (MQ=255)               ctCCACGGGCGCTATTATGACGAATAGTCGGGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCggg                                                                                                                                                                                                                     <  1:347364/237‑1 (MQ=255)                       gcgcGATTAAGGCGAATAGTCGCGGCGAGGCGTTGATTACGCGTACCTTCCAGGAAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATTGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                             <  1:412843/238‑1 (MQ=255)                        cgggATTAAGAACAATAGTCGCGGGTATGGGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGGGTAATTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                             <  1:1061632/234‑1 (MQ=255)                         gcgATTATGACGAATAGGCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACg                                                                                                                                                                                                           <  1:368577/238‑1 (MQ=255)                          cgATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGGAAACACAT‑CACGAGTTCCATGTGACTTTTTTTCTCAT‑TGTGTACC‑C‑CAAAGTGTTGGGGGAAACCCCCTCTTTTTA‑TGTGGTTTTTCAGCAGTTTTTTTTTGGGGGGGGGCCCgtt                                                                                                                                                                                                         >  2:350974/1‑239 (MQ=255)                              tttggcGAATTGTCGGGGGGATTCGTTGATTTCGGGTTCCTTCCAGACCGGTCCCCGGGAAATTACCTTTTAAGCGTAATTTTAATGGTGATTGGCCTTTACCAATTTCATGGCGGTGAAAAC‑CAACCACGATTTCCAAGGCATTTTTATCCGTATGTGTT‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTTAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                     <  1:350974/234‑1 (MQ=255)                              tttGACGAATAGTCGCGGCGATGCTTTGATTCCGCGTTCCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATTATGGGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                     <  1:900960/237‑1 (MQ=255)                              tATGACGAATAGTCGCTGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                     <  1:853980/239‑1 (MQ=255)                                          tCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAA‑CATCCCCCATTTTCACGGCCTTTTTTTCCCTATGTGTG‑GAC‑A‑CAATTGAGTGGGTTAAAACCTCCATTTCAA‑TTGATTTTTTCAGCAGGTTTTTTTTTGGGGGGGCCCCCTTGTTTGGGttttttgaa                                                                                                                                                                                         >  2:260880/1‑232 (MQ=255)                                           gggggggttgggTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑At                                                                                                                                                                                       <  2:725721/228‑1 (MQ=255)                                                         gtttACGCGGATCTGCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGGGCATTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTgg                                                                                                                                                                         <  1:482581/237‑1 (MQ=255)                                                          aTTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CAACCACAAATTTCATGGCAGTTTTTTTCGTATGTGTG‑ACC‑ACCATTGGGTTGTG‑GAAACCATCTTTTTTAT‑TGGATTTTTTTAGCAGGTTTTTTTTGGGTGGGGGCGCCTTGTTTTGGTTTTGATC‑ATTAGGCCTAAtttt                                                                                                                                                                          >  1:461901/1‑236 (MQ=255)                                                                    ccTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACCATTTTCCTGGCCTTTTTTTCCCTATGTTTG‑ACC‑A‑CAATGGGTTGGGGAAAACCATCTATTTTAA‑TTGTTTTTTTTACAAGTTTTTTTTTTGGGGTGGGCCCGTTTTTTTGGTTGTAAAT‑ATTAGGAAAAAATTTGTT‑G‑tttttta                                                                                                                                                             >  2:372463/1‑236 (MQ=255)                                                                    ccTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAA‑CATCCACGATTTTCCTGGCATTTTTATCCGTATGTGTG‑GCC‑A‑CAATGGGGTGGGTAAAACCATCTATTTCAA‑TGGGTTTTTTCAGCAGTTTTTTTTCGGGGGGTGCCACGTTGTTTTGGTTTTTATC‑ATTAGGCTTAAGTTTGCT‑G‑TTTTTaa                                                                                                                                                             >  1:409979/1‑238 (MQ=255)                                                                    ccTTCCAAACAGTGCACCGGAAAATAGCCTTTTAAGCGTAATTTGAATGATGATGAGCCTTTAGCAATTGCATGGCGGTGAAAAC‑CATCCACGATTTCCAGGGCATTTTTATCCGTAGGTGTG‑ACC‑A‑CAATGGTTTGGGTAAAACCATCTATTTCAA‑TGTATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAAc                                                                                                                                                            <  1:246050/239‑1 (MQ=255)                                                                         cAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCcaca                                                                                                                                                       <  1:430788/239‑1 (MQ=255)                                                                               ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTTCATGGCATTTTTTTCCGTATGTTTG‑GCC‑A‑CAATGGGTTGGGGAAAACCAACTATTTCAA‑TGGGTTTTTTCAGCAGTTTTTTTTCGGGGGGTGGCCCGTTTTTTTGGTTTTAATC‑ATTAGGGATAAATTGGCT‑G‑TTTTTTAACcaccaaaagg                                                                                                                                                 >  1:734141/1‑232 (MQ=255)                                                                               ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAA‑CATCCACCGTTTTCATGGCAATTTTTTTCCTATTTTTG‑GCCAA‑CAATGGG‑TGGGGAAAACCATCTATTTTTA‑TTGGTTTTTTTAGCAGATTTTTTTTTGGGGGGGCCCAGTGGTTTGGGTTGTGTAC‑ATTTTGCTTAAATTGGGTGG‑TTTTTAACCCCCAAAAca                                                                                                                                                  >  1:16907/1‑238 (MQ=255)                                                                               ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAA‑CATCACCGCTTTTCCTTGCATTTTTTATTGTATGTGTG‑AAC‑A‑CAACAGGTTTGGTAAAAACATCCATTTTAC‑TTTAATTTTTTACCAGTTTTTTTTTTGGGGTGGTCACGTTTTTTTGGTTTTGATA‑TATAGGGGTAGATTGGTT‑GTTTTTTAACCcaacaaaag                                                                                                                                                  >  2:772577/1‑232 (MQ=255)                                                                                 ttcccGGCAAATAACCTTTTAAGCTTAATTTTAATGGTTATGAGCCGTTCCCAATTTCCTGGCGGTGAAAAC‑CATCCCCTGTTTCCATGGCATTTTTAACCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCa                                                                                                                                                <  2:17896/235‑1 (MQ=255)                                                                                 cccccGGCAAATTACCTTTTAAGCTTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGTTTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCa                                                                                                                                                <  2:757131/235‑1 (MQ=255)                                                                                      ggCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTTAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCAc                                                                                                                                          <  1:326015/239‑1 (MQ=255)                                                                                           aTAACCTTTTAAGCGTAATTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAAC‑CCTCCCAGATTTCCATGGCATTTTTATCCGTATGTGTT‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGct                                                                                                                                     <  1:491211/239‑1 (MQ=255)                                                                                               ccTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAA‑CATCCCCGATTTCCCTGGCATTTTTTTCCGTATGTGTG‑ACC‑A‑CAAAAGGTTGGGGAAAACCATCCATTTTAATTGGTTTTTTTCA‑CAGGTTTTTTTTGGGGGGGGCCCCCTTTTTTGGGTTGTTTTCTTTTGGGCAAAGGTTGG‑T‑G‑GTTTTTACCCACAGAAAGCAAACCAccaggccttc                                                                                                                                 >  2:725129/1‑230 (MQ=255)                                                                                                           aaTTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGc                                                                                                                      <  1:294539/238‑1 (MQ=255)                                                                                                                      tgatgaGCCTTTACCAATTTCATGGCGGTGAAAAC‑CAGCCACGATTTCCATGGCATTTTTATCCGTATGGGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGa                                                                                                          <  1:1054434/239‑1 (MQ=255)                                                                                                                              cTTTACCAATTTCAGGGCGGTGAAAAC‑CGTTCACGATTTCCATGGCATTTTTATTCGTATGTTTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATg                                                                                                  <  1:379280/239‑1 (MQ=255)                                                                                                                                   ccAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                             <  1:573894/130‑1 (MQ=255)                                                                                                                                        ttCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCt                                                                                                                     <  1:830094/210‑1 (MQ=255)                                                                                                                                                            tCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCt                                                                                                                                                                       >  2:57690/1‑143 (MQ=255)                                                                                                                                                            tCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCt                                                                                                                                                                       <  1:57690/143‑1 (MQ=255)                                                                                                                                                                    ttCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATCCGGATTTCATTAATCGATTTGAGGGAATCATCTTCCg                                                             >  2:379504/1‑239 (MQ=255)                                                                                                                                                                        aTGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGTTCATCTTCCGGCa                                                          >  1:247786/1‑238 (MQ=255)                                                                                                                                                                           gCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTACTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                                         >  1:112735/1‑236 (MQ=255)                                                                                                                                                                               ttttttCCGTATGTTTG‑ACC‑A‑CCATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAAGC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGGAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                                  <  1:1013704/239‑1 (MQ=255)                                                                                                                                                                               ttttATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑GTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATCCGGATTTCATTAATAGCTTTGATGCGATAATATTCCGCTAAGCTggtg                                                  >  1:1058178/1‑237 (MQ=255)                                                                                                                                                                                    tCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGaaa                                             <  1:799017/239‑1 (MQ=255)                                                                                                                                                                                       gTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAA‑GGGCTAACCTGATGGATCGCCATGCGGATTTCTTTAATCGCTTTGATGCGATCACTTTCGGGCAAGCTGGAGGGCAAttcc                                         >  2:985044/1‑237 (MQ=255)                                                                                                                                                                                        tATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTAACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGCTGCGATCATCTTCCGGCAAGCTGGCGAAAAATTGc                                         >  2:246505/1‑239 (MQ=255)                                                                                                                                                                                        tATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCAATGGCTCTTCACGAAAGGGGTTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAAtt                                           >  1:609450/1‑237 (MQ=255)                                                                                                                                                                                          tgtgtg‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGt                                       >  1:19164/1‑239 (MQ=255)                                                                                                                                                                                                 aCC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                                      >  2:452953/1‑202 (MQ=255)                                                                                                                                                                                                 aCC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                                      <  1:452953/202‑1 (MQ=255)                                                                                                                                                                                                 aCC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                                  >  1:316732/1‑238 (MQ=255)                                                                                                                                                                                                 aCC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCAATGCGGATTTCATTAATCGCTTTGAGGCGATCATCTTCCGGCAAGCTGGGGAGAAAATGGGTTaaaa                                  >  2:1034749/1‑237 (MQ=255)                                                                                                                                                                                                   c‑a‑caATTGTTTGGGTAAAAACATTTTTTTCAA‑TGGATTTTTGCGGAAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCt                               <  2:16907/239‑1 (MQ=255)                                                                                                                                                                                                              ttgggtAAAACCTTCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTGTGGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCTACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAAt                       <  1:367486/238‑1 (MQ=255)                                                                                                                                                                                                              ttgggtAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCAAGAAATTGCGTTAAATCCTGCGTTaa                        >  1:210650/1‑237 (MQ=255)                                                                                                                                                                                                               ttggtTAAACCTTTTTTTTTAA‑TTGATTTTTTGAGGAGTTTTTTTTCGGGGTGTGCCACGTTGGTTGGGTTGTAATC‑CTTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCg                     <  2:461901/237‑1 (MQ=255)                                                                                                                                                                                                               tgggtAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCg                     >  1:9463/1‑239 (MQ=255)                                                                                                                                                                                                                 ggtAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt                   >  1:161676/1‑239 (MQ=255)                                                                                                                                                                                                                  gtAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTGTGTCGGGTGGTGCACAGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAA‑GGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                  <  1:985044/238‑1 (MQ=255)                                                                                                                                                                                                                   tAAAACCTTCTATTTCAA‑TGGTGTTTTGCAGGAGTTTTTTGTCGGGTGGTGCCACGTTGGTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                  <  2:118491/238‑1 (MQ=255)                                                                                                                                                                                                                   tAAAACCAGTTGTTTCAA‑TGGATTTTTGCAGCAGTTTTTTTGCGGGTGGTGCCAAGTTGTTAGGGGTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAAGCGCTTTGATGCGATCATCTTCCGGCAAGCGGGCGAGAAATTGCGTTAAATCCTGCGTTAAGCgttg                  <  1:361992/238‑1 (MQ=255)                                                                                                                                                                                                                    aaaaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                <  1:725129/239‑1 (MQ=255)                                                                                                                                                                                                                     aaaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                <  1:195796/238‑1 (MQ=255)                                                                                                                                                                                                                      aaCCTTCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCt                                                                                                                     <  1:463576/136‑1 (MQ=255)                                                                                                                                                                                                                      aaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGCAAGGGGCTAACCTGATGGATCGCCATGCGGCTTTCATTAATCGCTTTGATGCGTTCATCTTCCGGCAAGCTGGCGGGAAATTGCGTTAAATCCTGCGTTAATCctttttca              >  2:244853/1‑237 (MQ=255)                                                                                                                                                                                                                      aaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGTTGCGATCATCTTCCGGCAACCTTGCGAGAAATTGCGTTAAATCCTTGGTTAATCgttgtggg              >  2:563450/1‑236 (MQ=255)                                                                                                                                                                                                                      aaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCAAGAAAGGGGCt                                                                                                                     >  2:463576/1‑136 (MQ=255)                                                                                                                                                                                                                      aaCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCCCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAATTCGGTTTAATCTTGCGTTAATCGTCGTTGc              >  1:196702/1‑239 (MQ=255)                                                                                                                                                                                                                       aCCATTTTTTTTAA‑TGGATTTTTGCAGGAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAAGC‑CTTAGGCATAAGTTGGAT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATTGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCa             <  2:409979/239‑1 (MQ=255)                                                                                                                                                                                                                       aCCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCa             <  1:749444/239‑1 (MQ=255)                                                                                                                                                                                                                         cgacTATTTCAA‑TTGATTCTTACATCATTTTTGTGGCGGTTGGGGCCACTTTGGTTGGGGTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑GTGTTAACCCACAGAACTCAATCCACTGGCTCTTCACGAAAGGGGCTAACCGGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTTCATa           <  1:777909/236‑1 (MQ=255)                                                                                                                                                                                                                           tCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCTTTGTTCCATAGt         >  1:370362/1‑239 (MQ=255)                                                                                                                                                                                                                           tCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATatt         >  1:419214/1‑237 (MQ=255)                                                                                                                                                                                                                           tCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGt         >  1:305016/1‑239 (MQ=255)                                                                                                                                                                                                                           tCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTAAAACGGTGTTGCAAAGt         >  2:297400/1‑239 (MQ=255)                                                                                                                                                                                                                                  cAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAATTGCGTTAAATCCTGCGTTAATCGTTGTGCAAAATTATccccc  >  2:514988/1‑239 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   GGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAC‑CATCCACGATTTCCATGGCATTTTTATCCGTATGTGTG‑ACC‑A‑CAATGGGTTGGGTAAAACCATCTATTTCAA‑TGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATC‑ATTAGGCATAAGTTGGCT‑G‑TTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCC  >  NC_000913/634915‑635374

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