Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,358	C→T	G252G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,358	0	C	T	97.0%	187.3 / NA	70	G252G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/37);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGAT  >  NC_000913/4620114‑4620600                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTGCAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtg                                                                                                                                                                                                                                              <  1:235200/251‑1 (MQ=255)  cGCTGACTCCGTGTTCCAGTTTGCTGGACACGAAGGAACCGTTGGCCTGGATAAACACTACGAAATGTGCGAAATCGCCCGTGAAGAGCGGCCCAACGGCGGCTACAACATCGGTCGTGTTATCGCTGGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtgg                                                                                                                                                                                                                                             <  2:172989/251‑1 (MQ=255)    cgtactccgtttgtcAGATTGCCTGCCCAGAAGGAACCTTTGGCCTGGATAAACTCTACGTACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGtt                                                                                                                                                                                                                                           <  2:564585/237‑1 (MQ=255)     tGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGATATCGCTCGTCCGTTTATCGGCGACAAAGCCGGACACTTCCAGCGTACCGGGAAACGGCACGACCTGGCTGGTTG‑GCCGCCAGCCCCGACACGGCCTAC‑GAAACCTGTTGATGAAAAACACGGGCAggagggtt                                                                                                                                                                                                                                          >  2:332042/1‑245 (MQ=255)          ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGTTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                                                                                                                                                                     <  1:587606/251‑1 (MQ=255)               ttCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGGGCCGCCAGCAACGAC‑CGTGCTGCAGAAACCGGTTGATGAAAAAAACGGTAAGGGGGGTTCTGTCGGTa                                                                                                                                                                                                                                 >  1:371230/1‑250 (MQ=255)                 ccAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTTTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaaa                                                                                                                                                                                                                               <  2:94773/250‑1 (MQ=255)                   aGATTGCCTGCCACGAAGGAACCTTGGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACGTCCAGCGTACCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAtt                                                                                                                                                                                                                            <  1:319198/251‑1 (MQ=255)                        gcctgccACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGcaga                                                                                                                                                                                                                       <  1:626593/251‑4 (MQ=255)                                      aacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGTCGGCTACAACATCGGTCAAGTTATCGCTCGACCGTTTATCGGAGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGATGCATAAACTGGTTGATGAAAAAAAAGGACAGGGGGTTTCTGTCGGTAAAATTGAAGACATATACGCCAACt                                                                                                                                                                                                         >  1:397843/11‑251 (MQ=255)                                          tttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAAGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGGGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCACACTGCgg                                                                                                                                                                                                     >  1:208218/7‑251 (MQ=255)                                            tggccTGGATAAACCCTACGAACTGGGCGAAATCGCCCGTGTAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGt                                                                                                                                                                                                    <  2:345062/246‑1 (MQ=255)                                              gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGACGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCACCACCGCC‑CGTGCTGCAGAAACTGGTTGATGGAAAACACGGTCAGGTGGTTTCTGTAGGTAAAATTTCAGACATCTAAGACAACTCCGGtgtc                                                                                                                                                                                                 >  1:105480/3‑248 (MQ=255)                                               ccTGGATTAACTCTCCGAACTGTGCGAAATCGCGCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGATAAAGCCGGTAATGTCCAGCGTACCGGTAACCGTCACGACCTGGTTG‑TTGAGCCGCCATCACCGAC‑CGTGCTTCAGAAACTGGTTGATGAAAAACACTGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                                                                                                                                <  1:617375/250‑1 (MQ=255)                                               ccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                                                                                                                                <  1:446750/250‑1 (MQ=255)                                                   gtaaaaccctccGAACTGTTCGAAACTGCCAGTGAAGTGCTGCCCAACGGCGGCTACAACTTCGTCCGTTTGATCGCTCGTCCGTTTATCGGTTCCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                                                                                                                            <  2:172842/240‑1 (MQ=255)                                                   gATAAACTCTACGAACTGGGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                                                                                                                            <  1:501122/251‑1 (MQ=255)                                                      aaaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaaga                                                                                                                                                                                          >  1:513447/1‑248 (MQ=255)                                                                    gtgCGAAATCGCCCGTGAAGAGCTGTCCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTg                                                                                                                                                                           <  2:363562/251‑1 (MQ=255)                                                                    gtgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACt                                                                                                                                                                            <  1:231169/250‑1 (MQ=255)                                                                             cGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAc                                                                                                                                                                   <  1:168044/250‑1 (MQ=255)                                                                              gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                                                                                                                                                                 >  1:156971/1‑251 (MQ=255)                                                                                   tGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                                                                                                                                                                >  1:235136/1‑183 (MQ=255)                                                                                   ggAAGAGCTGACCAACGGCGGCTACAACATCTGTCGTGTTATCGCTCGTCCGTTTACCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                                                                                                                                                                <  2:235136/182‑1 (MQ=255)                                                                                     aaGAGCTGACCAACGGCGGCTACAACATCGGTCGGGTTATCGCTCGTCTGTGTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGTGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                                                                                                                                                          <  2:388344/251‑1 (MQ=255)                                                                                           tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                    >  1:386476/1‑251 (MQ=255)                                                                                                aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                               <  1:294656/251‑1 (MQ=255)                                                                                                     cggcTACAACATCGGCCGTGTTATCGCTCGCCCGTTTATCGGCGACAAAGCCGGTAAATTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCCCCGAC‑AGCGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTACCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAAGCCACCAGCaa                                                                                                                                            >  2:293571/1‑249 (MQ=255)                                                                                                       gcTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                          >  1:315824/1‑249 (MQ=255)                                                                                                        cTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                                                                                                                          >  1:34968/1‑248 (MQ=255)                                                                                                            aacatCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTAACGGTAACCGTCACGAACTGGCTG‑TTGAGCCGCCAGCACCGCC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGATTCTGTCGGCAAACTTGCAGACATCTACGCCACCTGCGGCATCACCACGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACACCACCATCAAAGAGCTGa                                                                                                                                    >  1:286307/4‑250 (MQ=255)                                                                                                                 cgggcgttttATCGCGCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGACACGTCCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAaagaag                                                                                                                              <  2:46657/243‑1 (MQ=255)                                                                                                                         ttATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGAt                                                                                                                      >  2:167342/1‑251 (MQ=255)                                                                                                                           atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                                                                                                                      >  1:343004/5‑169 (MQ=255)                                                                                                                           atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                                                                                                                      <  2:343004/165‑1 (MQ=255)                                                                                                                           agcGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCAGCTAGCACCGAC‑CGTGCTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTTGTTTCTTTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAATAAAGTGAAAGCAACTGGCCTGGACGCGTTTTTTGACGCCTCCATCAAAGAGATGAAAGAAGCGGGTGATaa                                                                                                                    <  1:605667/249‑1 (MQ=255)                                                                                                                             cgcctgtcctttaaacggtgACAAAGCCGGTAACTTCCAGCTTACTGGTAACCATCACGCCCTGGCTG‑TTGATCCGCCAGCACCGAT‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCATGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACCTCGGTATCACCCAGAAAGTGAAAGTAACTGGCCTGGACGCGCTGTTTGACGCCACCTTCAAAGAGATGAAATAAGCCGGTGATatca                                                                                                                  <  2:78078/232‑4 (MQ=255)                                                                                                                             cGCTCGTCCGTTTATCGGTGACAAAGCCGTTAACTTCCAGCGTACCTGTAACCTTCACGACCTGGCTT‑GTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCAAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGTGATGAAAGAAGCGGGTGATAac                                                                                                                   <  1:400470/250‑1 (MQ=255)                                                                                                                             cGCTCGTCCGTTTACCGGTATCAAAGCCGGTAGCTTCCAGCGTGCCGGAAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGGAACTGGTTGATGAAAAACACGGTCAGGTGTTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTTGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAaca                                                                                                                  <  1:285304/251‑1 (MQ=255)                                                                                                                                 cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAt                                                                                                              >  1:562635/1‑251 (MQ=255)                                                                                                                                   tCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCg                                                                                                            >  1:591204/1‑251 (MQ=255)                                                                                                                                       tttATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCATACATCTACGCCAACTGCGGCATCACCAAGAAAGAGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCCCCATCAAAGAGATGAAAGAATCggtgaataaccccaacgcctt                                                                                                        >  2:256456/1‑232 (MQ=255)                                                                                                                                         tATCGGTGAAAAAGCCGGAAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTc                                                                                                       <  2:156971/250‑1 (MQ=255)                                                                                                                                            cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                    >  1:387987/5‑250 (MQ=255)                                                                                                                                            cggtgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCAcc                                                                                                    >  1:573123/5‑250 (MQ=255)                                                                                                                                                       ccGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACACCACCCTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGAc                                                                                        >  2:261468/1‑251 (MQ=255)                                                                                                                                                        cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTAATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACt                                                                                       <  1:169618/251‑1 (MQ=255)                                                                                                                                                         ggTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACtt                                                                                      >  1:267888/1‑251 (MQ=255)                                                                                                                                                          gTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTc                                                                                     <  1:48938/251‑1 (MQ=255)                                                                                                                                                                cccctcgtcccGGTAACCGTCACGACCTGGCTG‑TGGAGCCGCCAGCACCGTC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGGCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGActct                                                                               <  2:387987/242‑1 (MQ=255)                                                                                                                                                                  catcGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCACCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGCTGAAAAACACGATCAGGAGGCTTCTGTCGGTAACAATGCAGACATCTACGCCAACTGCGGTATCACCACGAAAGGGAAAGCAACTGGCCTGTACGCGCTGTTTGACGCCACCATCAAAGAGATACAAGAAGCTCGAAATAACACCCACGACTTCCCCCAATTCGTTGACTTCGActctac                                                                             >  2:60062/4‑249 (MQ=255)                                                                                                                                                                     cGTACCGGTAACCGTCACGACCGGGCTG‑TTGAGCCGCCAGCACCGAC‑AGTGCAGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGGCGGAAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTg                                                                          <  2:265870/251‑1 (MQ=255)                                                                                                                                                                       tACCGATAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGCAGCGGGTGATAACACCATCGTCTTCACCAACATCGTTGACTTCGCCTCTTCCTggg                                                                        >  2:311708/1‑251 (MQ=255)                                                                                                                                                                                       gACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGt                                                        >  1:546951/1‑251 (MQ=255)                                                                                                                                                                                         ccTGGCTG‑TTGAGCCGCCAGCACCGAA‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCCGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACATCg                                                      >  2:457009/1‑251 (MQ=255)                                                                                                                                                                                                   tGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                      <  1:640468/130‑1 (MQ=255)                                                                                                                                                                                                   tGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTg                                                                                                                                                                      >  2:640468/1‑130 (MQ=255)                                                                                                                                                                                                             gCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                   <  1:399864/251‑1 (MQ=255)                                                                                                                                                                                                               accgac‑cgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                      >  1:189881/1‑246 (MQ=255)                                                                                                                                                                                                               accgac‑cgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                 <  1:251500/251‑1 (MQ=255)                                                                                                                                                                                                               accgac‑cgTGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                      <  2:189881/246‑1 (MQ=255)                                                                                                                                                                                                                      cgTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGAGATACCACCATCGTCTTCACCAACTTCGTTGACTACGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGagcttt                           >  2:505988/1‑246 (MQ=255)                                                                                                                                                                                                                            gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                      <  1:256456/250‑1 (MQ=255)                                                                                                                                                                                                                            gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                      <  1:440735/250‑1 (MQ=255)                                                                                                                                                                                                                                     ggTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTg             <  1:277375/250‑1 (MQ=255)                                                                                                                                                                                                                                      gTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGcc           <  2:90776/251‑1 (MQ=255)                                                                                                                                                                                                                                         gatgaAAAACACGGACAGGTGGATTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACACTTCCTGGGGCCACCGTCGCTCCGTCGCCGGTTATGCTGCGGGCCTGGAGTTGTACCACCGCCGCCTGCCggg        >  2:276264/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTg    >  1:408635/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCTGGTGATAACACCATCGTCTTCACCAACGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGa   <  1:98925/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGACAGTGAAAGCAACTGGCCTGGCCGAGCTGTTTGACGCCACCAACAAAGAGATGAAAGAAGCGGGTGTTAACACCATCGTCTTCCCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCAGCCGCCGTCTGCCGGAGCTGAt  >  2:601797/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGAT  >  NC_000913/4620114‑4620600

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