Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	93.4%	169.2 / 4.0	61	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (37/20);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCC  >  NC_000913/3815634‑3816108                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  gtCGCGCCCGTTATTAATCAGCCCAGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATGTCGGCGTAAACTCCGCAAACTTTAACACCGGCTTGCTAAGCGCAATTTCAATAAACTGGCGCTGATATGTGTTCATGGCTTCGCTCCCCATCTTCCTTTGCTACAGACAAAAAAAAGGCGTCTCAGCAGTCGTCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGTAGTCGA‑G‑T‑ccccc                                                                                                                                                                                                                                   <  2:360322/250‑2 (MQ=255)              ttAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                      >  1:102553/1‑251 (MQ=255)              ttAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                      >  2:144854/1‑251 (MQ=255)               tAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                      >  2:421522/1‑250 (MQ=255)                         ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGCTTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAggtg                                                                                                                                                                                                             >  1:288956/1‑246 (MQ=255)                          gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                          <  1:224317/251‑1 (MQ=255)                                  gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGACAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTc                                                                                                                                                                                                  >  1:395653/1‑251 (MQ=255)                                       aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGCAGCTCTTCATggg                                                                                                                                                                                             >  2:349891/1‑251 (MQ=255)                                        ggggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATgggg                                                                                                                                                                                            >  1:175350/1‑250 (MQ=255)                                          ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                          <  1:296712/251‑1 (MQ=255)                                             ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                       >  1:417956/1‑251 (MQ=255)                                              tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                        <  2:102553/249‑1 (MQ=255)                                                   gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                                                                                                                                 >  1:384951/1‑251 (MQ=255)                                                      cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGNAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                              >  2:246107/1‑251 (MQ=255)                                                       ggATTTCAGCGTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                                                                                                                             <  2:262865/251‑1 (MQ=255)                                                        aaggtccggcgaaaCTCGACAAACTTTAACAGCCGCTCGCTGAGCGGAAAATCAATAATCTTGCGCTGATAGGGTGTCATGCCTTCGCTGCTTATCTTACCTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGGCTTAAAAATCTGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGCGCCAGCGCCAACAAGGTGGTTAGCTCTTCATGGGTGGACGTCGCGCCtact                                                                                                                                                                            <  2:385157/240‑4 (MQ=255)                                                         aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                           >  2:120919/1‑251 (MQ=255)                                                         aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                           <  1:48753/251‑1 (MQ=255)                                                                             aaCTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAatg                                                                                                                                                       <  2:175350/251‑1 (MQ=255)                                                                                             tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                         <  2:186771/249‑1 (MQ=255)                                                                                             tgctAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCAACAGCCAAAAAAAAGGCGACTAATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCTGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCGGCCGTCCCCTGCACTTCGATGATGCGCCCGTCTTCg                                                                                                                                        <  2:184425/250‑1 (MQ=255)                                                                                                 aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                      <  1:445219/248‑1 (MQ=255)                                                                                                    cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                <  1:349891/251‑1 (MQ=255)                                                                                                    cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                <  1:421522/251‑1 (MQ=255)                                                                                                     gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                               >  1:199131/1‑251 (MQ=255)                                                                                                       aaaTTCACTAAACAGGCGCTGATATGGTTTCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAACAAGGCGCATCCTCAGTCGCCTTAAAAATAAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAAACAGGTGAGAAGCCCTTAATGGGTGAACGGCCCGCCTTCAGCCGTCCCATGCACTCAAATGATGCCACCGTCTTCGGtcatccccac                                                                                                                             >  2:436232/1‑246 (MQ=255)                                                                                                            cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                         >  1:407867/1‑250 (MQ=255)                                                                                                                 aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGCCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCATTCTGCCGTCCCCTGCACTTCAATGAAGCGCCCGTCTTACGTCATCCCTACGTTCATGGCgg                                                                                                                  >  2:463087/1‑251 (MQ=255)                                                                                                                          tGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                           <  1:482581/250‑1 (MQ=255)                                                                                                                           gATATAGTTTCTTGCCGTCTCTCCACATCTTGTTTTTCTACAGTCTAAAAAAAGGCGACTCCTCAGTCTCCTTAAAAATCCGTTTGCCCGCGCCGCCTGCTGTGTCCCTCCAATTCATTCGA‑GGC‑CCCC‑TCGCGCCAGCGCCAACTAGGTGAGTAGCTCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctcttt                                                                                                           <  2:350526/249‑3 (MQ=255)                                                                                                                              atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGACGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTACTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCgtctg                                                                                                      >  2:289857/1‑247 (MQ=255)                                                                                                                              atGGTTACATGCCTTCGCTTCTCAGCTTCCCTTTCTCCAGACAAAAAAAAGGCGATTCATCAGTCGCTTTAAAAACCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTAAAATGTATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCGTTATGGGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCATGCGGTCATCCCTACTCTCATGTCGGTCTCTTCGGct                                                                                                       <  2:381285/250‑2 (MQ=255)                                                                                                                                        gccctCGCTCCTCATGTTACTTTTCTACAAACAAAAAAAAGGCGACACATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTTCAATGGATTCGC‑T‑T‑CCCCCTCGGGCCAGCGTCAACAAGGTGCTTAGTGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCaat                                                                                            <  2:54483/247‑2 (MQ=255)                                                                                                                                              gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTATCAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCCACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATTCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCTTGGTCTTCAACGTCTTc                                                                                      >  2:156377/1‑251 (MQ=255)                                                                                                                                                 cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                     >  1:366602/1‑249 (MQ=255)                                                                                                                                                    cATCTGACTTCTCTAACGACAAAAAAAAGGCGACACAGCAGTCGCCTTAAAAATCACTTTGCCGGCGCCGCCTTCTGCGTCGCTACAATGGATTCGACT‑C‑CCCC‑GCTGGCCAGCGCCAACCAGTTGTGTAGCTCTTCATTGTTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGTCTGCGTCTTCAACGTATTCCAGagc                                                                                <  2:20196/251‑3 (MQ=255)                                                                                                                                                         tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaa                                                                            >  2:252826/1‑250 (MQ=255)                                                                                                                                                          aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaa                                                                           <  1:137367/250‑1 (MQ=255)                                                                                                                                                             tttCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaaattgc                                                                       >  2:161382/1‑247 (MQ=255)                                                                                                                                                                        aaataaaaGGCGACTAATCAGTCGCCTTTCACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TTCACA‑TCAGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTCCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCATGTCGGTCTCGGCGGCGGCGTCTGCCACGTATCCCAGATCGCaatttgcttcaacgttt                                                             >  2:28341/5‑235 (MQ=255)                                                                                                                                                                           aaaaaGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTAATCACTACGTTCATGTCGGTCTCTGCGCCTGAGTCTTCAACGTATTCCAGCTCGctaatttcttcaccgttttct                                                          >  2:256339/1‑230 (MQ=255)                                                                                                                                                                              atggCGATTCATCAGTAACCTTAAAGATCAGTTTGACAGCCACGAATTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCTGAGTCTTCAACGTATTTCAGCTCGCAAATTGCTTCCCCGTTTACAAtac                                                      >  2:273046/3‑249 (MQ=255)                                                                                                                                                                                 gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                    >  1:350878/1‑250 (MQ=255)                                                                                                                                                                                 gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTAAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCCTCACTACGTTCATTTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                    >  1:245359/1‑250 (MQ=255)                                                                                                                                                                                       catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGa                                               <  2:377798/249‑1 (MQ=255)                                                                                                                                                                                          caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                          >  1:355915/1‑251 (MQ=255)                                                                                                                                                                                                 cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTTCGGCTAc                                   >  1:413564/1‑251 (MQ=255)                                                                                                                                                                                                     aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATc                               >  1:284774/1‑251 (MQ=255)                                                                                                                                                                                                          cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCCCCGTTTACAATACCCACAGAAACTGCGGCTACCATCCCttt                          >  2:41603/1‑251 (MQ=255)                                                                                                                                                                                                                gccacccccgccTTCTGCGTCTCTACAATGGATTCGAGT‑C‑CCCC‑TCGGGCCAGCGCCAACAAGGTGGGGAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGGCCGGCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTGCACCGTTTACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                     <  1:273046/243‑1 (MQ=255)                                                                                                                                                                                                                  cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                            >  2:344277/1‑97 (MQ=255)                                                                                                                                                                                                                  cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                            <  1:344277/97‑1 (MQ=255)                                                                                                                                                                                                                     cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCACGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTgg               >  1:102560/1‑251 (MQ=255)                                                                                                                                                                                                                      gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt              >  1:451433/1‑251 (MQ=255)                                                                                                                                                                                                                      gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt              <  1:3097/251‑1 (MQ=255)                                                                                                                                                                                                                      accgccTTCTGCGTCGATACAATGGATTCGA‑T‑T‑CCCCCTCGGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATTATGCGCCCATCTTCGGTCATCACTACGTTCATGTCGATCTCTTCGGCTGCGCCTTCAACGTATTCCCGCTCGCACATTGCTTCACCTTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCggttttgt              >  2:63/2‑245 (MQ=255)                                                                                                                                                                                                                         gccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGtttt           >  1:58825/1‑251 (MQ=255)                                                                                                                                                                                                                             tcggCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGGGCCAGCGTCAACAAGGTGAGTAGTTCGTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAg        <  2:333823/247‑1 (MQ=255)                                                                                                                                                                                                                              cTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTAAGCt      >  2:484868/1‑251 (MQ=255)                                                                                                                                                                                                                                gCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATAGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTg    >  2:54780/1‑251 (MQ=255)                                                                                                                                                                                                                                 cGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACTGCGGCTACCATCCCCTTCATCGGATTGGTGTTAAGCTTGc   >  2:441607/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTAAGCtttcc  >  2:336014/1‑248 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCC  >  NC_000913/3815634‑3816108

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