Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,519	C→T	T48T (ACG→ACA)	dcuS ←	sensor histidine kinase DcuS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,519	0	C	T	95.3%	120.2 / ‑5.8	43	T48T (ACG→ACA)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (20/21);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACT  >  NC_000913/4351272‑4351758                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 cATTCAGCGCTTTAAGGATGTCATCACTTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACGTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCTCGATGGCCGGGATGCCACTCTCCTGCGTTTTTTTCTGCAAGCCCTGACGTATTTCCGGCGAGTCGTCGAGGGTACGCTCCACTGCCAGTGCCTTGTTGGCTCGCCCAGCTCGTGTc                                                                                                                                                                                                                                               <  2:36144/251‑1 (MQ=255)     cAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCatat                                                                                                                                                                                                                                           >  1:402785/1‑251 (MQ=255)           tttAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGAATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGGTGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTg                                                                                                                                                                                                                                      >  1:20149/1‑250 (MQ=255)            ttaaggatgtaccccccTTTAAATGGCTGACCAATACGCTGGGCTTAAGGATGCGAGTAGCGAAGTCTTTGCCACTCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCCGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTg                                                                                                                                                                                                                                      <  2:137994/234‑1 (MQ=255)             tAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGAt                                                                                                                                                                                                                                    <  2:223122/250‑1 (MQ=255)               aGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCCTCCCGTGTCATATCTATGATTTg                                                                                                                                                                                                                                 >  2:218830/1‑251 (MQ=255)                ggATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTTCGGCGAGTCGGCGAGGGGACGCGCCACTGCCAGCACCTTGTTGGCTAGCCCAGCTCGTGTCATATCACtagtttgc                                                                                                                                                                                                                                >  2:326449/1‑244 (MQ=255)                     tcatcaCCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCgaga                                                                                                                                                                                                                            >  1:347809/1‑250 (MQ=255)                       atcaCCTTTAAATGGCTGCCCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTgagatcggg                                                                                                                                                                                                                         >  2:162586/1‑243 (MQ=255)                       atcaCCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACGTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTg                                                                                                                                                                                                                                 <  1:162586/243‑1 (MQ=255)                            ttttAAATGGCTGACAAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGTTAACGACACTAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTTACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCAATGCCAGTGCCTTGTTAGCTAGCCCATCTAGTGTCATATCACTGCTTTGCGGGAAGTAAATc                                                                                                                                                                                                                    >  2:213441/2‑251 (MQ=255)                             tttAAACGGCTGACAAATCCGCAAGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATCTAGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                                    >  2:360660/1‑250 (MQ=255)                                ataTGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCTAAGACTTTGCATAACGGTAACGGCAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGTATTTCCGGCGAGTCGGCGAGGTTACGCGCCACTGCCAGTGCCTTGTTGGCCAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGAt                                                                                                                                                                                                                <  2:319494/249‑1 (MQ=255)                                aaaTGGCGGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGTTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGAt                                                                                                                                                                                                                <  2:319479/251‑1 (MQ=255)                                        gACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCAGTAACGACAATAAACAGCCGATCGTTGCGTTTGCGTACGGCTTCCGAGATGCACTGGGTGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGTGTCGGCGAGGGCACGCGCCACTGCCAGTGCCTCGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGCGAAGTAAATCAGATGCaccacc                                                                                                                                                                                                        >  1:195766/1‑251 (MQ=255)                                            aaTACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAAt                                                                                                                                                                                                     >  1:97103/1‑250 (MQ=255)                                                                    gTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGa                                                                                                                                                                            <  1:325050/251‑1 (MQ=255)                                                                            gACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAg                                                                                                                                                                    <  2:178808/251‑1 (MQ=255)                                                                                  gcttaccGGTAACGTCAATAAACATCAGATCGTTGCCTTTGCGTACGGCTTCCTCGATTGCCTGGATGCCACTCTCCTGCGGTTTTTTCTTCAAGCCCTGACGGAGTTCCGGCGAGTCGGCGAGGGTACGCGCCCCTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGGAGTAAATCAGATGCCCCACCAATAGCACCGAGAACAGTACCGCACTGAGTATTAAGATCACt                                                                                                                                                              <  2:177964/245‑1 (MQ=255)                                                                                      atCGGTAACGACAATAACCAGCAGACCGTTGCGTTTTCTTACGGTTGCCTCGATGGCCGGGATGCCACTCTCCGGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGATAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgtg                                                                                                                                                           <  2:222818/250‑1 (MQ=255)                                                                                        cGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                         >  1:320506/1‑250 (MQ=255)                                                                                              cGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAAt                                                                                                                                                 >  1:96358/1‑251 (MQ=255)                                                                                                       aCAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTAATCgg                                                                                                                                         >  2:280467/1‑250 (MQ=255)                                                                                                                cGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTCCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGt                                                                                                                               <  2:319044/251‑1 (MQ=255)                                                                                                                 gTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGt                                                                                                                               >  1:117428/1‑250 (MQ=255)                                                                                                                       tttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCACTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTATTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACAGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGc                                                                                                                        >  1:29005/1‑251 (MQ=255)                                                                                                                            gatcGGCTTGCACGATGGCGTGGAGTCCACTCGGCTGCGGTTTGTTCTGCAAGCCCTGACGGTTTTCCGGCGAGTCGGCGAGGGTACGAGCCACTGCCAGTGCCTTGTTGGCTAGATCATCTCGTGTCAGATCCCTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCAACGAGAACAGTACCGCACGGATCATTAAGATCACTGTG‑TAGATCAATTGCAGCGGACGTGTGTGTAACATGCGGTa                                                                                                                    <  2:82160/247‑1 (MQ=255)                                                                                                                              aCGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAgg                                                                                                                  <  2:117428/250‑1 (MQ=255)                                                                                                                                    tCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatg                                                                                                            <  1:379970/250‑1 (MQ=255)                                                                                                                                         gATGGCCTGGATGCCACTCTCCTGCGGTTTTTTATGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGCGAGTACGCGCCACAGCCAGTAACCTGTTGGCTAGCCAATCTAGCGTCACATCACTGATTTGCGAGAAGTAAATCAGATGCCCCGCCAATCGCACCGAGAAAAGTCCCGCACTGCTCATTCAGATAACTGTGGAA‑CTCAATTTCCTCGGCAGTTTGTGTAACATGCGGTAGGGCAATGAATGtc                                                                                                      >  1:394018/1‑251 (MQ=255)                                                                                                                                                 ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAAGTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGc                                                                                              >  1:307923/1‑251 (MQ=255)                                                                                                                                                  gATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCt                                                                                             <  1:60954/251‑1 (MQ=255)                                                                                                                                                     tccACTCTCCTGCCGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGACTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCACGCGTCCTATCAATGATTAGCGAGACGAAAATCAGATGCACCCCCACTAGCACCGAGAACAATACCGCACTGATCATTAAGATCACTGTAGTA‑CTCAATTTACTCGGCCGTATGCGAAACATGCGCAAGGGCAATGAATGTCTCATCAgtttac                                                                                          >  2:215266/2‑246 (MQ=255)                                                                                                                                                      ccACTCTCCTGCGGTTTTTTCTGCAAGCCATGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCt                                                                                         >  2:241064/1‑251 (MQ=255)                                                                                                                                                         ctctCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgt                                                                                      >  2:286341/1‑251 (MQ=255)                                                                                                                                                               tGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGGA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGAc                                                                                 >  1:77136/1‑250 (MQ=255)                                                                                                                                                                  ggTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATAGCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCATCGTCCTTGTGTGACAAAt                                                                             <  2:77136/251‑1 (MQ=255)                                                                                                                                                                      ttttCTGCAATCCCTTACGGATTTCCTGCGAGCCGGCGAGTGTACGCGCCCCTGCCAGTGTCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCTCCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCAGCAGCTTCCTTGTGTGACAAATTTc                                                                          <  2:12515/250‑1 (MQ=255)                                                                                                                                                                                            ttCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGc                                                   >  1:402617/1‑251 (MQ=255)                                                                                                                                                                                                                      ccACTGCCAGTGTCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGAATTATCTCTGATGAGGCGTGTAATTCAAAGGGAGTAAGAATGa                         <  1:257953/251‑1 (MQ=255)                                                                                                                                                                                                                                 gCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCtatata              <  1:286341/251‑1 (MQ=255)                                                                                                                                                                                                                                      gttgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAgggg          <  1:218830/250‑1 (MQ=255)                                                                                                                                                                                                                                             aGCCCATCTCGTGTCATATCACTGATGTGCGAGAAGTAAATCAGATGCACCACCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGGGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGAct  <  1:251681/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 CATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTA‑CTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACT  >  NC_000913/4351272‑4351758

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