Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	93.5%	118.6 / 2.8	46	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (21/22);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTG  >  NC_000913/3815636‑3816126                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTGGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCGCATCAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCTTCTGCTTCGCTACAATGTATTCGA‑T‑T‑ccccctc                                                                                                                                                                                                                                                   <  2:23123/251‑4 (MQ=255)      ccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGcc                                                                                                                                                                                                                                              >  1:200874/1‑251 (MQ=255)       cGGTATTAAACAGCCCGGCGTTGAATAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAATCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCa                                                                                                                                                                                                                                             <  2:408255/251‑1 (MQ=255)                           ttGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGc                                                                                                                                                                                                                         >  1:209562/1‑251 (MQ=255)                                         gCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                            >  1:203132/1‑250 (MQ=255)                                                       aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                             >  1:341906/1‑251 (MQ=255)                                                       aTTTCAGCGTAAACTCGCCAAACTTTAACACATGCTTGCTACGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATACCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                             >  1:65528/1‑251 (MQ=255)                                                         ttCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTCCTAAGCGCAAATTCAAGAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGATTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGAGTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                            <  1:270787/250‑1 (MQ=255)                                                         ttCAGCGTAAACTCGCCAAACTTTAACCTCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTTCAGACAAAAAAAAGTCTACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGTGCCAGGTCCTACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                            <  1:396464/250‑1 (MQ=255)                                                         ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAACCTGCGCTGATATGTTTTCATTCCTTCTCTCCTCATCTTTCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAGTCGA‑T‑T‑CCCCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                            <  2:336158/250‑1 (MQ=255)                                                                      cGCCAAACTTTAACACCTGCTTGCTAAGCTTAAATTAAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACGTTTCTACAGACAAAAAAAAGGCGACTCTTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                              <  1:188221/251‑1 (MQ=255)                                                                             cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                                                                                                       >  2:104414/1‑251 (MQ=255)                                                                                     ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                                                                                                               >  2:94338/1‑251 (MQ=255)                                                                                                    cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                                 >  1:273924/1‑250 (MQ=255)                                                                                                      aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                                                              >  2:102127/1‑251 (MQ=255)                                                                                                          cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                                                                                                          <  1:300919/251‑1 (MQ=255)                                                                                                              aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTTCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                                                                      <  2:300574/251‑1 (MQ=255)                                                                                                                         gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAACAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                                                                           <  2:409767/251‑1 (MQ=255)                                                                                                                                  tCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                   <  2:336176/250‑1 (MQ=255)                                                                                                                                   cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                                                                 >  2:381376/1‑251 (MQ=255)                                                                                                                                     tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCCCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGGGTCTTCa                                                                                                                >  2:162060/1‑250 (MQ=255)                                                                                                                                      gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGACAGAGTCTTCAAc                                                                                                              >  1:155725/1‑251 (MQ=255)                                                                                                                                        cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                            >  2:78602/1‑251 (MQ=255)                                                                                                                                          tCGCTCTTCATCTTACGTTTCTACAGACAACAAAAAGGCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAActt                                                                                                            <  2:319862/249‑3 (MQ=255)                                                                                                                                              ccctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGACAGAGTCTGCAACGTATTCCa                                                                                                      <  2:155725/250‑1 (MQ=255)                                                                                                                                               cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGACTTAAAAACCAGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGATTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCGGTCTCTTCGGCAGAGTCTTCAACGTATTcc                                                                                                       <  2:299709/249‑1 (MQ=255)                                                                                                                                                              cTACAGACAAAAAAAAGGCGACTCATCAATCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCTCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGACGAGTAGCTCTTCATGGGTTAACGGCTCTCCTTCTGCCGTCCCCTGCACTCCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                                                                                                               <  1:440637/194‑1 (MQ=255)                                                                                                                                                                  aGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATTTCGGTCTCGGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                                                                  <  1:185926/251‑1 (MQ=255)                                                                                                                                                                  aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                                                                  <  2:326265/251‑1 (MQ=255)                                                                                                                                                                    aCAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCTCTCCAAGGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                                                                <  2:335206/251‑1 (MQ=255)                                                                                                                                                                                    ccatcaGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGTTTCGA‑TGC‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGTTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGTTCGCACACCGCGTCGCCGTTCACAATTCCGACAGaa                                                                <  2:293980/250‑1 (MQ=255)                                                                                                                                                                                     catcaGTCGCTCTAAAAAGCAATTTGCCAGCGCCGCAGTCGGCGTCGCTACCATGTAGTCGA‑T‑T‑CCCCCGCGGGCCAGAGCCATCAAGATGAGGAGCTCTTCATGGGTGATCGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATTCGCCCGTCTGCGTGCCTCTCTACGTTCATTTCTGTCTCGGCTTCAGAGTTTTCAACGGAATCCATATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaaa                                                               <  2:228346/251‑1 (MQ=255)                                                                                                                                                                                     catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaa                                                                >  2:46957/1‑250 (MQ=255)                                                                                                                                                                                         aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTAGCTACAATGGATTCAA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc                                                           >  1:351988/1‑251 (MQ=255)                                                                                                                                                                                             gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCGTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCGTGTCGGTCTCTGCGGCAGAGTCTTCAACGGATTCCAGATCGCAAACTTCGTCGCCGTTCACAATTCCGACAGAAGCTGCGGCt                                                       <  2:209562/251‑1 (MQ=255)                                                                                                                                                                                                    aaagcAGTTTGCCGTCGCCGCCTTCTGCGTCGCTACAAGTGATACGA‑T‑T‑CCCCCGCGGGCCAGAGCCAACAAGAGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCTGTGCGCGGTTCACAGTTCCGGCATAAACTGCGGCTTCCATcc                                                <  2:57816/247‑1 (MQ=255)                                                                                                                                                                                                    aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc                                                >  1:270111/1‑251 (MQ=255)                                                                                                                                                                                                          gttggCCAGCGCCGCCTTCTGCGTCGCTACAATGGAGTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATGCCTTTCa                                          <  1:322586/247‑1 (MQ=255)                                                                                                                                                                                                                       cccTTCTGCGACGCTACCATGGATTCGAGT‑C‑CCCC‑GCGGTCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCTGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAGCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                             <  2:7031/250‑1 (MQ=255)                                                                                                                                                                                                                             tGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtt                       >  1:15637/1‑251 (MQ=255)                                                                                                                                                                                                                               cGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                     >  1:195546/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTTAGCTTGCCg                   >  1:373257/1‑251 (MQ=255)                                                                                                                                                                                                                                    cTACACTGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGCGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTAATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGttt                >  2:413624/1‑251 (MQ=255)                                                                                                                                                                                                                                       cAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTACGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcc             <  1:368544/251‑1 (MQ=255)                                                                                                                                                                                                                                            gattcga‑t‑t‑cCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg        <  2:200874/251‑1 (MQ=255)                                                                                                                                                                                                                                             attcga‑t‑taCCCC‑TCGGGAACGAGCCAACAAGATGAGTAGCTCTTAATGGGTGACCGGCTCGCCTTCTGCAGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGGCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCATTTCATCGGATTGGTTTTCCGCTTGCCGtttgcccacagc       >  2:298225/1‑242 (MQ=255)                                                                                                                                                                                                                                             attcga‑t‑t‑cCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGAACTTCAATGATGCGCCCGTCTTCGGTCATCACTCAGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGTTCACAATTCCGACAGAACCTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccaccaac       >  1:399822/1‑249 (MQ=255)                                                                                                                                                                                                                                                   gattccccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCt   >  1:150246/6‑251 (MQ=255)                                                                                                                                                                                                                                                    attccccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTg  >  1:338459/5‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTG  >  NC_000913/3815636‑3816126

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