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	89.7%	147.5 / 11.8	58	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base C (29/23);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGC  >  NC_000913/3815633‑3816130                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         gATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCAAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGAGCTGATATGGTTTCATGCCTTCGCTCCTCATCATACTTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑cc‑c‑c                                                                                                                                                                                                                                                           >  1:514959/1‑251 (MQ=255)  aTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGCAATAGGGGCTTTTGAGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑cc‑c‑cc                                                                                                                                                                                                                                                          >  2:360974/1‑250 (MQ=255)    cgcgCCAGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑cc‑c‑cctc                                                                                                                                                                                                                                                        >  2:378428/1‑248 (MQ=255)        cccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGc                                                                                                                                                                                                                                                   >  2:17828/1‑251 (MQ=255)           ggTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGcc                                                                                                                                                                                                                                                  >  1:422839/1‑249 (MQ=255)              aTTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCagag                                                                                                                                                                                                                                              <  1:177817/250‑1 (MQ=255)                     aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCcaaca                                                                                                                                                                                                                                        >  1:385349/1‑249 (MQ=255)                      gCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCTCTGATATGGTTTCAGGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCATGCGTCTTAAAAATCAGTTTGCCAGCGGCGCCGTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCGC‑GTGCCAGAGCCAACAAGa                                                                                                                                                                                                                                     <  2:445298/251‑1 (MQ=255)                            cGTTGAAGAAAAAGGGGCTTTTGCGACCGGATTTCAGCGTACACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCACATTCAATAAACTGGCGCTGATATGGTTTCATACCTTCGCTCCTCATCTTACTTTTCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                >  2:457164/1‑250 (MQ=255)                                gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGct                                                                                                                                                                                                                            <  1:131758/250‑1 (MQ=255)                                   gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                                                                                                                                                                                        >  2:484086/1‑251 (MQ=255)                                               tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                             <  1:326030/250‑1 (MQ=255)                                                            ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                                <  2:135522/250‑1 (MQ=255)                                                             tCAGCGAAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCCCGCCTTCTGCCg                                                                                                                                                                                              >  2:33737/1‑251 (MQ=255)                                                                  gTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                         >  2:25771/1‑251 (MQ=255)                                                                   gaaactcaacAAACTTTACCTCCTGCTTGCTAAGCGCAAATTCAATAAACTGTGGCGGATATTGTCTCTGGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCGCCTCATCAGTCTCTTTAAAAATCATTTTGCCAGCCCCGCCTTCTGCGTCGCTTCAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTTAACGGCTCGCCCTCTGCCGTCCCCt                                                                                                                                                                                        <  2:363418/242‑1 (MQ=255)                                                                    aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                                                                                                                                       <  1:114146/251‑1 (MQ=255)                                                                            cAAACTTTAACACCTGCTTGCTAAGCGCAAATTAAATAAACTGGCGATGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACAATGGATTCGAATC‑CC‑C‑C‑TC‑GGGCCAGAGCACACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGCCCACTGCACtgca                                                                                                                                                                                >  2:452811/1‑247 (MQ=255)                                                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                              >  1:164834/1‑251 (MQ=255)                                                                                                  aatttcAAATTCAAGAAAGTGGTTCTGATATGTTTTCATTCATTCGCTCCTCATCTTACTTTTCTTCAGACAAAAAAAAGGCGTCTCATCAGTAGTCTAAAAAAGCAGTTAGCCAGCGGAGCCTTCTGCGTCGCTACAATGTATTCGA‑TC‑CC‑CTC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTGCATGGGTGAACGGCTGGTCTTCTGCACTCCCCTTCACGTCAATGATGTGGCTGTCTTCGTtca                                                                                                                                                          <  2:22710/245‑1 (MQ=255)                                                                                                            tCAATAAACTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAAAACAGGCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTTATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                                                                <  2:57097/250‑1 (MQ=255)                                                                                                             cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCACCCGTCTTCGGTCATCACTCCGTTc                                                                                                                                               >  2:381888/1‑250 (MQ=255)                                                                                                              aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTTCAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCGCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                             <  1:465952/251‑1 (MQ=255)                                                                                                                      ggCGCTGATTTGGTTTCATGCCTTCGCTCGTCATCTTACTATTCTCCATACTAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCCTCGCCGCCTACTGCGTAGCTACACTGGCTTCGA‑TT‑CC‑C‑CCTC‑GGGCCTGAGCCAACAAGGTGGGTAGCTCTTCAGGGGTTAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGttactctcggtct                                                                                                                                     >  2:348269/1‑240 (MQ=255)                                                                                                                                  gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                                                                                         >  1:239690/1‑251 (MQ=255)                                                                                                                                       aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                                                                                                                    <  1:236314/251‑1 (MQ=255)                                                                                                                                           cTTCGCTACTCATCTTACTTTTATACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGACTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGt                                                                                                                >  2:331563/1‑251 (MQ=255)                                                                                                                                                      aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                           >  1:339311/1‑245 (MQ=255)                                                                                                                                                           aCTTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAGTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                 <  2:403525/250‑1 (MQ=255)                                                                                                                                                                tCTACAGACACAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                                                                           >  2:236446/1‑251 (MQ=255)                                                                                                                                                                   aCAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAATTCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTGCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTTATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGc                                                                                        <  1:255951/251‑1 (MQ=255)                                                                                                                                                                          aaaaaaagaagaataatcaGTCGACTTACAAATCAGTTTACCAGAGACGACTTCTGCGTCGCTACAATGGAATCGA‑TT‑CC‑C‑CCTC‑AGGCCAGAGGCAACAAGATGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATAACTACGTTCATGTCGGTCTCTGAGGCAGAGTCGTCAACGTATTCAAGATCGCAAACCGCTTCGCCGttacc                                                                                 >  2:174612/16‑248 (MQ=255)                                                                                                                                                                          aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACTGCTCTCCTTCTGCCTTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATTACTACGTTCATGTCGGTCTCTGCTGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                                 <  1:496357/251‑1 (MQ=255)                                                                                                                                                                              aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtt                                                                             >  1:227155/1‑251 (MQ=255)                                                                                                                                                                                aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcacactacc                                                                           >  2:426451/1‑246 (MQ=255)                                                                                                                                                                                   cGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGATCTTCATGGGGGAACGGATCGCCTTCTGCCGTCCCCTGCACTTCAATGAAGCGCCCTGCTACGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGAATTCCAGATCGCAAACCGCTTAGCCGTTACCAATTCCGAc                                                                        >  2:269176/1‑251 (MQ=255)                                                                                                                                                                                      ctcatccgacaccTTAAAAATCATTTTGCCCGCGCCGCCTTCTGCGTCGCTACAAAGGATTCGA‑TT‑CC‑C‑C‑CCCGGGCCAGAGCCAACAAGATGCGTAGCTCTTCATGGGTGAACGGCCCGCCTTCCGCCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGTCATCACCACGTTCCTGTCGGTCTCTGCGGCAGAGGCTTCCACGTCTTCCAGATCGCCAACAGCTCCGCCGTTCACAGTTCCGAc                                                                        >  2:221451/12‑248 (MQ=255)                                                                                                                                                                                      cTCATCAGTCGCCTTCAAAATCAGATTGCCCACGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGACATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGCTCACAATTCCGAc                                                                        >  2:311886/1‑248 (MQ=255)                                                                                                                                                                                      cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                        >  1:370416/1‑248 (MQ=255)                                                                                                                                                                                        catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACCAGGGATTCGA‑TTACC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGATCTTCATGGGAGAACGGCTCGCATTATGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAACGCTTCGCCGTTCACACTTCCGACAgca                                                                    >  1:24446/1‑248 (MQ=255)                                                                                                                                                                                           caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                                                                <  2:108452/251‑1 (MQ=255)                                                                                                                                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAcc                                                        >  1:160653/1‑251 (MQ=255)                                                                                                                                                                                                     aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCAGCACTCCGTTCATGTCGGTCTCTGCTGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt                                                      <  1:66686/251‑1 (MQ=255)                                                                                                                                                                                                           cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCttt                                                <  2:254010/251‑1 (MQ=255)                                                                                                                                                                                                           cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt                                                 <  1:484086/250‑1 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                                                                          <  2:271252/187‑1 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                                                                          >  1:271252/1‑187 (MQ=255)                                                                                                                                                                                                                 gCCAGCGCCGCCTTCTGCGTCGCTACAAGGGATTCGA‑TT‑CCGC‑C‑GC‑GGGCCAGAGCCAACGAGATGAGTAGCTTTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACGTCGGTCATGGCGGTCTCTGCGTCGGAGGCTTCAACGTATTCCATTTCGCAACCCGCTGCTTCGTTGACAATGACGACAGTAACTGGGGCTACCTTCGCTGTCATCgg                                          <  2:216052/251‑1 (MQ=255)                                                                                                                                                                                                                      cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCATATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCTGATTgg                                     <  1:406988/251‑1 (MQ=255)                                                                                                                                                                                                                       gccgccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                                    >  1:422081/1‑251 (MQ=255)                                                                                                                                                                                                                         cgccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                                  >  1:274063/1‑251 (MQ=255)                                                                                                                                                                                                                          gccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                 <  2:385349/251‑1 (MQ=255)                                                                                                                                                                                                                               cTGCGTCGCTACAATGGATTCGC‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                            >  2:185044/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                         >  2:413567/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                         <  1:465863/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTCGCTACAATGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTCGCTCTTCATAGGTGAACGGCTCGCCTTATGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTACATGTCGGTCTCTGCGGCAGAGTCTTCACCGTCTTCCAGCTCGCAACCCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCCTCCCTTTCATCGGATTGGTTTTCAGCTTGcc                        >  1:205722/1‑251 (MQ=255)                                                                                                                                                                                                                                        tACAATGGATTCGA‑TT‑CCAC‑C‑TC‑GGGCCAGGGGCAACAAGATGAATCGCTCATCATGGGAGAACGGCTCGCCTTCTGCAGACCCCTGCACATCAATGATGCGCCCGTAATCGGTCATCACTACGTACCTGACGGTCTCTGCGGAAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGtttt                   >  2:311733/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTGGATTCGA‑TT‑CC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGCAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccac               >  1:278649/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ttaggttccc‑c‑c‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTATTCATGCGTGAACGGCTCGCCATCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTAATCACTACGGTCATTTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCCCCGTTCGTAATTCCGACAGAACCTGCGGCTACCATCCCTTTCATACGATTGGGTTTCAACTTGCCGCTTTcccccaact          >  2:19147/9‑248 (MQ=255)                                                                                                                                                                                                                                                  tcga‑tt‑cC‑C‑CCTC‑GTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGTCAGAGTCTTCACCGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCtt         >  2:5991/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tcga‑tt‑cC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTGCCATCCCTTTCAGCGGCTTGGTTTTCAGCTTGCCGTTTTCCACCAGCt          <  2:423035/250‑1 (MQ=255)                                                                                                                                                                                                                                                   cga‑tt‑cC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTc        <  2:126041/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cgagtccc‑c‑c‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTCATCACTGCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCTCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTc        <  1:494933/245‑1 (MQ=255)                                                                                                                                                                                                                                                          cc‑c‑c‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTGTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgc  <  2:357549/251‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGC  >  NC_000913/3815633‑3816130

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