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	95.0%	177.5 / ‑0.6	61	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (35/22);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTT  >  NC_000913/3815633‑3816123                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  ggcggcccccGGTATTAAACAGCCCGGCGTTAAAAAAATAGGGGCTTTTCCGCCCGGATTTCAGAGTAACCTCGCCAACCTTTAACCCCTGCTGGCTAAGCGCAAATGAAAAAACCTGGCGCTTATAGGGTTCCAGGCCTCCCCCCCTCACCTTACTTTCCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTc‑ccc                                                                                                                                                                                                                                                    <  2:214732/244‑1 (MQ=255)                     aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCcaaca                                                                                                                                                                                                                                 >  1:35436/1‑249 (MQ=255)                                gaagaaATAGGTGCTTTTTCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGct                                                                                                                                                                                                                     <  2:446974/250‑1 (MQ=255)                                        aGGGGCTTTTGCGCCCGGNTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                         <  2:238927/238‑1 (MQ=255)                                        aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                         >  1:238927/1‑238 (MQ=255)                                        aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCCTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATgg                                                                                                                                                                                                             <  1:183833/250‑1 (MQ=255)                                           ggCTTTTGCGCCCGGATTTCAGCGTAAACACNCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                         <  2:249792/251‑1 (MQ=255)                                             cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                       >  1:472566/1‑251 (MQ=255)                                               tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                      <  1:428267/250‑1 (MQ=255)                                               tttGCGACCGGATTTCCGCGTAAACTCGCAAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACgg                                                                                                                                                                                                     >  1:429917/1‑251 (MQ=255)                                                 ttcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCAC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                    <  2:35436/248‑1 (MQ=255)                                                 tGCGCCCGTATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAGAGGCGGCTCATCAGGCGCATTAAAAAGCAGTTTGCCAGCGCCGCCTTGTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGGCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                    <  2:143015/250‑1 (MQ=255)                                                                  gTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                  >  1:91726/1‑251 (MQ=255)                                                                  gTAAACCCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                  <  2:91726/251‑1 (MQ=255)                                                                        ccgcgtaacgttaccacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt                                                                                                                                                                            <  1:253298/237‑1 (MQ=255)                                                                          gCCAAACTTTAACACCTGCTTGCTAAGCGCCAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACACGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                          >  2:480647/1‑251 (MQ=255)                                                                            cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTAGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGTGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCa                                                                                                                                                                         >  2:402158/1‑250 (MQ=255)                                                                                   tAACACCTGCTTGCTAAGCGCAAATTCAATTAACGGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAACGTACTTTTCGACAGACAAAAAAAAGGCGACACATCAGTCGCCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCGCCAATGTATTCGAGTCC‑CCC‑TC‑GGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTGCAATGATgc                                                                                                                                                                  <  2:335057/250‑1 (MQ=255)                                                                                    aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                                                                                                                                                                <  1:372187/251‑1 (MQ=255)                                                                                               gcaacgcgaaaattcaagatccttgcgcTGAAATGTTTTCATGCCTTCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACGCAGCAGTCTCCGTAAAAATCAGTTTGCCATCGCCGCCATCTGCTTCGCTACAATGGAGTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTATCTCTTCATGGGTGAACGGCTCGCCTTCTGTCGTCCCCTGCACTTCTATGATGCGCCCGTCTTCGGt                                                                                                                                                     <  2:282294/227‑1 (MQ=255)                                                                                                                 aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTAAAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATAACTACGTTCATGTc                                                                                                                                   >  2:178271/1‑251 (MQ=255)                                                                                                                 aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCCCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                                                                   <  1:276112/251‑1 (MQ=255)                                                                                                                      ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                                                              >  1:336019/1‑251 (MQ=255)                                                                                                                         gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                                                           >  1:480168/1‑251 (MQ=255)                                                                                                                           tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                         >  1:267505/1‑251 (MQ=255)                                                                                                                                  gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGt                                                                                                                  >  1:322144/1‑251 (MQ=255)                                                                                                                                     tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCCCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCtt                                                                                                               >  1:260646/1‑251 (MQ=255)                                                                                                                                           cTTCGCTACTCATCTTACTTTTCTACAGACAACAAACAGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGAGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGACAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGGCTTCAACgg                                                                                                         >  2:206813/1‑250 (MQ=255)                                                                                                                                              cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTAtt                                                                                                      >  1:414893/1‑251 (MQ=255)                                                                                                                                                     cATCTTACTTTTCTACAGACAAAAAAAAGNCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCGTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACTTATTCCAGa                                                                                                 <  2:249094/249‑1 (MQ=255)                                                                                                                                                              tttCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatt                                                                                        <  1:117613/249‑3 (MQ=255)                                                                                                                                                               ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatttct                                                                                     >  2:322001/1‑246 (MQ=255)                                                                                                                                                                 cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                   >  2:511074/1‑244 (MQ=255)                                                                                                                                                                   aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttcac                                                                                 <  1:389571/251‑10 (MQ=255)                                                                                                                                                                   aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                   <  1:471834/249‑8 (MQ=255)                                                                                                                                                                      gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                              >  2:182970/1‑251 (MQ=255)                                                                                                                                                                      gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                              >  1:155674/1‑251 (MQ=255)                                                                                                                                                                         aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGGCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGttt                                                                            >  2:422116/1‑249 (MQ=255)                                                                                                                                                                          aaacaaaGGCGCCGCATCAGCCGCCTTAAAAATCAGTTGGCCAGCGCCGTCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCGC‑GGGCCAGCTCCAACAAGGTGAGTAGCTCTTCATGGGTTCACGGCTCGCCGTCTGCCATCCCCTGCACTGCAATTATGCGCCCGGCTTCGGTCATCACTCCGTGCATGTCGGTCTCTGCGGCTTATTCTTCAACGGATTCCAGATCGCAAAGTGCTTCACCGtttac                                                                          <  2:210915/247‑4 (MQ=255)                                                                                                                                                                              aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTCACCC‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCAATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAt                                                                       >  1:183388/1‑250 (MQ=255)                                                                                                                                                                               aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtac                                                                     >  1:262752/1‑249 (MQ=255)                                                                                                                                                                                 ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                   <  1:440499/251‑1 (MQ=255)                                                                                                                                                                                 ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                   <  1:344055/251‑1 (MQ=255)                                                                                                                                                                                   cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                 >  1:197973/1‑251 (MQ=255)                                                                                                                                                                                    gaatcataaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCATC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                 >  1:293667/4‑250 (MQ=255)                                                                                                                                                                                             gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCg                                                       >  1:510088/1‑251 (MQ=255)                                                                                                                                                                                               cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGc                                                     >  2:384908/1‑251 (MQ=255)                                                                                                                                                                                               cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                         >  2:229719/1‑183 (MQ=255)                                                                                                                                                                                               cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                         <  1:229719/183‑1 (MQ=255)                                                                                                                                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAATCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTaac                                                 >  2:303507/1‑249 (MQ=255)                                                                                                                                                                                                     aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGCTTACAATACCGACAGACACTGCGGCTAccct                                               >  2:142792/1‑249 (MQ=255)                                                                                                                                                                                                       aaaTCAGTTTGCAAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCACAAAGGCGAGTAGCTCTTGATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCCCTTTAATTGTGCCCCCATCGATGGTGTTCACTACCTTCATGTCGGTCTATTCTGCTGACTCTTCAATGTATCGCAGGTCCCAAATTGCTTCAGCGTTNACAATCCCGCCAGAAACTGCGGccaccaaac                                             >  2:255360/1‑248 (MQ=255)                                                                                                                                                                                                                gtcccgcgcCGCATTCTGCGTCGCTACAATGTATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCAAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCTAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCCTTCATCg                                    <  2:488228/246‑1 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGGCTTCACCGTATTCCAGATCGCAACTTGCTTCCCCGTTTACACTACCTACAGAAACTGCGGCTACCACCCCTTTCCTCgg                                   >  2:66177/1‑251 (MQ=255)                                                                                                                                                                                                                    aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAAAGAAACTGCCGCCACCATCCCTTTCATCggggt                                >  2:285705/1‑248 (MQ=255)                                                                                                                                                                                                                    aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCgggtt                                <  1:322001/251‑4 (MQ=255)                                                                                                                                                                                                                        ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCCTCCCTTTTATCCggttggtt                            >  2:37546/1‑251 (MQ=255)                                                                                                                                                                                                                        ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCGggttggtt                            >  1:308149/1‑251 (MQ=255)                                                                                                                                                                                                                               cTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                     >  2:378524/1‑251 (MQ=255)                                                                                                                                                                                                                                         acaagagctgcgatgcc‑ccc‑TC‑GGGCCAGCGCCAACATGGGTAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTTTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCTTCCCTTTCAACGGGTTGGTTTTCAGCTTGCCGGTTTc           <  2:366404/235‑1 (MQ=255)                                                                                                                                                                                                                                         aCAAGGGATTCGA‑TTC‑CCC‑CCCGGGCCAGCGCCAACCAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCGCTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATATCGACAGAAACTGCGGCTACAATCCCTTTCATCGGGTTGGTGTTCAGCTTGCCGtttt            <  2:193/250‑1 (MQ=255)                                                                                                                                                                                                                                          cAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTTCCGTTTTcc          >  2:19336/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tcga‑ttc‑CCCCTC‑GGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGACACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCtt  >  2:64700/1‑251 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTT  >  NC_000913/3815633‑3816123

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