Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,142,505	G→C	7.1%	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.1%	336.8 / 14.9	98	R148R (CGC→CGG)	udk	uridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (42/49);  new base C (2/5);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

CCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGC  >  NC_000913/2142364‑2142646                                                                                                                                              |                                                                                                                                              cccGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCa                                                                                                                                        >  4:160585/1‑149 (MQ=255) cccGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCa                                                                                                                                        <  2:756022/149‑1 (MQ=255)   cgccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCCTGATGCGCCGCATCaag                                                                                                                                      <  2:600823/149‑3 (MQ=255)   cgccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCCTGATGCGCCGCATCaag                                                                                                                                      <  4:109073/149‑3 (MQ=255)    gccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGc                                                                                                                                     <  2:76114/149‑1 (MQ=255)    gccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCCTGATGCGCCGCATCAAGc                                                                                                                                     <  3:319451/149‑1 (MQ=255)     ccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCa                                                                                                                                    <  2:545557/149‑1 (MQ=255)         gcgGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGata                                                                                                                                <  2:43774/149‑1 (MQ=255)            gCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATcc                                                                                                                             >  4:286581/1‑149 (MQ=255)                gATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCg                                                                                                                         >  2:329526/1‑149 (MQ=255)                 aTAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                        <  2:449032/149‑1 (MQ=255)                 aTAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                        >  2:573777/1‑149 (MQ=255)                  tAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCggg                                                                                                                       <  1:343530/149‑1 (MQ=255)                       tCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATc                                                                                                                  >  1:238782/1‑149 (MQ=255)                       tCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATc                                                                                                                  >  4:113317/1‑149 (MQ=255)                         cGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCaa                                                                                                                <  3:695342/149‑1 (MQ=255)                                 gTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATg                                                                                                        >  4:740316/1‑149 (MQ=255)                                 gTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATg                                                                                                        <  1:128810/149‑1 (MQ=255)                                  tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                                                                       >  2:233942/1‑149 (MQ=255)                                   ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                      <  4:149751/149‑1 (MQ=255)                                   ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                      <  2:356921/149‑1 (MQ=255)                                   ttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                      <  2:726253/149‑1 (MQ=255)                                    tAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgag                                                                                                     <  1:245914/149‑1 (MQ=255)                                      gAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   <  3:279895/149‑1 (MQ=255)                                      gAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   <  1:329526/149‑1 (MQ=255)                                            tCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                             >  4:587669/1‑149 (MQ=255)                                            tCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                             >  2:701483/1‑149 (MQ=255)                                              aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagtt                                                                                           >  4:175155/1‑149 (MQ=255)                                               aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTc                                                                                          >  1:683532/1‑149 (MQ=255)                                               aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTc                                                                                          <  4:676938/149‑1 (MQ=255)                                                tGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCg                                                                                         >  1:493487/1‑149 (MQ=255)                                                   aTTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCa                                                                                      <  1:276259/149‑1 (MQ=255)                                                      gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   <  1:686941/133‑1 (MQ=255)                                                      gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   >  2:686941/1‑133 (MQ=255)                                                      gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGc                                                                                   <  2:261452/149‑1 (MQ=255)                                                       cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCa                                                                                  >  1:636244/1‑149 (MQ=255)                                                        aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                                                              <  4:236204/136‑1 (MQ=255)                                                        aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                                                              >  3:236204/1‑136 (MQ=255)                                                        aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   >  3:514031/1‑131 (MQ=255)                                                        aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                                                                   <  4:514031/131‑1 (MQ=255)                                                        aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCaa                                                                                 >  1:201881/1‑149 (MQ=255)                                                          gAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                               <  3:160585/149‑1 (MQ=255)                                                              aTCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACgcgc                                                                           >  3:724227/1‑149 (MQ=255)                                                               tCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCa                                                                          <  3:530720/149‑1 (MQ=255)                                                                   gcgcACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCg                                                                      <  2:669417/149‑1 (MQ=255)                                                                     gcACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTc                                                                    >  3:471202/1‑149 (MQ=255)                                                                      cACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCa                                                                   <  1:233942/149‑1 (MQ=255)                                                                           ttttttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGcaac                                                              >  1:126221/1‑149 (MQ=255)                                                                             ttttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGcaacaa                                                            >  2:712906/1‑149 (MQ=255)                                                                              tttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACaaa                                                           >  2:446411/1‑149 (MQ=255)                                                                               ttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACaaaa                                                          <  4:683852/149‑1 (MQ=255)                                                                                  atatTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGc                                                       <  2:176846/149‑1 (MQ=255)                                                                                     ttGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGATATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCtt                                                    <  4:13118/149‑1 (MQ=255)                                                                                       gcgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCg                                                  >  4:136590/1‑149 (MQ=255)                                                                                       gcgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCg                                                  <  3:630625/149‑1 (MQ=255)                                                                                        cgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCga                                                 <  4:724227/149‑1 (MQ=255)                                                                                           cATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGaa                                              >  3:648707/1‑149 (MQ=255)                                                                                            aTCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCTTCCGTCAGCAACAAAATGCCTTCGAGAat                                             <  3:451951/149‑1 (MQ=255)                                                                                               aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgat                                          <  4:631040/149‑1 (MQ=255)                                                                                                  gAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGAc                                       <  1:174129/149‑1 (MQ=255)                                                                                                   aaTCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGAcc                                      <  1:1247/149‑1 (MQ=255)                                                                                                        aTTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACcttctt                                 <  3:175155/149‑1 (MQ=255)                                                                                                            aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             <  3:587634/149‑1 (MQ=255)                                                                                                            aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             >  4:168528/1‑149 (MQ=255)                                                                                                            aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             >  4:587634/1‑149 (MQ=255)                                                                                                            aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             >  1:442849/1‑149 (MQ=255)                                                                                                            aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             <  2:493487/149‑1 (MQ=255)                                                                                                            aaCGCCCACGCTCGCTAACGTCACGCTTGATCCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAATACGGAGAAGTTAAGTTCGTCACGCAAACGCGCATTCGTGAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                             <  3:468615/149‑1 (MQ=255)                                                                                                              cGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTc                           <  3:168528/149‑1 (MQ=255)                                                                                                                 ccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                               <  4:356611/94‑1 (MQ=255)                                                                                                                 ccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                               >  3:356611/1‑94 (MQ=255)                                                                                                                 ccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAAc                        >  4:518378/1‑149 (MQ=255)                                                                                                                   aCGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCg                      >  4:374515/1‑149 (MQ=255)                                                                                                                      cTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCa                   <  2:266318/149‑1 (MQ=255)                                                                                                                       tCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCAc                  >  3:40616/1‑149 (MQ=255)                                                                                                                       tCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCAc                  >  2:219363/1‑149 (MQ=255)                                                                                                                        cGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGGTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCAcc                 <  3:489778/149‑1 (MQ=255)                                                                                                                         gTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCg                >  1:32709/1‑149 (MQ=255)                                                                                                                          ttAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGt               >  4:524815/1‑149 (MQ=255)                                                                                                                            aaCGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGttt             >  4:297966/1‑149 (MQ=255)                                                                                                                             aCGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCGACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttc            <  2:636244/149‑1 (MQ=255)                                                                                                                               gTCTCGCTTGATGCGGCGCATCAGGAAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGt                                                                     <  2:434670/90‑1 (MQ=255)                                                                                                                               gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGt                                                                     >  1:434670/1‑90 (MQ=255)                                                                                                                               gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt          >  4:131837/1‑149 (MQ=255)                                                                                                                               gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt          <  4:129971/149‑1 (MQ=255)                                                                                                                               gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt          <  3:297966/149‑1 (MQ=255)                                                                                                                               gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt          <  3:286581/149‑1 (MQ=255)                                                                                                                                ttctgcctGATCCGCCGCATCAGGCAGATACCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACACGCATCCGTCAGCAACAAATTGCCTTCAAGAATGATGACCTTCTTCGGCTCCACCGTCACCGtttcttt         <  4:633855/142‑1 (MQ=255)                                                                                                                                tCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttcttt         >  3:683190/1‑149 (MQ=255)                                                                                                                                 cACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctttc        <  3:102834/149‑1 (MQ=255)                                                                                                                                  ctgcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCATCGAAAATGGAGATGTTAACTTCGTCACGCAAAAGAGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACGTTCTTCGGCTCAACCGACACCGCTTCTTTCa       >  3:633855/6‑149 (MQ=255)                                                                                                                                  aCGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCAc                  >  4:533816/1‑138 (MQ=255)                                                                                                                                  aCGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCAc                  <  3:533816/138‑1 (MQ=255)                                                                                                                                  aCGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa       >  1:324231/1‑149 (MQ=255)                                                                                                                                    gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa     >  4:241448/4‑149 (MQ=255)                                                                                                                                    gCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa     <  1:75158/149‑1 (MQ=255)                                                                                                                                      ttGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcg   >  3:583823/1‑149 (MQ=255)                                                                                                                                       tGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcgc  <  4:221243/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              CCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGC  >  NC_000913/2142364‑2142646

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