Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,616,022	T→C	V348V (GTA→GTG)	yjjI ←	DUF3029 domain‑containing protein YjjI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,616,022	0	T	C	96.4%	160.0 / ‑2.5	56	V348V (GTA→GTG)	yjjI	DUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (28/26);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAG  >  NC_000913/4615777‑4616265                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTc                                                                                                                                                                                                                                                 >  1:167390/1‑251 (MQ=255)              cgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGct                                                                                                                                                                                                                                    >  1:129894/1‑251 (MQ=255)                       gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCg                                                                                                                                                                                                                           >  1:333292/1‑251 (MQ=255)                                    gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCa                                                                                                                                                                                                               <  1:618890/250‑1 (MQ=255)                                                     cATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacac                                                                                                                                                                                             >  1:280213/1‑251 (MQ=255)                                                         atgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGc                                                                                                                                                                                          >  1:494443/1‑250 (MQ=255)                                                           gCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAg                                                                                                                                                                                       <  1:454002/251‑1 (MQ=255)                                                             ggggCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGtt                                                                                                                                                                                     >  1:633029/1‑251 (MQ=255)                                                               ggCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTaa                                                                                                                                                                                   <  1:155359/251‑1 (MQ=255)                                                                  gACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCg                                                                                                                                                                                >  1:64576/1‑251 (MQ=255)                                                                      gTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                                                                                                                                                            <  1:544542/251‑1 (MQ=255)                                                                           cAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCa                                                                                                                                                                       <  1:12685/251‑1 (MQ=255)                                                                              aTGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGAAATCGGAACTGATCCCCGACTGTGCGTGTAACATGCCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGGTAACCCACTTCATTTGCGGCGGCTTCTTTACCGTATCGCGGGGCAATCCCTTTTTTCTCACACAGAAAGTTAAACGCATCCGCCAgac                                                                                                                                                                    >  2:312075/1‑249 (MQ=255)                                                                                   tGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑Ac                                                                                                                                                              >  1:373890/1‑251 (MQ=255)                                                                                    gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑Ac                                                                                                                                                              <  1:116855/250‑1 (MQ=255)                                                                                    gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACa                                                                                                                                                             <  1:393804/251‑1 (MQ=255)                                                                                       cGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGCACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCCCGGGGGGATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGc                                                                                                                                                          >  1:26758/1‑251 (MQ=255)                                                                                                    ccGCCATACGGCAAACGGGGTCCCGGGGTGGTGCCGATATCGGACTTGTTCCCGGACTGGGGGGGTAACTTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGcac                                                                                                                                             <  2:117295/250‑1 (MQ=255)                                                                                                      gCCATACGGCAAACGCGCGCCCGGCGTGGAGCCGATATCGGAACTGAACCCCGACTGTGCGTGTAACCTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTAGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCCTAGCGCGGGGCAATACATTCTTTCTCAAAAAGCAATTTAACCCCTTCCGCCCTCCCATCACTTGC‑AAACATTGGcaac                                                                                                                                            >  2:591484/1‑248 (MQ=255)                                                                                                        cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑CCATGCCAAACATTGGCACaaa                                                                                                                                          >  1:393597/1‑250 (MQ=255)                                                                                                        cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCCAACATTGGCACaaa                                                                                                                                          >  1:35307/1‑250 (MQ=255)                                                                                                              gCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACACACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGACAATCCCTTCTTTCTCAAACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACAGTGGCCCAACACgatc                                                                                                                                    >  1:399293/1‑247 (MQ=255)                                                                                                                        gcCCGGCGGGGTGCCGATATCGGACTTGACCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTTGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCa                                                                                                                         <  2:483957/251‑1 (MQ=255)                                                                                                                         cccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATc                                                                                                                          >  1:259790/1‑249 (MQ=255)                                                                                                                          ccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGc                                                                                                                       <  1:503573/251‑1 (MQ=255)                                                                                                                             gCGTGGGGGCGATATCGGACCTGAGCCCCGACTGGGCGTGTAACATGGCGAGTGTGTGCCATCCATATTGCGCGGGGGTATTGGCGACAAACGCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCGTCTTTACCGTAGCGCGCGGCAATCCCTTCATTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAG‑ACATGCCAAACATTGGCACAAAGCGTTCAGGGTTAATCAGCCct                                                                                                                    <  2:107184/251‑1 (MQ=255)                                                                                                                              cGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTAACGGGGGAATTGGCGCCAGACCCCGCCAGTTGCGCGCTGCTGCGATAAACCGCTTCATTTGCGGCGGCTTATTTAACGGAGAGCGCGGCAATCCCTTAATTCGCCCACAGCCAGTTCACCGCTTCCGGCCGCCCAT‑ACATGCCAAACATTGGAACAAAACATTCAAGGTTAATCAAcccct                                                                                                                   >  2:520898/1‑249 (MQ=255)                                                                                                                               ggggtgccgagagtgaaaCTGACCCCCGACTGGGCGTGGAACATGGCGCGTTGTTGCCGGCCATAGTTCACGGGGGTATTGGCGTCAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCAATTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttc                                                                                                                  <  2:44725/235‑1 (MQ=255)                                                                                                                                gggtgCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttct                                                                                                                 <  1:166918/250‑1 (MQ=255)                                                                                                                                   tgCCGATATCGGAACTGATACCCGACTGTGCGTGTAACATGGCGCGTTTTTACCAGCCATATTTCACGGGGGTATTGGCGACAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGACTTCTTTACCGAAGCGCGCGGCAAACCCTTATTTATCACACAGAAAGTTAACCGCATCCGCCCGACCAT‑ACCTGCCAAACATTGGCACAAAACTGTCAGGGGTAATCAGcactacttac                                                                                                              >  2:568712/1‑242 (MQ=255)                                                                                                                                    cccGATATCGGAACCGATCCCCGTCTTGGCGTGTAACATGGCGCGTTTTTTCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGTCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCa                                                                                                             <  1:40706/250‑1 (MQ=255)                                                                                                                                                     tCCCCGACTGTGCGTGTCACATGGCGCGTTTTTGCCAGCCATATTTCCCGGGGGTCTTGGCGACAAACTCCGCCAGTTGCGCGATGATGCGCTAACCCACtt                                                                                                                                                                                                                                                  >  2:287513/1‑102 (MQ=255)                                                                                                                                                     tCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACtt                                                                                                                                                                                                                                                  <  1:287513/102‑1 (MQ=255)                                                                                                                                                                 cGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGCTAACCCACTTCATTTGCGGCGGCTTCTTTACCGAAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACCTGCCAAACATTGGCACAAAACGTTACGGGTTCATCCGCCCTTCTTTCACCAGGAAGATATTCTCAAAGAAGGGTGAt                                                                                >  2:284419/1‑251 (MQ=255)                                                                                                                                                                  gtgtAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTTGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCATGAAGCTATTCTCAAAGAAGTTTGAtt                                                                               >  2:651683/1‑251 (MQ=255)                                                                                                                                                                      accagggcccgTTTTTCCCAGCCATATTTCACGGGGGTATTGGCGACAATCTCCGCCAGTTGCGCGCTGATGCGATAACCCACGTCATGTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAttgttg                                                                           <  2:360620/242‑1 (MQ=255)                                                                                                                                                                         aTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATa                                                                        <  1:549201/251‑1 (MQ=255)                                                                                                                                                                                             atatTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCAt                                                    <  1:137667/251‑1 (MQ=255)                                                                                                                                                                                               atTTCCCGGGGGTATTGGCGACAAACTCCGCCAGTGGCGCGCGGATGCGATAACCCACTTCATTGGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCg                                                  <  2:590987/251‑1 (MQ=255)                                                                                                                                                                                                  tCACGTGGGTATTGGCGCAAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAAGCCATGCTTTCACACACAGCAAGGTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCGTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCgat                                                <  2:232451/250‑1 (MQ=255)                                                                                                                                                                                                      gggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATgg                                           >  1:301759/1‑251 (MQ=255)                                                                                                                                                                                                       ggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGc                                          <  1:651683/251‑1 (MQ=255)                                                                                                                                                                                                              tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGatgatg                                            <  1:440842/242‑1 (MQ=255)                                                                                                                                                                                                              tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGAATCGatgatg                                            >  2:440842/1‑242 (MQ=255)                                                                                                                                                                                                                  cacaaacatcGCCAGTTGCTCTCTGATGCTTTAACCCCTTTCATTTGCGGCGCCTTCTTTACCGTACATCGCGGCAAACTCTTCTTTTTCTCTCGTCAAGGTTACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGTTTAATCAGCCCTGCTGTTACCAGGAAGCTATTCTTAAAGAAGTGTGATTGTTGCTAGTGGAATTCTCACCGCGCAGCGATGATGGCGAGctgttgct                               <  2:595524/242‑1 (MQ=255)                                                                                                                                                                                                                  atcAAACTCCCCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTGCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATGCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                               <  2:637630/249‑1 (MQ=255)                                                                                                                                                                                                                   aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                              >  1:300515/1‑251 (MQ=255)                                                                                                                                                                                                                   aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                              >  1:91505/1‑251 (MQ=255)                                                                                                                                                                                                                      acctacGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGACGGCTTCTTTCCCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACACTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCCCCAGGAAGCTATTCTCACAGAAGTGTGATTGTTGATAGATGGATTAACACCGCGCATCGATGATGGCAAGCTGCTGCAGGcc                           >  2:526363/6‑250 (MQ=255)                                                                                                                                                                                                                      aaCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                               >  1:637630/1‑247 (MQ=255)                                                                                                                                                                                                                       aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGc                            >  1:589550/1‑249 (MQ=255)                                                                                                                                                                                                                       aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCag                          <  1:532155/251‑1 (MQ=255)                                                                                                                                                                                                                        cTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGAGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGAACAACACGTTCAAGGTTAATCAACCCTTCTTTCACCCGGAAGCTATTCTCACAGCAGTGTGATTGTTGCTAGGGGAATTCACACCGCGCATCGATGATGGCGTTCGTCTGCTGGccgg                         >  2:352538/1‑248 (MQ=255)                                                                                                                                                                                                                         tCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                         >  1:577185/1‑250 (MQ=255)                                                                                                                                                                                                                            gCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTg                     <  2:333292/251‑1 (MQ=255)                                                                                                                                                                                                                                               gATAACCCACTTCATTTGCGGCTGTTTCTTTACCGTAGCGCGGGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATGTGCACAAAACGTTCAGGGTTAAGCAGCCCTGCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGTCAGTAGTGCGGTAGAGTGCGCGTAaag  <  2:393597/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     TTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAG  >  NC_000913/4615777‑4616265

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