Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,142,505	G→C	7.5%	R148R (CGC→CGG)	udk ←	uridine/cytidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,142,505	0	G	C	7.5%	316.0 / 15.0	93	R148R (CGC→CGG)	udk	uridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (41/45);  new base C (3/4);  total (44/49)
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.54e-01

GGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTA  >  NC_000913/2142357‑2142649                                                                                                                                                     |                                                                                                                                                 ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcgg                                                                                                                                                  <  3:271222/149‑1 (MQ=255)  gTTTTTCCCGCCGCGCGGCACGATAATGACCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                 >  3:65290/1‑149 (MQ=255)        cccGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATc                                                                                                                                            <  1:219000/148‑1 (MQ=255)         ccgccgCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCa                                                                                                                                           >  2:385234/1‑148 (MQ=255)          cgccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaa                                                                                                                                          <  1:748316/148‑2 (MQ=255)           gccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGc                                                                                                                                        <  4:876808/149‑1 (MQ=255)            ccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                                                                        <  2:710659/148‑1 (MQ=255)            ccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                                                                        <  2:630395/148‑1 (MQ=255)            ccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                                                                        >  1:778877/1‑148 (MQ=255)            ccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCa                                                                                                                                       <  4:460249/149‑1 (MQ=255)             cgcgcgGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAg                                                                                                                                      <  3:841000/149‑1 (MQ=255)                  ggCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATc                                                                                                                                 >  3:87440/1‑149 (MQ=255)                         tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCggg                                                                                                                          >  4:935353/1‑149 (MQ=255)                         tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCggg                                                                                                                          <  4:403115/149‑1 (MQ=255)                         tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                           <  2:941738/148‑1 (MQ=255)                         tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                           <  2:141112/148‑1 (MQ=255)                         tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                           <  2:366951/148‑1 (MQ=255)                            tGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTa                                                                                                                       <  4:800933/149‑1 (MQ=255)                             gTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTAt                                                                                                                      >  3:487242/1‑149 (MQ=255)                                cGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCa                                                                                                                    <  2:161750/148‑1 (MQ=255)                                     atTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGaaa                                                                                                              <  3:326719/149‑1 (MQ=255)                                     atTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATAGTGCGCGATGACTGAATCCATGGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGaa                                                                                                               >  2:187198/1‑148 (MQ=255)                                       tGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGaaaa                                                                                                             >  1:575399/1‑148 (MQ=255)                                        gTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAAt                                                                                                            <  2:731270/148‑1 (MQ=255)                                         tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                                                                          <  4:166529/149‑1 (MQ=255)                                          ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                                                                          >  2:274785/1‑148 (MQ=255)                                          ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                                                                          >  1:409898/1‑148 (MQ=255)                                            aGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgag                                                                                                        <  1:45152/148‑1 (MQ=255)                                             gAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                      >  3:678206/1‑149 (MQ=255)                                               aGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaag                                                                                                     <  1:3923/148‑1 (MQ=255)                                                  cTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                                                                 <  3:126194/149‑1 (MQ=255)                                                   tCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCACATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                >  4:658275/1‑149 (MQ=255)                                                    cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                >  1:916372/1‑148 (MQ=255)                                                    cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                >  1:357844/1‑148 (MQ=255)                                                    cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaag                                                                                                                                         >  1:159108/1‑105 (MQ=255)                                                    cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaag                                                                                                                                         <  2:159108/107‑3 (MQ=255)                                                     aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagt                                                                                               >  2:402888/1‑148 (MQ=255)                                                      aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagtt                                                                                              >  2:397709/1‑148 (MQ=255)                                                           ttGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCAc                                                                                        <  3:393777/149‑1 (MQ=255)                                                               aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGat                                                                                                                                    >  3:257705/1‑101 (MQ=255)                                                               aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGat                                                                                                                                    <  4:257705/101‑1 (MQ=255)                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCaaa                                                                                   <  4:102821/149‑1 (MQ=255)                                                                 gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCaaa                                                                                   <  2:580166/148‑1 (MQ=255)                                                                 gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  <  4:612344/149‑1 (MQ=255)                                                                 gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  >  3:151375/1‑149 (MQ=255)                                                                 gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  >  3:738314/1‑149 (MQ=255)                                                                 gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCACATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  <  3:658275/149‑1 (MQ=255)                                                                  aaCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCg                                                                                                                            <  2:291440/106‑1 (MQ=255)                                                                  aaCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCg                                                                                                                            >  1:291440/1‑106 (MQ=255)                                                                  aaCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaa                                                                                                                                          <  1:495745/92‑2 (MQ=255)                                                                  aaCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaa                                                                                                                                          >  2:495745/1‑91 (MQ=255)                                                                       cGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCa                                                                             >  2:923274/1‑148 (MQ=255)                                                                        gggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCAt                                                                            <  2:926632/148‑1 (MQ=255)                                                                         ggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATc                                                                           <  1:385234/148‑1 (MQ=255)                                                                               cGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgag                                                                                                        >  1:806659/1‑113 (MQ=255)                                                                               cGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgag                                                                                                        <  2:806659/113‑1 (MQ=255)                                                                                            ttGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  >  3:612211/1‑122 (MQ=255)                                                                                            ttGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                  <  4:612211/122‑1 (MQ=255)                                                                                             tGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCtt                                                       <  2:387023/148‑1 (MQ=255)                                                                                               cgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCg                                                     <  2:357844/148‑1 (MQ=255)                                                                                                  cATCCCTGAATCCATTGAACGCCCACGCCGGTTAACGTCACGCTTGATGAGGCGCATCAGGCAGATATCCAGCGCGGTATCAACGAAAATGGAGAAGTCAAGTTCGTCACGCAAATGCGCCTCCGTCAGCAAAACAATGCCTTCgagac                                                 <  4:910370/149‑2 (MQ=255)                                                                                                  cATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGaa                                                 <  3:60069/149‑1 (MQ=255)                                                                                                   aTCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAat                                                <  3:889273/149‑1 (MQ=255)                                                                                                    tCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAat                                                >  2:594777/1‑148 (MQ=255)                                                                                                     cACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatg                                               >  1:7582/1‑148 (MQ=255)                                                                                                      aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatga                                              <  1:594777/148‑1 (MQ=255)                                                                                                      aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatga                                              <  2:623982/148‑1 (MQ=255)                                                                                                       cTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgatg                                            >  4:943325/1‑149 (MQ=255)                                                                                                       cTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgat                                             <  2:200979/148‑1 (MQ=255)                                                                                                         gAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGAc                                          >  4:230119/1‑149 (MQ=255)                                                                                                         gAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATACAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAAtgatga                                           <  1:39155/148‑1 (MQ=255)                                                                                                           aTCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                <  3:279315/93‑1 (MQ=255)                                                                                                           aTCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                >  4:279315/1‑93 (MQ=255)                                                                                                               aTTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTc                                                      <  1:906461/131‑1 (MQ=255)                                                                                                               aTTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTc                                                      >  2:906461/1‑131 (MQ=255)                                                                                                                ttGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACcttctt                                    >  1:362294/1‑148 (MQ=255)                                                                                                                  gAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCg                                  >  1:635601/1‑148 (MQ=255)                                                                                                                   aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCgg                                 >  2:375375/1‑148 (MQ=255)                                                                                                                        ccACGCTCGTTAACGTCACGCTTGATGCGGCGCTTCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAAc                           >  3:812309/1‑149 (MQ=255)                                                                                                                         cACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAAc                           >  2:75032/1‑148 (MQ=255)                                                                                                                          aCGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAAcc                          >  2:558257/1‑148 (MQ=255)                                                                                                                           cGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCg                         >  1:158172/1‑148 (MQ=255)                                                                                                                            gCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTc                       <  4:678206/149‑1 (MQ=255)                                                                                                                             cTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTc                       >  2:775442/1‑148 (MQ=255)                                                                                                                                   aaCGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtt                 <  1:274785/148‑1 (MQ=255)                                                                                                                                     cGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttct              <  3:398940/149‑1 (MQ=255)                                                                                                                                      gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt             >  3:756716/1‑149 (MQ=255)                                                                                                                                      gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt             <  4:72885/149‑1 (MQ=255)                                                                                                                                        cACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCAAGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttcttt            <  1:75032/148‑1 (MQ=255)                                                                                                                                         aCGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa          >  3:912053/1‑149 (MQ=255)                                                                                                                                          cGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa          >  1:883426/1‑148 (MQ=255)                                                                                                                                          cGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa          >  1:569566/1‑148 (MQ=255)                                                                                                                                           gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGAAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa        >  4:643669/4‑149 (MQ=255)                                                                                                                                            cTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa        >  1:173008/1‑148 (MQ=255)                                                                                                                                                 tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTa  <  4:677018/149‑1 (MQ=255)                                                                                                                                                  gcggcgCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTa  <  2:916372/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTA  >  NC_000913/2142357‑2142649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑