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.4%	169.4 / 7.0	61	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (34/23);  total (36/25)
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.14e-01

ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815634‑3816132                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            aTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTCTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑Tcc‑ccc                                                                                                                                                                                                                                                             <  2:69811/251‑2 (MQ=255)   cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑Tcc‑ccct                                                                                                                                                                                                                                                            <  2:270711/250‑3 (MQ=255)         cGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGcc                                                                                                                                                                                                                                                     >  2:12518/1‑250 (MQ=255)          ggTATTAAACAGCCCGGCGGTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCGACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCag                                                                                                                                                                                                                                                   <  2:211331/251‑1 (MQ=255)            tATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCAGGATTTCAGCGTAAACTCGCCACACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAATGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGATACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCag                                                                                                                                                                                                                                                   >  2:445768/1‑249 (MQ=255)                  aCAGCCAGGCGTTGCAGAGATAGGGGATTATGCGAACGGCTTTAAGCGTAAACTAGCCAAACTTTACCACCTGCTTGCCAAGCGCAAATTCCATAAACTAGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTACACATCGGTTTGCCAGCCCCGCCTTCTGAGTCACTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCcaac                                                                                                                                                                                                                                            >  2:200833/1‑250 (MQ=255)                    aGCCCGGCGTTGAAGAAATAGGGGGTTTTGCGCCCGGATTTAAGCGTACACTCGCCAAACTTTACCACCTGCTTGCTAAGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGCT‑TCC‑CCCTC‑GGGCAAGAGCcacca                                                                                                                                                                                                                                           >  2:464705/1‑246 (MQ=255)                       ccGGCGTTGAAGAACTAGGGGCTTTTGCGCCCGGATTTCAGCGTAACCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                                      >  2:234232/1‑251 (MQ=255)                            gTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAg                                                                                                                                                                                                                                 <  1:5361/251‑1 (MQ=255)                                   aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                                                                                                                                                                                          <  1:476071/251‑1 (MQ=255)                                    aaTAGGGGCTTTTGCGCCCGGATTTCAGCGGAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                                                                                                                                                                                                         <  1:323518/251‑1 (MQ=255)                                         gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                                                                                                                                                    >  2:345435/1‑251 (MQ=255)                                            cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                 >  2:407450/1‑251 (MQ=255)                                                                 gTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                            >  1:37715/1‑251 (MQ=255)                                                                     aCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAATGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGc                                                                                                                                                                                         >  1:75509/1‑250 (MQ=255)                                                                       ccGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt                                                                                                                                                                                      <  1:143428/250‑1 (MQ=255)                                                                         gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATCGTTTCATGCCGTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGAT‑TCC‑CC‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                    <  2:321222/251‑1 (MQ=255)                                                                          ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCACC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCCCTTCATGGGTGAACGGCACGCCTTCTGACGTCCCCTGCACTTc                                                                                                                                                                                    >  1:345339/1‑250 (MQ=255)                                                                              aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTAGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATACGCT‑TCC‑CCCTC‑GGGCCAGAGACAACAGGGTGATGAGATCTACATGGGTGCACGGCTCGACTACTGCCGTCCCCTGCACCTCAatg                                                                                                                                                                                >  2:353988/1‑250 (MQ=255)                                                                                 ttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgc                                                                                                                                                                            >  2:400361/1‑251 (MQ=255)                                                                                   aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                                                                                                           <  1:12518/250‑1 (MQ=255)                                                                                   aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                                                                                                           <  1:431983/250‑1 (MQ=255)                                                                                     cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                                                                        <  1:150974/251‑1 (MQ=255)                                                                                        cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                                                                     >  2:91830/1‑251 (MQ=255)                                                                                           cttgctAAGCGCGAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                                  >  1:466530/1‑251 (MQ=255)                                                                                              gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CC‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                               >  2:490725/1‑251 (MQ=255)                                                                                                                       cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                                                                      <  1:295069/251‑1 (MQ=255)                                                                                                                              atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                                                                               <  1:413192/251‑1 (MQ=255)                                                                                                                                     cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑c                                                                                                                       <  1:400361/251‑1 (MQ=255)                                                                                                                                     cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGACCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑c                                                                                                                       >  1:5851/1‑251 (MQ=255)                                                                                                                                              gCTCCTCATCTTACTTGGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTc                                                                                                              <  2:330810/251‑1 (MQ=255)                                                                                                                                               ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTACAA‑GTATTcc                                                                                                             >  2:412698/1‑251 (MQ=255)                                                                                                                                                    cATCTTACTTTTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCATCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCCGAGCCAACAAGATGAGTAGCTCTTCATGTGTGAACGGCGCGTCGTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCTGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATc                                                                                                        <  1:64647/251‑1 (MQ=255)                                                                                                                                                    cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATGCCC‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGAt                                                                                                         <  2:281292/250‑1 (MQ=255)                                                                                                                                                    cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAgctc                                                                                                        >  2:147441/1‑248 (MQ=255)                                                                                                                                                     aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTcc                                                                                                             >  2:127235/1‑245 (MQ=255)                                                                                                                                                         tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCAGAGTCTT‑CAACGTATTCCAGATCg                                                                                                       <  2:492600/247‑1 (MQ=255)                                                                                                                                                               tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCACCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCtt                                                                                             >  2:104374/1‑251 (MQ=255)                                                                                                                                                                    aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGACGCCTTCTGCGTCGCTACAATGGATCAGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAACATGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTATGCCATCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCg                                                                                        >  1:335874/1‑251 (MQ=255)                                                                                                                                                                         aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                                   <  1:260902/251‑1 (MQ=255)                                                                                                                                                                               aggtgACTCATCAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCGC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                                                                             <  1:482006/247‑1 (MQ=255)                                                                                                                                                                               aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                                                                             <  1:127235/251‑1 (MQ=255)                                                                                                                                                                                          caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                                                                  >  1:53609/1‑251 (MQ=255)                                                                                                                                                                                              cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGACACTGCGGc                                                              >  2:289281/1‑251 (MQ=255)                                                                                                                                                                                               gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGACACTGCGGCt                                                             >  2:317519/1‑251 (MQ=255)                                                                                                                                                                                                ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                             >  1:173074/1‑250 (MQ=255)                                                                                                                                                                                                 cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc                                                           >  1:58030/1‑251 (MQ=255)                                                                                                                                                                                                   tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAcc                                                          >  1:129407/1‑250 (MQ=255)                                                                                                                                                                                                          cAGTTTGCCAGCGCCGTCTTCTGCGTCGCTACAATGGCTTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGTTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCACCCCGCTGCGCCGTTCACAATTCCGACCGAAACTCCGGCTACCATCTCttt                                                  >  2:123317/1‑251 (MQ=255)                                                                                                                                                                                                           aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATTTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGACAGAAACTGCGGCTACCATCCCTTTc                                                 >  2:383288/1‑251 (MQ=255)                                                                                                                                                                                                                         gcattCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGCGCCAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGCCAGAGTCTT‑CAACGTATTCCAGATCGCCAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCCTCCCTTTCATCGGCTTGGtttt                                   >  2:226322/4‑251 (MQ=255)                                                                                                                                                                                                                          ccTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                                  >  1:92240/1‑251 (MQ=255)                                                                                                                                                                                                                                tcGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGCCGTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCTGATTGGTTTTCAGc                               <  2:50262/247‑1 (MQ=255)                                                                                                                                                                                                                                gCGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtt                             <  1:55033/250‑1 (MQ=255)                                                                                                                                                                                                                                 cgtcgtTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGACCTTCTTCGGTCATCTCTACGTCCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAACCCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                                   <  1:409278/211‑1 (MQ=255)                                                                                                                                                                                                                                 cGTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                           >  1:95638/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGTACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGATTGGTTTTCAGCTTGcc                          >  2:333240/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTCGCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                          >  2:422330/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCTACAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtt                       >  1:165665/1‑251 (MQ=255)                                                                                                                                                                                                                                        aCAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATTTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtttt                     <  2:336439/250‑1 (MQ=255)                                                                                                                                                                                                                                        aCAATGGATTCGAT‑TCC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGttt                      <  2:56811/249‑1 (MQ=255)                                                                                                                                                                                                                                        aCAATGGATTCGAT‑TCC‑CCCGC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGGCCCCTGCACTTCAATGATGCTCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTG‑CAACGTATGCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtttt                     <  1:464705/250‑1 (MQ=255)                                                                                                                                                                                                                                              gattcgat‑tcC‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacc                >  1:457901/1‑249 (MQ=255)                                                                                                                                                                                                                                                         c‑cc‑TC‑GGGCCATAGCCAACAAGATGCGTAGCTCTTCATGGGAGAACGGCTCACCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTCATCACTAAGATCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCCACCCGCTTCGCCGTTCACAATTCAGACCGAAACTGCGGCTAACCTCCCTTTCCTCGGGTTGGTTTTCCGCTTGCCGTTTTCCAACAGCATCTGTAgcg   >  2:5789/1‑251 (MQ=255)                                                                                                                                                                                                                                                            c‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CCACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  2:204026/1‑250 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTT‑CAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815634‑3816132

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: ‑