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	88.7%	149.3 / 11.7	53	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base C (33/14);  total (37/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCAC  >  NC_000913/3815643‑3816117                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            ggTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCag                                                                                                                                                                                                                                   <  2:1595/251‑1 (MQ=255)          cAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCcaaca                                                                                                                                                                                                                           >  2:93245/1‑250 (MQ=255)                  cGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTGGCCAGCGCCGCCTTCTGCGGCTCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTa                                                                                                                                                                                                                  <  1:302824/251‑1 (MQ=255)                            aTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTGAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGGCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                                                                                                                                                                        <  1:23103/251‑1 (MQ=255)                                 ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                   <  1:8160/251‑1 (MQ=255)                                          gcCCGGCTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTACGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCAAGCGCCGCATTCTGCGTCGCTCCAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGCGAACGGCTCGc                                                                                                                                                                                          >  1:235785/1‑251 (MQ=255)                                              ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                                                                                                                                      >  1:70354/1‑251 (MQ=255)                                                                 ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCCATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                    >  2:58883/1‑250 (MQ=255)                                                                       tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGTTATAGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCAC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                             >  2:332047/1‑251 (MQ=255)                                                                        ttAACCCCAGCTAGCTAACCCCAACTTCAATACACTAGCGCTGATATGGATTCATGCCTTCGCTCCTCATCTTACTTTTCTCAAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAATAATCAGTTTGCAAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CAC‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCACTGAACttacatgattc                                                                                                                                                            >  2:228759/1‑242 (MQ=255)                                                                                 gcttgctACGCGCAAATTCAATAAACTGTCGCTCATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTTAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCt                                                                                                                                                    <  2:232763/250‑1 (MQ=255)                                                                                                    aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                >  1:145829/1‑251 (MQ=255)                                                                                                    aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                <  1:126406/251‑1 (MQ=255)                                                                                                     aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                                                                                               <  1:235209/251‑1 (MQ=255)                                                                                                         aCTGGCGCTGATATGGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                                                                                           >  2:41004/1‑251 (MQ=255)                                                                                                                    tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGcc                                                                                                                >  1:264877/1‑250 (MQ=255)                                                                                                                     atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTgagttc                                                                                                                >  2:134320/1‑245 (MQ=255)                                                                                                                            cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                                                        >  2:119588/1‑251 (MQ=255)                                                                                                                               gCCTTCGCTCCTCATCTTACTTTACTACAGACAAAAAAAAGTCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCGGCGTCGCTACAATGGATTCGATTC‑C‑CGC‑TC‑GGGCCAGAGCCAACAAGATTAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAAGAATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                                                      <  2:123127/250‑1 (MQ=255)                                                                                                                                   tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                    >  1:156907/1‑232 (MQ=255)                                                                                                                                   tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAAGGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                    <  2:156907/232‑1 (MQ=255)                                                                                                                                   tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTACCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGATTAGCTCTTCATGGGTGAACGACTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGACTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTAt                                                                                                 >  1:246446/1‑251 (MQ=255)                                                                                                                                     gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCCACGTATTc                                                                                               >  2:326464/1‑251 (MQ=255)                                                                                                                                       tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                              >  2:185480/1‑250 (MQ=255)                                                                                                                                                  cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                                                                  >  2:102221/1‑251 (MQ=255)                                                                                                                                                             aCAAAAAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTAGGTCATCACTACGTTCATGTCGGTCTCTGAGGCCGAGTCTTCAACGTATTCCAGCTCGCCAACAGCTTCGACGtt                                                                       >  2:180303/1‑251 (MQ=255)                                                                                                                                                               aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCAAACCGCTTCGCCGTTc                                                                      >  1:159889/1‑250 (MQ=255)                                                                                                                                                                 aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                    >  1:94544/1‑250 (MQ=255)                                                                                                                                                                     aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGATCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                                                               >  1:147655/1‑251 (MQ=255)                                                                                                                                                                         cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                           >  1:89184/1‑251 (MQ=255)                                                                                                                                                                          gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                           >  1:108414/1‑250 (MQ=255)                                                                                                                                                                          caCTCATCAGTCGCCTTAAAAATACGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGACTTCGATTC‑C‑C‑CCTC‑GGGCCTGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCCGATCGCACACCGCTTCTCCGTTCCCAATTCCgcc                                                           >  2:42049/2‑248 (MQ=255)                                                                                                                                                                             ccatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaa                                                       <  2:102033/250‑1 (MQ=255)                                                                                                                                                                             ccatcaGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GTGCCAGATCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTCCATGGCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCTACAGaa                                                       <  2:145829/250‑1 (MQ=255)                                                                                                                                                                               atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAAc                                                     >  2:93942/1‑251 (MQ=255)                                                                                                                                                                                 caGTCGCCTTCAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTACTCCAGATCGCAATCCGCTTCGCAGTTCACAATTCCGACAGaacctt                                                   >  2:1655/1‑247 (MQ=255)                                                                                                                                                                                   gTCGCCTTAAAACTCAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGCTTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg                                                 >  2:325103/1‑251 (MQ=255)                                                                                                                                                                                     cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                <  2:108414/250‑1 (MQ=255)                                                                                                                                                                                         ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATACGATTCAC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACCGCACCTGCGGCTAcc                                           >  2:227413/1‑251 (MQ=255)                                                                                                                                                                                          tacaaaTCAGTTTGCCAGCGCCACCGTCTGAGTAGCTCCAATAGATACGCTTC‑C‑C‑CCTC‑GGGCCAGAGCCCACAAGATGAGTAGCTCTTCATGGGTGACCGGCACGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAAGTTCATGTCGGTCTCTGCGGAAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACCGCGGCtgcc                                           >  2:300311/4‑247 (MQ=255)                                                                                                                                                                                             aaaTCAGTTTGCCAGCGCCGCCTTCCGCGTCGCTACAAGGGATTCGATTC‑C‑C‑C‑CCCGGTCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc                                       <  1:5122/251‑1 (MQ=255)                                                                                                                                                                                              aaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATccc                                      <  2:252855/251‑1 (MQ=255)                                                                                                                                                                                               aTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGCTCACAATTCCGACAGAAACTGCGGCTACCATCCCt                                     >  2:180263/1‑251 (MQ=255)                                                                                                                                                                                                  aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                                                         <  2:83896/196‑1 (MQ=255)                                                                                                                                                                                                  aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                                                         >  1:83896/1‑196 (MQ=255)                                                                                                                                                                                                    tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGAAACTGCGGCTACCATCCCTTTCAt                                >  1:314732/1‑251 (MQ=255)                                                                                                                                                                                                      tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                              >  2:89861/1‑251 (MQ=255)                                                                                                                                                                                                         cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAAAGAAACTGCGGCTACCATCCCTTTCCTCGGAt                           >  1:185236/1‑251 (MQ=255)                                                                                                                                                                                                          aGCGCCGCCTTCTGCGTCGCTACAATGCATCCGATTC‑CGC‑C‑TC‑GGGCCCGAGCCAACAAGATGAGTAGCTCTTCCTGGGAGCCCGGCTCGCCTTCTGCCGGCCCCTGCCCTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCGTGTCGGTCTCTGCGGCAGAGCCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGCAACTGCGGCTAACATCCCTTTCATCGGAtt                          >  2:122707/1‑251 (MQ=255)                                                                                                                                                                                                            cgccgcCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACAATCCCTTTCATCGGATTgg                        >  1:174214/1‑251 (MQ=255)                                                                                                                                                                                                                   ttCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                   >  2:314916/1‑249 (MQ=255)                                                                                                                                                                                                                           cGCTACAATGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGt         >  1:111387/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTGGATTCGATTC‑C‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccac  <  2:226109/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            GGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCAC  >  NC_000913/3815643‑3816117

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