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	91.7%	121.3 / 7.5	49	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base C (31/13);  total (32/17)
Fisher's exact test for biased strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCG  >  NC_000913/3815633‑3816131                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          gATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCCCTACAATGGATTCGA‑T‑Tcccc                                                                                                                                                                                                                                                            >  2:483654/1‑251 (MQ=255)    cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑Tccccctc                                                                                                                                                                                                                                                         >  1:218465/1‑248 (MQ=255)                   aCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAAGTAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGAGT‑CCCCC‑TC‑GGGCCAGAGCcaaca                                                                                                                                                                                                                                         <  1:384161/251‑1 (MQ=255)                       cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTAAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAgct                                                                                                                                                                                                                                     >  2:257620/1‑249 (MQ=255)                            cGTTGAAGAAAGAGGGGCTTTTGCGCCCGGATTTCAGCTTAAACTCGCCAAACTTTAACACCTGCGTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCGTGCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGA‑TGCCCCC‑TC‑GGTCCAGAGCCAACAAGATGAGTa                                                                                                                                                                                                                                <  1:490696/251‑1 (MQ=255)                                               tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                              <  2:68246/250‑1 (MQ=255)                                               gtgtcgaccGGATTTCAGCGTAACCTCGCCGAACTTTAACCCCGGCTTGCTAAGCGCAAATTTAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCCCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAAGAAGGTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGAGTCGA‑T‑TCCCCCTC‑GGGCAAGAGCCAACAAGATGAGTAGCTCTTCATGGTTGAAc                                                                                                                                                                                                               <  2:109349/242‑1 (MQ=255)                                                  tcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCCTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                                                                                                                                           <  1:286120/249‑1 (MQ=255)                                                  gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGCCTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCGTCAGGGGTGAACGtctc                                                                                                                                                                                                          >  2:188293/1‑251 (MQ=255)                                                     cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGCGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                                                                                                                                       >  2:387940/1‑251 (MQ=255)                                                      ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                                                                                                                                                                                                      >  1:27354/1‑251 (MQ=255)                                                         gATTTCAGCGTAAACGCGCCAAACTTTAATACCTCCTTGCGATGCGCAAATTCAATAATCTGGCGCTGATATGGTTTCATGACGTCGCTCCTCACCTTACTTTGCTACAGCCAACAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGTCTTTTGCGGCGCTAGAGGGGATTCGA‑T‑TCCCCCGC‑GTGCGAGATCGAACAATATGAGTATCTCATCAGGGGTGAGCGGCTCGCCGTCt                                                                                                                                                                                                   <  2:306748/251‑1 (MQ=255)                                                                               aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                                                                                                                                                                              >  1:484523/1‑250 (MQ=255)                                                                                cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                                                                                                            >  1:360201/1‑251 (MQ=255)                                                                                                  aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCCCGCCTTCTGCCGTCCCCTTCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                          >  2:293904/1‑251 (MQ=255)                                                                                                         aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                                                                   >  2:205095/1‑251 (MQ=255)                                                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                                 >  1:462524/1‑251 (MQ=255)                                                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTATCTGCGGCAAAGTCttacacgt                                                                                                                 >  1:272961/1‑245 (MQ=255)                                                                                                                                                ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTcc                                                                                                            >  2:390709/1‑251 (MQ=255)                                                                                                                                                              tttCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg                                                                                               <  2:308388/250‑1 (MQ=255)                                                                                                                                                              tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCACACCGc                                                                                              >  1:218165/1‑251 (MQ=255)                                                                                                                                                                 cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                                                                           >  1:423353/1‑251 (MQ=255)                                                                                                                                                                 cTACAGACAAAAAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAGCGCCGACTTCTGCGTAGCTACCATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGACTTCGGTCATCACTCAGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                                                                           >  1:302955/1‑251 (MQ=255)                                                                                                                                                                   aCAGACAAAAAAAAGGCGACTCATCAGTCGCCATAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATGCGAGT‑CCCCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCATGTCGGGCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACTGCTTc                                                                                           <  1:171348/249‑1 (MQ=255)                                                                                                                                                                          aaaaaacGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCGGCGGAGTCTGCAACGTAGTCCAGATCGCAAACCGCTTCGCCTTTcac                                                                                  <  2:470884/251‑1 (MQ=255)                                                                                                                                                                                  gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCAGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCAGCTTCGCCGTTcacactcacg                                                                           >  2:491990/1‑244 (MQ=255)                                                                                                                                                                                       tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGAAGCTCTTCATGGGTGAACGGCGCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCGGTCTCTGCGGCCGCGTCTTCCACGTATTCCAGATTGCAAACTGCTTCGCCGTTccacattccgccaga                                                                      >  2:187249/1‑236 (MQ=255)                                                                                                                                                                                            aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGATCTTCATGGGTGAAAGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGAGGCAGAGGCTTCAACGTATTCCAACGCGACAACCGCTTCGCCGTTCAAAAATCCGACAGAAACtcc                                                                >  2:283870/1‑249 (MQ=255)                                                                                                                                                                                              tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                              >  1:220889/1‑251 (MQ=255)                                                                                                                                                                                                gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGC‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATTCGCCCGTCTTCGGTCCTCCCTACGGTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACATATCCGACAGAAACTGCGGCt                                                            >  2:88822/1‑251 (MQ=255)                                                                                                                                                                                                gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                            <  1:50655/251‑1 (MQ=255)                                                                                                                                                                                                gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                            <  1:483654/251‑1 (MQ=255)                                                                                                                                                                                                 ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                            >  1:285512/1‑250 (MQ=255)                                                                                                                                                                                                  cgtAAAAATCAGGTGCCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGA‑T‑TCCCC‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCGTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGGATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc                                                          <  2:286981/249‑1 (MQ=255)                                                                                                                                                                                                     aaaaatcatatagCCAGCGCCGCCTTCTGCGTCGCTACAATGGATCCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTTAACGGCTCGCCTTCGGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCGTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACATTTCCGACAGAACCTTCGGCTACCa                                                        <  2:126527/238‑1 (MQ=255)                                                                                                                                                                                                     aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt                                                       >  1:17739/1‑251 (MQ=255)                                                                                                                                                                                                              tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAt                                              >  1:194426/1‑251 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTGCTCCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTGTCAt                                              <  1:470831/248‑1 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCgg                                           <  1:368536/251‑1 (MQ=255)                                                                                                                                                                                                                      cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                      >  1:83631/1‑251 (MQ=255)                                                                                                                                                                                                                       gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                                     >  1:387904/1‑251 (MQ=255)                                                                                                                                                                                                                        ccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCTGCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtt                                    <  1:468055/251‑1 (MQ=255)                                                                                                                                                                                                                         cgccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGACAGAAACTGCGGACACCATCCCTTTCATCGGATTGGttt                                   >  2:374849/1‑251 (MQ=255)                                                                                                                                                                                                                          gccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                  >  1:202941/1‑251 (MQ=255)                                                                                                                                                                                                                                   gTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                         >  1:345082/1‑251 (MQ=255)                                                                                                                                                                                                                                   gTCGCTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                         >  1:165220/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTACAATGGATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCCACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGAAGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAATGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGATTGGGTTTCAGCGTGCCGttt                     >  2:244591/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTGGATTCGA‑T‑TCCCCAGC‑GGGCAAGAGCCAACAAGATGAGTAGATCTTCATGAGTGAACGGCTCGCCTTCTGACGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTAAACGTATTCCAGATCTCAAAACGCTTCGCCGTTCAAAATTCAGACAGAAACTGCGGCTACCCTCCCTTTCATCGACTTGGTTTTCAGCTTGACGttttaca                 >  2:446403/1‑247 (MQ=255)                                                                                                                                                                                                                                              ggATTCGA‑T‑TCCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacca              <  1:435696/251‑1 (MQ=255)                                                                                                                                                                                                                                               gattcga‑t‑tcCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTACGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacc               >  2:25742/1‑249 (MQ=255)                                                                                                                                                                                                                                                 ttcga‑t‑tcCCCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCCGCTTGCCGTTTTccacccgcc           >  2:359860/1‑247 (MQ=255)                                                                                                                                                                                                                                                       ttccccc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCCGCTTCTGt      >  2:411785/4‑251 (MQ=255)                                                                                                                                                                                                                                                        tccccc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt      >  2:430935/3‑250 (MQ=255)                                                                                                                                                                                                                                                        tccccc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt      >  1:486742/3‑250 (MQ=255)                                                                                                                                                                                                                                                           ccc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcg  >  1:272372/1‑251 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCG  >  NC_000913/3815633‑3816131

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