Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,348,366	C→A	T328K (ACG→AAG)	nrdB →	ribonucleoside‑diphosphate reductase 1 subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,348,366	0	C	A	97.7%	127.9 / ‑4.1	43	T328K (ACG→AAG)	nrdB	ribonucleoside‑diphosphate reductase 1 subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (28/14);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

CCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGC  >  NC_000913/2348120‑2348608                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    cccAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAgcgc                                                                                                                                                                                                                                                 >  2:265629/1‑251 (MQ=255)               aaTCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATccc                                                                                                                                                                                                                                   <  1:460418/251‑1 (MQ=255)                     ctgcGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGAt                                                                                                                                                                                                                             >  1:178751/1‑251 (MQ=255)                           aGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAacac                                                                                                                                                                                                                       >  1:228604/1‑251 (MQ=255)                                cgcgGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGc                                                                                                                                                                                                                  >  1:364816/1‑251 (MQ=255)                                      cGATCCTGAGTTGGCGGAAATTGCCGAAGAGTGTAAGCAGGATTGCTTTGACCTGTTTGTTCAGGCAGCTCAACCGGAGAAAGCCTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATCCAGACATTCTCTGCCAGTACGTTGAATCCCTCCCCATTATCCTTATGCAGTCAGTCGGTTTGGATCTGCCGTTCCAGAAGCCCTCCAACCCGGTCCCGTGGCTCACCACTTGCCTggggt                                                                                                                                                                                                            >  2:53510/1‑248 (MQ=255)                                           cTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGTGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGCTCTGCAGTTCCAGAAGCGCTCCAACACGATCCCGTGGATCAACACTTGGCTGATGTCTGAt                                                                                                                                                                                                       >  1:268231/1‑251 (MQ=255)                                                    cGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAg                                                                                                                                                                                              <  1:265629/251‑1 (MQ=255)                                                      gAAATTGCCGAAGCGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGAGCTCCAACCCGATCACGTGGATCACCCCTTGGCTGGTGTCTGCTAACGTGCAGGt                                                                                                                                                                                            >  2:391007/1‑251 (MQ=255)                                                       aagtagCCGATGAATGTAAGCAGAGGTGCTATGACTTTTTTGTACAGGCAACTCAACAGTAGAAAGACTGTGCGGATTATCTGGTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACACCACCAATCTCCGTAGGCAGGCAGTCGGTTTGGATCTTCCGTTCCAGAAGCGCTCCAACCCGTTCCCTTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAgggt                                                                                                                                                                                           <  2:402602/246‑3 (MQ=255)                                                        aaTTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGTATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTg                                                                                                                                                                                          <  1:415559/251‑1 (MQ=255)                                                            gCCGAAGAGTGTAAGCAGGAGTGCTTTGCCCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTcc                                                                                                                                                                                      <  2:364816/251‑1 (MQ=255)                                                             ccgcagagTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCTTGTATCAACACTTGGCTGGTGTCTGATAACGTCCCGGTTGCTccc                                                                                                                                                                                     >  1:192399/5‑250 (MQ=255)                                                               gAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAAGCAGCTCAACAGGAGAACGACTGGGCGGATTATCTGTTCCGCGACGGTTAGATGGTTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGTTCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACAATTGGCTGGTGTCTGATAACGTGCATGTTGCTCCGcc                                                                                                                                                                                   >  2:46790/1‑250 (MQ=255)                                                                            aGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGGGTCTGATAACGAGCAGGTTGCTCCGCAGGAAGTGGAAGTc                                                                                                                                                                      >  1:468003/1‑251 (MQ=255)                                                                             ggAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACt                                                                                                                                                                                                                      >  1:215832/1‑202 (MQ=255)                                                                             ggAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACt                                                                                                                                                                                                                      <  2:215832/202‑1 (MQ=255)                                                                                      tGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATc                                                                                                                                                            >  1:29595/1‑251 (MQ=255)                                                                                            gtttgttCAGGCAGCGCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGt                                                                                                                                                        <  1:95620/249‑1 (MQ=255)                                                                                                gttCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGc                                                                                                                                                   >  1:329076/1‑250 (MQ=255)                                                                                                           tCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATATGCCGTTCCAGAAGCGCACAAACCCGATCCCGTGGATCAACACTTGGCTGGGGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCCGTTCTTATCTGGTCGGGACGATTG‑ACTCgg                                                                                                                                      >  1:428882/1‑251 (MQ=255)                                                                                                               cAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGACTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCCGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGGTTG‑ACTCGGAAgg                                                                                                                                  >  2:461166/1‑250 (MQ=255)                                                                                                                      aaGACTGGGCGGCTTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACAcc                                                                                                                           >  1:203840/1‑251 (MQ=255)                                                                                                                                         gTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGTTAACGTGCAGGTTGCTCCGCATGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGCTTTGAGTAACTTcc                                                                                                        >  2:91961/1‑251 (MQ=255)                                                                                                                                            ccGCGACGGTTGTATGATTGGTCTGAATAAAGACATTCCCTGCCAGTACGTGGAATACATCACCAATTGCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCTGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGc                                                                                                     <  2:415399/251‑1 (MQ=255)                                                                                                                                                       cGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGccc                                                                                          >  1:450085/1‑251 (MQ=255)                                                                                                                                                                gTCTGAATAAAGACCTTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGACCCGCGTTAccc                                                                                 >  1:352874/1‑251 (MQ=255)                                                                                                                                                                      aTAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTgcgc                                                                            >  1:390027/1‑250 (MQ=255)                                                                                                                                                                               ttCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGCTTTGAGTAACTTCCAGCTCTGATGCCCCGCGTTACCCTGCGCATCACTGGc                                                                   >  2:326226/1‑250 (MQ=255)                                                                                                                                                                               ttCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTa                                                                                                              >  2:418489/1‑207 (MQ=255)                                                                                                                                                                               ttCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTa                                                                                                              <  1:418489/207‑1 (MQ=255)                                                                                                                                                                                tctcTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGcac                                                                 >  1:231123/1‑251 (MQ=255)                                                                                                                                                                                      ccaggacgggtaaTACATCCCCAAGATCCGTATGCAGGCAGGCGGTTTGGATCTGCCGGTCCAGAAGCGCGCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGAGGGCCCGCGTTACCCTGCGCATCACTGGCACACAAct                                                           <  2:468003/240‑1 (MQ=255)                                                                                                                                                                                       cAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCtt                                                                                                                                                               <  1:185066/151‑1 (MQ=255)                                                                                                                                                                                       cAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCtt                                                                                                                                                               >  2:185066/1‑151 (MQ=255)                                                                                                                                                                                              ttGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTTCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCACCACTTGGCTGGCGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGCGGAAGTCTGTTCTTATCTGGTCGGGCAGTTTGCAC‑CGGGAGTGGACCCCGACGCTTTGAGTAACTTCCAGCTCTGATGGCACGGGTTACCCTGCGCATCACTGGCCCGCAACTGCTGTGcc                                                   >  2:133490/1‑251 (MQ=255)                                                                                                                                                                                                       tCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTAGCACACAACTCCTTTGCCAGGATGAac                                          >  2:104540/1‑251 (MQ=255)                                                                                                                                                                                                        cACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAac                                          <  2:272216/250‑1 (MQ=255)                                                                                                                                                                                                        cACCAATAGCCGTATGCAGGCAGTCGGTGTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACGTCCAGCTCTGATGGCCCGCGTTACCCTGCGCAACACTGGCACACAACTGCTGTGCCAGGATGAac                                          <  1:469378/250‑1 (MQ=255)                                                                                                                                                                                                               atCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACAcccttcc                                  >  2:164896/1‑251 (MQ=255)                                                                                                                                                                                                                                          cTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAAt       >  2:173149/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGGCAGTTCTTATCTGGTCGGGCAGATTG‑ACGCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCcaca         <  1:461166/247‑1 (MQ=255)                                                                                                                                                                                                                                               gTTCAAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATACCGTGCAGGTTGCTCCGCAAGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGCAGTGGACACCGACGATTTGAGTAACTTACAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAAAACCCTTCCCTTCTGGCGGCGCTGGagaccaaaaatgtggc  >  1:32959/1‑236 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTG‑ACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGC  >  NC_000913/2348120‑2348608

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