Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,620,569:1	+A	9.0%	coding (967/1224 nt)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,569	1	.	A	9.0%	240.8 / 21.9	67	coding (967/1224 nt)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base . (27/34);  new base A (3/3);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

AAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCG  >  NC_000913/4620422‑4620715                                                                                                                                                     |                                                                                                                                                   aaaGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTgg                                                                                                                                                     >  2:65051/1‑148 (MQ=255)       aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑aaactgt                                                                                                                                             <  2:224402/149‑6 (MQ=255)       aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGtt                                                                                                                                             <  4:312885/149‑1 (MQ=255)            ccTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGAcc                                                                                                                                        <  1:65051/149‑1 (MQ=255)             cTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AAc                                                                                                                                                 >  4:188165/1‑139 (MQ=255)             cTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AAc                                                                                                                                                 <  3:188165/139‑1 (MQ=255)                 aCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGt                                                                                                                                   >  4:148627/1‑149 (MQ=255)                  cgcgCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTc                                                                                                                                  <  8:145413/149‑1 (MQ=255)                  cgcgCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTc                                                                                                                                  <  6:247199/149‑1 (MQ=255)                   gcgcTGTTTGACTCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCt                                                                                                                                 <  1:306203/149‑1 (MQ=255)                     gcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGc                                                                                                                               >  3:344142/1‑149 (MQ=255)                        gTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATATCACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGACGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCg                                                                                                                             >  1:293006/1‑148 (MQ=255)                              cGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCTTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGa                                                                                                                      <  8:135055/149‑1 (MQ=255)                               gCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGAt                                                                                                                     >  3:349524/1‑149 (MQ=255)                                 caccaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGt                                                                                                                   >  4:296511/1‑149 (MQ=255)                                 caccaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGt                                                                                                                   >  7:283647/1‑149 (MQ=255)                                  accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGtc                                                                                                                  >  1:393638/1‑149 (MQ=255)                                        aaaGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                                                                            >  4:356367/1‑149 (MQ=255)                                         aaGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGATGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgctg                                                                                                           >  6:22484/1‑149 (MQ=255)                                          agagATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCtgctgc                                                                                                          >  5:281245/1‑149 (MQ=255)                                            ataTGAAAGAAGCGGGTGATAAAACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcgc                                                                                                        <  2:293006/147‑1 (MQ=255)                                                   agaagCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATgac                                                                                                  <  5:296196/149‑1 (MQ=255)                                                         gggTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATcc                                                                                           <  3:119981/149‑1 (MQ=255)                                                           gTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCt                                                                                          >  1:163982/1‑148 (MQ=255)                                                           gTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCt                                                                                          >  3:356542/1‑148 (MQ=255)                                                                aaCACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATgacgac                                                                                               <  8:279401/138‑1 (MQ=255)                                                                aaCACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATgacgac                                                                                               >  7:279401/1‑138 (MQ=255)                                                                  cacCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCa                                                                                  >  1:342253/1‑149 (MQ=255)                                                                    ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCAcc                                                                                <  4:29327/149‑1 (MQ=255)                                                                        cGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTg                                                                            <  3:296511/149‑1 (MQ=255)                                                                         gTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCT‑‑‑‑‑CTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCa                                                                        <  3:148790/148‑1 (MQ=255)                                                                              cACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCa                                                                        <  4:344142/147‑1 (MQ=255)                                                                                ccAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGt                                                                    >  4:329147/1‑149 (MQ=255)                                                                                 cAACTTCGTTGACTTCGAGTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGtt                                                                   <  3:78256/149‑1 (MQ=255)                                                                                  aaCTGCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTg                                                                  <  7:191042/149‑1 (MQ=255)                                                                                  aaCGTCGTTGACTTCGACTCTTCCTGGGGCGACCGGCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGTTGACCACGGTTg                                                                  <  5:237340/149‑1 (MQ=255)                                                                                    cTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGc                                                                 >  3:261161/1‑148 (MQ=255)                                                                                            aCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAcc                                                        <  7:2744/149‑1 (MQ=255)                                                                                            aCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAcc                                                        <  3:329147/149‑1 (MQ=255)                                                                                            aCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGAc                                                         <  7:129282/148‑1 (MQ=255)                                                                                             cTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCt                                                       <  3:148627/149‑1 (MQ=255)                                                                                                cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTg                                                      <  2:393638/147‑1 (MQ=255)                                                                                                cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGa                                                    >  8:280453/1‑149 (MQ=255)                                                                                                cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGa                                                    >  5:271822/1‑149 (MQ=255)                                                                                                cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCAACGCTGACCACGGTTGAGATCCGACCTg                                                      <  8:283647/147‑1 (MQ=255)                                                                                                 gACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCTCCGCTGACCACGGTTGCGATCCGACCTGGa                                                    <  4:124551/148‑1 (MQ=255)                                                                                                      ttCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCgg                                                <  6:328568/147‑1 (MQ=255)                                                                                                      ttCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTa                                              >  2:390695/1‑149 (MQ=255)                                                                                                      ttCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTa                                              >  4:301992/1‑149 (MQ=255)                                                                                                            gggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCTGAGCTGATGTCTCTGCTGCGCTATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGAc                                         <  6:99603/148‑1 (MQ=255)                                                                                                            gggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACt                                            >  3:306191/1‑145 (MQ=255)                                                                                                            gggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACt                                            <  4:306191/145‑1 (MQ=255)                                                                                                            gggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGAc                                         <  4:140136/148‑1 (MQ=255)                                                                                                             ggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGAc                                         >  3:53281/1‑148 (MQ=255)                                                                                                              gCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACca                                       >  3:170130/1‑149 (MQ=255)                                                                                                                cACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGCCATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCacac                                    >  3:206992/1‑149 (MQ=255)                                                                                                                cACCGTCGCGACGTCGCAGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGAGATCCGACCTGGACCGGTACTGACCacac                                    >  5:135921/1‑149 (MQ=255)                                                                                                                   cGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATgacgac                                                                                               >  6:282241/1‑88 (MQ=255)                                                                                                                   cGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATgacgac                                                                                               <  5:282241/88‑1 (MQ=255)                                                                                                                    gTCGCGATGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGt                                <  5:22484/149‑1 (MQ=255)                                                                                                                    gTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGt                                <  6:271822/149‑1 (MQ=255)                                                                                                                    gTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTGCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGt                                <  2:260363/149‑1 (MQ=255)                                                                                                                             tCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACAtt                       >  4:245124/1‑149 (MQ=255)                                                                                                                               gCCGGTTATGCTGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTc                      <  3:325591/148‑1 (MQ=255)                                                                                                                                caggTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCg                    <  6:135921/147‑1 (MQ=255)                                                                                                                                  ggTTATGCCGCGGGTCTGGAAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCgg                   <  4:170130/149‑1 (MQ=255)                                                                                                                                           gcgGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtata         <  1:390695/149‑1 (MQ=255)                                                                                                                                               gTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGc     >  2:26429/1‑149 (MQ=255)                                                                                                                                                  tGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGccc   >  5:48068/1‑148 (MQ=255)                                                                                                                                                   gg‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCg  >  3:227360/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGG‑AACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCG  >  NC_000913/4620422‑4620715

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