Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,652	C→T	C350C (TGC→TGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,652	0	C	T	96.2%	156.1 / NA	54	C350C (TGC→TGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (17/33);  minor base G (0/2);  total (19/35)
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

GTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCT  >  NC_000913/4620405‑4620888                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               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                                                                                                                                                                                               >  2:61594/1‑170 (MQ=255)                                                                                                                               gCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTa                                                                                                                                                                                               <  1:61594/170‑1 (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)                                                                                                                                                                                                                                          gcgtgccacggtTGGGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGTGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTTGTCAGACTCGGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagatggtgaaaacatg  <  2:83605‑M1/246‑62 (MQ=255)                                                                                                                                                                                                                                              aCCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCaaaa                                                                                                             >  2:100898/1‑140 (MQ=255)                                                                                                                                                                                                                                              aCCACGGTTGTGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCaaaa                                                                                                             <  1:100898/140‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGA‑TCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAC‑TGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCT  >  NC_000913/4620405‑4620888

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