Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,111,259	Δ1 bp	coding (718/900 nt)	rfbD ←	dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,111,259	0	G	.	95.1%	153.5 / 2.6	42	coding (718/900 nt)	rfbD	dTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (21/18);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTTC  >  NC_000913/2111021‑2111507                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         tcatcaTTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTg                                                                                                                                                                                                                                             >  1:18365/1‑251 (MQ=255)   atcatcacTTTTAATTCATCTTGCTCCACCAGCACGAACAAGATGCAAAAACTATTAAAGTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGACAACCTTCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTTAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGAT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGTG‑AATGCCTGCTTTgcg                                                                                                                                                                                                                                           <  2:16831/244‑1 (MQ=255)   atcatTCCTTTTAATTCATATTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGGGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTATGTG‑AAAGGGAATGCCTGCTTTg                                                                                                                                                                                                                                             >  2:221540/1‑249 (MQ=255)     aattacgtttaagtaaTCTTGCGACACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAAGAATTCATGGAGCAGTCGTGTCCCTCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCAGTAGTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGCTACTGCGGTGAGCTGGTTGAGTGCAAGGG‑AATGCCTTCTTTgcgcg                                                                                                                                                                                                                                         <  2:61508/236‑1 (MQ=255)         ccTTTTAATTCATCTGGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGGAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCctc                                                                                                                                                                                                                                     <  2:458769/251‑1 (MQ=255)               aaTTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGGAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCaa                                                                                                                                                                                                                                <  2:156965/250‑1 (MQ=255)                          cTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAgcga                                                                                                                                                                                                                      >  2:456046/1‑248 (MQ=255)                            ccaccaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgcagca                                                                                                                                                                                                                  <  1:311718/251‑1 (MQ=255)                                     cGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTAc                                                                                                                                                                                                         >  1:546287/1‑251 (MQ=255)                                       aaCAAGATGCACAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTAcc                                                                                                                                                                                                        >  2:336218/1‑250 (MQ=255)                                               gCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCCACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTCTTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTAATGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCGGGTTGTgg                                                                                                                                                                                               >  2:513122/1‑251 (MQ=255)                                                      cTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGCGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACt                                                                                                                                                                                        >  2:129526/1‑251 (MQ=255)                                                       tATTAAATAGCTGTATTCGTAAATAATTCATTGCGCATTCGTTTCACGCCAAACTGCCAGTCAGGCAAAACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGGGCTGGTGTAGGATAGCCT‑GTTGTTTGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑ACTGCCTGCTTTGCACGCCTCTTCAAAACCCAGCTAAGCATAATCTTACCAGGTTGTGGTACCCCTg                                                                                                                                                                                       >  2:257659/1‑251 (MQ=255)                                                                  tGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATg                                                                                                                                                                            >  1:130055/1‑251 (MQ=255)                                                                      gtcgtcAATAATTCATTGAGCATTCGTTTAACGCCCACCTGCCAGTCAGGCAAGACACGCGCAAAGTTCTGCTGAAATTTTTCTGTATTCCGGCGAGAGTTATCTGG‑ACGCCGAGCTGGTGTAGGATACGCT‑GTTGTTGGTACTGCGTTGCGCTTGTTGCGCG‑AAAGGGCATGCCTGCTTTGCGCGCATATTCAAAAAACGGCGCAGCATCATCATACCAGGTTGTGGTACCAACGGCTACCACATGGTAc                                                                                                                                                                        >  2:504493/1‑251 (MQ=255)                                                                                        aGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑CATGCCTGCTTTGCGCGCCTCTTCCAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAACTGGTACAAGCCTGCAACATCCGGt                                                                                                                                                      >  2:283184/1‑251 (MQ=255)                                                                                            ttCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑ACTGCCTGCTTTGCGCGCCTCTTCAAAAAACAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTTGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTAt                                                                                                                                                  >  2:380322/1‑251 (MQ=255)                                                                                               ggtacACGCCAACCTGCCAGTCAGGCAAGACAAGCGAAAAGTTCTGCTGAAATGTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGTTTTATTCa                                                                                                                                               <  2:263300/247‑1 (MQ=255)                                                                                                        cAACCTGCCGGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGGGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGCTAGGCT‑GTTGTTGGTACTGCGGTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAGCCAGCGCAGCCTAATc                                                                                                                                                                                                             >  1:359048/1‑180 (MQ=255)                                                                                                                                   cAAAGTTCTGGTGAAATTTTTTTGTATTAAGGCGGAAGTTA‑GTGGGACGACGAGCTGGTGTAGGATAGGCT‑GTGGTTGGGTCTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTGGCGGGCCGCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCAATGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATcagcca                                                                                                           <  2:269869/251‑1 (MQ=255)                                                                                                                                    aaaGTTCTGCTGAAATTTTTCTGTATTAAGGCGCGAGTTATGTGG‑CCGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGc                                                                                                                                                                >  2:250466/1‑197 (MQ=255)                                                                                                                                    aaaGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGc                                                                                                                                                                <  1:250466/197‑1 (MQ=255)                                                                                                                                            gctgAAATTTTTTTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGGGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTTCCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCa                                                                                                   <  1:110209/250‑1 (MQ=255)                                                                                                                                                aaaTTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCAcc                                                                                              >  1:283360/1‑251 (MQ=255)                                                                                                                                                 aaTTTTTCGGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTGGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTTCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCt                                                                                             <  2:510385/251‑1 (MQ=255)                                                                                                                                                  aTTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAAACAGCGCAGCATAATCGTACCAGGTTGTGGAACCACTGGCTACCAACTGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGCATTGCATGGGCCTTACAATCAGCCAACAGTTAAGCACCTg                                                                                            >  2:154771/1‑251 (MQ=255)                                                                                                                                                     tttCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGtt                                                                                          <  2:18365/250‑1 (MQ=255)                                                                                                                                                        ttgttgtaaggcaggtgTTATGTGG‑ACGGCGTTGTGGTGTACGATAGGCT‑GGTGTTGTTACGGGGTTGATCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCGTCAAAAACCAGCGGAGCATAATCGTACCAGGTTGTGGTAACACTGGCTACCAAATTGTGCAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGATGTACAATCAGCCAGCAGTTCAGCACCTGTTGgc                                                                                       <  2:8939/234‑1 (MQ=255)                                                                                                                                                            aatAAGGCGTGCGTTTTGTGG‑ACGACGAGATGGTGTATGATAGGCTTGTGGTGGGT‑CTGCGTTGAGCTTGTTGAGCG‑AAAGGAAATGCCTGCGTTGCGCGCCACTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTTCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTGTATTCAGTGCGACATGAATGGCATGTGCTGTACAATCATCCAGCAGTTCAGCACCTGTTGGCGCAc                                                                                   <  2:117848/248‑1 (MQ=255)                                                                                                                                                                     agagTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCCAACTGATc                                                                         >  2:371010/1‑251 (MQ=255)                                                                                                                                                                                     acgaGCTGTGGTAGGATAGGCG‑GTTGTGGGTACGGTGCTTATCTTGGTGAGTGCAAGGG‑AAGGCCTGTCTTGCTCGCCTCTGCAAAAACCAGCGCAGCAGAATCGGACCAGGTTGTGGTTCCTCTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTTCGTCACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCTCCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTa                                                           <  2:384725/250‑1 (MQ=255)                                                                                                                                                                                       ggggtggggtAGGATTGGCC‑GGTGTTGGTACTGGGTTGGGGTTGTTGAGTGCCAGGG‑AATTGCTGCTTTGCGCGCCCTTTCAAAAACCCGCCCAGCATAATCCTTCCAGGGTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAtt                                                        <  1:186728/243‑1 (MQ=255)                                                                                                                                                                                                            gttgttGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCa                                      <  1:290695/249‑1 (MQ=255)                                                                                                                                                                                                             ttgttgGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGAc                                   >  2:410972/1‑251 (MQ=255)                                                                                                                                                                                                                   gTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTCATACCCGCTAATTCTTCCCGCTCTTTTGCCAGCCGTAAc                              >  2:449659/1‑250 (MQ=255)                                                                                                                                                                                                                    tACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACAt                            >  1:243996/1‑251 (MQ=255)                                                                                                                                                                                                                            tGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTg                     <  1:456046/250‑1 (MQ=255)                                                                                                                                                                                                                               gCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGa                 <  1:318394/251‑1 (MQ=255)                                                                                                                                                                                                                                   gttgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCCGTTCAGCACCTGTTGGCGCACCAAACTGATCGCTAATAACCGCTAATTCTTCACGCCCTTTTGCCAGACGTAACATCGTTTTGGCGAAGtt             >  1:514088/1‑251 (MQ=255)                                                                                                                                                                                                                                     tgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTa            <  1:231797/250‑1 (MQ=255)                                                                                                                                                                                                                                     tgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCAAGACGTAACATCGtttttgagaaggtat           >  2:115442/1‑240 (MQ=255)                                                                                                                                                                                                                                       aGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAttt         >  1:1005/1‑251 (MQ=255)                                                                                                                                                                                                                                               gggATTGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACTAGGTTGTGGTTCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttccttttc  <  1:301772/251‑1 (MQ=255)                                                                                                                                                                                                                                               gggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTCCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttccttttc  <  2:203752/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGG‑ACGACGAGCTGGTGTAGGATAGGCT‑GTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTTC  >  NC_000913/2111021‑2111507

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