Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,637	C→T	T345T (ACC→ACT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,637	0	C	T	98.2%	163.7 / ‑4.3	55	T345T (ACC→ACT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/36);  total (19/36)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAA  >  NC_000913/4620394‑4620879                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cGCCAAGTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACTATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGAcc                                                                                                                                                                                                                                               <  2:353786/251‑1 (MQ=255)   ccAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCAc                                                                                                                                                                                                                                             <  1:273171/251‑1 (MQ=255)      aCTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTGGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGt                                                                                                                                                                                                                                          <  1:335448/251‑1 (MQ=255)       cTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGtt                                                                                                                                                                                                                                         >  1:258980/1‑251 (MQ=255)            gtctCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGACTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGGCCACGGTTCTGCTGCGGGTCTGGCGCTGTTCGGCCGCCGTCTGCCGGAGCTGCTGTCTCTTCTGCGCGATTGCGCCATCCTGAACCTCAATGATGACCACGGTTGTGA‑t                                                                                                                                                                                                                                   >  2:342077/4‑251 (MQ=255)              aTCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑Tc                                                                                                                                                                                                                                  <  2:129709/250‑1 (MQ=255)                  ccAAGAAAGTGAAAGCAACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGAGAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCGGGTCGACCGCCGTCTGCCGGAGCGGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGAcc                                                                                                                                                                                                                             <  2:100455/251‑1 (MQ=255)                         agtgaaagcacacggCCTGGACGAGCTGTTTTACGCAACCATCAAAGAGATGAAAGAAGCGGGTGAGAACACCATCGTCTTCACCAACTTCGTTGACGTCGACTCTTCCTGGGGCCACCGTCGCGTCGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCg                                                                                                                                                                                                                      <  2:15262/238‑1 (MQ=255)                              aatcaaCTGGCCTGGACGCGCTGTTTGACGCCTCCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTTCCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGTGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACt                                                                                                                                                                                                                 <  1:55017/246‑1 (MQ=255)                                  aacgggCCGGGACTCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGTGTGATAACACCCTCGTCTTCACCATCGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGAcc                                                                                                                                                                                                             <  2:313040/247‑1 (MQ=255)                                   aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACca                                                                                                                                                                                                            <  1:420027/251‑1 (MQ=255)                                         cTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGCCGACATCCTGCTCCTCAAAGCTGACCAAGGATG‑GAGTCCGACATGGACCGGCACTGGCCACAcgcg                                                                                                                                                                                                      >  2:105403/1‑251 (MQ=255)                                                    gTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTc                                                                                                                                                                                           <  1:136298/251‑1 (MQ=255)                                                       tgccGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACCCCATCGTCTGCCCCTACGTCGTTGACTTCGACTCTTCCTGTGGCCTCCGTCGCGACGTCGCCGTTTATGCTGCGGGTCTGGAGCTGGTCGACCGCCGGCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCg                                                                                                                                                                                         <  2:86987/247‑1 (MQ=255)                                                            cccccctCAAATAGATGAAAGAAGCGGGGGATAACACCAGTGTCGTCACCAACCTCGTTGACATCGACTCTTCCGGGGGCTACCGTCGCGGCGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGTCCGGTACTGACCACACGCGTGAacagatcggaagatct                                                                                                                                                                                   <  2:136303/245‑14 (MQ=255)                                                             caccaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGCTGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGC‑TCCGACCTGGACCGGTACTGCCCACACGCGTGAACACATTACGGTACTgg                                                                                                                                                                                  >  1:26583/1‑251 (MQ=255)                                                                   cAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttac                                                                                                                                                                             <  1:327556/250‑5 (MQ=255)                                                                     aaGAGATGAAAGAAGCGGGTGATAACACCTTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttacggc                                                                                                                                                                          <  1:332951/251‑8 (MQ=255)                                                                     aaGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTACCGGAGCTGATGTCTCTGCTGCGCGATGGAGACATCCTGATCCTCACTGCTGACACCGGTTGTGA‑TCCGACCTGGACCGGTCCTGACCACACGCGTGAACAAAGTCAGGTACTGGtttacggc                                                                                                                                                                          >  2:459572/1‑244 (MQ=255)                                                                      agagATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGcc                                                                                                                                                                         >  1:436197/1‑251 (MQ=255)                                                                              aagaagCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTa                                                                                                                                                                 <  2:417186/251‑1 (MQ=255)                                                                                 aagCGGGTTATAACCCCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTaaaa                                                                                                                                                              <  1:321383/251‑1 (MQ=255)                                                                                         gctAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACCTTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGcac                                                                                                                                                      >  2:46224/3‑249 (MQ=255)                                                                                                ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGt                                                                                                                                              <  1:312130/251‑1 (MQ=255)                                                                                                  aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTcc                                                                                                                                            >  2:79883/1‑250 (MQ=255)                                                                                                    cGTCTTCACCAATGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCacc                                                                                                                                          <  1:222263/251‑3 (MQ=255)                                                                                                          cACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGaa                                                                                                                                     <  1:346761/250‑1 (MQ=255)                                                                                                          cACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGAAGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGGCCACACGCGTGACCACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGC‑CACCTGGGTCacagtgca                                                                                                                                     >  2:163195/1‑243 (MQ=255)                                                                                                               cctgcGTTGGCTTCGACTCTGCCTGGGGCCACCGTCGCGGCGTCGCCGGTTAGGCTGCGGGTCTGGGGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTTCTGACCACACGCGTGAACACATTCCGGTACTGGGTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGTGTCACCGTGAAACCTTc                                                                                                                               <  2:26583/247‑1 (MQ=255)                                                                                                                   cGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTcgcg                                                                                                                            <  1:135919/250‑1 (MQ=255)                                                                                                                   cGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCgg                                                                                                                           <  1:242039/251‑1 (MQ=255)                                                                                                                         cTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGATCACGGTTGCGA‑TCCGACCTGGAACGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAACGTAAAACCGGGATCAC‑TGGTTCATCGTGAAACCTTCGCGGATATCg                                                                                                                     >  2:210759/1‑251 (MQ=255)                                                                                                                            cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGatattggt                                                                                                                   <  1:150484/250‑5 (MQ=255)                                                                                                                                    ccTGGGGCCACCTCCGCGACGGCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGCCGGCGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGGTGGGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTg                                                                                                          <  2:109635/251‑1 (MQ=255)                                                                                                                                     cTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAcgcttg                                                                                                         >  2:295715/1‑246 (MQ=255)                                                                                                                                     cTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTgg                                                                                                         >  1:172214/1‑251 (MQ=255)                                                                                                                                      tGGGGCCACAGACGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGGTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGCCGACATCCTGATCCACACTGCTGCCCACGTTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACTGCACGAAAGTAAAACAGGTCTCCC‑TGGCACACCGTGAAACCCTCGCGGatattggttagactcatgc                                                                                                        >  2:259204/1‑236 (MQ=255)                                                                                                                                          gCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTa                                                                                                                                                                                      <  1:61594/170‑1 (MQ=255)                                                                                                                                          gCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTa                                                                                                                                                                                      >  2:61594/1‑170 (MQ=255)                                                                                                                                                gTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACCTCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGCCCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACATTCGCGGATATTGGACAGCCTCTGGCAAAATAtttt                                                                                              >  1:454689/1‑251 (MQ=255)                                                                                                                                                 gcgcgACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTg                                                                                             <  1:104009/250‑1 (MQ=255)                                                                                                                                                                    gctgcgGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAat                                                                           <  2:111667/247‑1 (MQ=255)                                                                                                                                                                              tGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGGTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGat                                                                                  <  1:446364/233‑1 (MQ=255)                                                                                                                                                                              tGGAGCTGTTCGACAGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCCCCGTGAAACCTTCGCGGATATTGGTCAGCCTCTGGCAAAATATTTTGGTACTTCTGat                                                                                  >  2:446364/1‑233 (MQ=255)                                                                                                                                                                              gggagcTGTTCGACCGCCGTCTGCCGGTGCTGATGTGTCTGCTGCGCGATGACGTCATCCTGATCCTCACTGCTGGCCAAGGGTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCTCCGTGAAACCTTCGCGGATATTGGGCAGACTCGGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGc                                                                 <  2:296977/245‑1 (MQ=255)                                                                                                                                                                                gagcTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCAt                                                              <  1:295715/248‑1 (MQ=255)                                                                                                                                                                                  gcTGTTCGACCGCCGTCTTCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAcc                                                                                                                                                            <  1:147678/155‑1 (MQ=255)                                                                                                                                                                                  gcTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGCCCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAcc                                                                                                                                                            >  2:147678/2‑156 (MQ=255)                                                                                                                                                                                   cTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCCAAATATTTTGGTACTTCTGATATGGAATATGACAAAGCCATGtt                                                           >  2:155883/1‑251 (MQ=255)                                                                                                                                                                                                       ggAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttg                                       <  2:258980‑M1/251‑16 (MQ=255)                                                                                                                                                                                                              aTGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgcccc                                >  2:469188‑M1/1‑229 (MQ=255)                                                                                                                                                                                                                gTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCGCGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgctgcttgatgccctc                               <  1:36075‑M1/250‑24 (MQ=255)                                                                                                                                                                                                                                gacgacATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGACAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctccc               <  1:155883‑M1/250‑40 (MQ=255)                                                                                                                                                                                                                                    acagcctgatCCGCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagg          <  2:172214‑M1/247‑45 (MQ=255)                                                                                                                                                                                                                                             tCCTCACTGCTGACCCCGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATGCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGTATATTGGTCAGACTCTGGCAAAATATTTTGGTACTGCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacatggagagggt  <  2:454689‑M1/250‑53 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAA  >  NC_000913/4620394‑4620879

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: ‑