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.1%	119.1 / 11.0	47	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base C (24/17);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTC  >  NC_000913/3815646‑3816114                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         aTTAAACAGCCCGGCGGTGAAGAAATAGGGGCTTTTGCGCCCGGATGTCAGCGTAAACTCGCCAAACTTTAACGCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTGAAAAATCAGGTGGGCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGc                                                                                                                                                                                                                             <  1:139834/251‑1 (MQ=255)                  tGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGAGAAAACTCGCCCAACTTTAACCCCTGCTTGCTAAGAGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGCGCCAACAAGATGAGTAGct                                                                                                                                                                                                            >  1:195238/1‑251 (MQ=255)                       aaaTAGGGGCTTTTGCGCCCGGCTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                                                                                                                                                                        >  2:125292/1‑250 (MQ=255)                           aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACACTGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACACGATGAGTAGCTCTTCATggg                                                                                                                                                                                                   >  2:165185/1‑251 (MQ=255)                              ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                                                                                                                                 >  2:88560/1‑250 (MQ=255)                                  tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCTTCGCTACAATGTATTCTA‑TGCCC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGTGTGAACg                                                                                                                                                                                             <  2:308987/250‑1 (MQ=255)                                     gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                                                                                                                          <  1:146264/250‑1 (MQ=255)                                     gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTGCATGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑TGCGACTCATCAGTCGCCTTAAAAATCAGTGTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                         <  1:269733/251‑1 (MQ=255)                                         ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                                                                                                                                                                                     >  1:187739/1‑251 (MQ=255)                                           ggATTTCAGAGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTATCTCTTCATGGGTGAACGGCTCGCCTTc                                                                                                                                                                                   >  2:288858/1‑251 (MQ=255)                                             aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCCACAGACAAAAAAAA‑GGCGACTCATCAGTCACCTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGACCAGAGCCAACAAGATGAGTAGCTCTTCATGGTCGAACGGCCCGCCttatg                                                                                                                                                                                 >  2:232715/1‑248 (MQ=255)                                                 cAGCGTAACCTCGCCAAACATTACCCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGCTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑C‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCg                                                                                                                                                                              >  2:315130/1‑251 (MQ=255)                                                            cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                  >  2:269671/1‑251 (MQ=255)                                                               cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCa                                                                                                                                                                >  2:175462/1‑250 (MQ=255)                                                                                 tgatAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                               <  2:234907/246‑1 (MQ=255)                                                                                     aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCTACAAGGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAGCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                                                           <  1:315130/250‑1 (MQ=255)                                                                                       gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                                                        >  1:247622/1‑250 (MQ=255)                                                                                           aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                                                                                                                                                                                     >  1:264431/1‑201 (MQ=255)                                                                                           aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCAGCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACTGCGCGCCt                                                                                                                                                                                     <  2:264431/201‑1 (MQ=255)                                                                                                 aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                             <  1:170377/251‑1 (MQ=255)                                                                                                 aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                             <  1:175462/251‑1 (MQ=255)                                                                                                                ttatGTTTTCATTCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTGAAAAATCACTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCATTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                                                               <  2:281709/249‑1 (MQ=255)                                                                                                                   tGGTTTCATGCCATCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGA‑T‑TCC‑CAC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGAACTTCAATGATGCGCCCGTCTTCGGTCATCAATACGTTCATGTCGATCTCTGAGGc                                                                                                              >  2:167841/1‑248 (MQ=255)                                                                                                                            gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATTGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                                                   <  1:125292/250‑1 (MQ=255)                                                                                                                                tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                                                               <  1:248292/250‑1 (MQ=255)                                                                                                                                   ctcctcATCTAACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTACAATGGATTCGA‑T‑TCCAC‑C‑TC‑GGGCCAGAGCCAACAAGAAGAGTAGCTATTCATGGGGGAACGGCTCGCCTTATGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCagcgtattacacgacttcc                                                                                           >  2:94370/1‑234 (MQ=255)                                                                                                                                          tCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                         >  1:102281/1‑150 (MQ=255)                                                                                                                                          tCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                         <  2:102281/150‑1 (MQ=255)                                                                                                                                                          aCAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTTCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                                                    >  2:49324/1‑251 (MQ=255)                                                                                                                                                           cgaaaaaaa‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                                                   >  2:168889/3‑251 (MQ=255)                                                                                                                                                               aaaaa‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcaca                                                                >  1:26821/1‑250 (MQ=255)                                                                                                                                                                aaaa‑GGCGACTCATCAGTCGAATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCACTGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcacact                                                              >  2:306585/1‑249 (MQ=255)                                                                                                                                                                             atcaGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCGTCACTACACTGGAATCGC‑T‑TCC‑C‑CATC‑GGGCCAGCGCCCACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTATTCCAGATc                                                                                      >  2:157928/1‑215 (MQ=255)                                                                                                                                                                             agcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCAGCACTACGTACATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                                                      <  1:157928/213‑1 (MQ=255)                                                                                                                                                                              tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                 >  1:217128/1‑251 (MQ=255)                                                                                                                                                                                 gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg                                              >  1:247693/1‑251 (MQ=255)                                                                                                                                                                                 gTCGCCTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATGCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCTTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGGCGTTCGCAATTCCGACATAAACTGCg                                              <  2:210319/251‑1 (MQ=255)                                                                                                                                                                                     ccTTAAAAATCAGTTTGCCAGCGACGCCTTCTGCGCCGCTAAAATGGATTCGACT‑CCC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAAGGCTCGCCTTCTGCCGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGAGTCAGGGTCTTCAACGTCTTCCAGATCCCAAACCGCTTCACCGGTCACAATTCCGACAGCCACTGCgcctt                                          >  2:228877/1‑247 (MQ=255)                                                                                                                                                                                      cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc                                         >  2:160650/1‑251 (MQ=255)                                                                                                                                                                                                     gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCgg                          <  1:160650/251‑1 (MQ=255)                                                                                                                                                                                                        aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACCGAAACTGCGGCTACCATCCCTTTCATCGGAtt                       >  2:112664/1‑251 (MQ=255)                                                                                                                                                                                                             cgccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                  <  1:49324/251‑1 (MQ=255)                                                                                                                                                                                                               ccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                <  2:247693/251‑1 (MQ=255)                                                                                                                                                                                                                cTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                >  1:155783/1‑250 (MQ=255)                                                                                                                                                                                                                       gTCGCTAAAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCAAGAGACAACAAGACGAGTAGCTCTTAATGGGTGAACGGCTCGCCTTCTGACGTCCCCAGCACTTCAATGTAGCGCCCGTCTTCGGTCCTCACTACGTTCATGCAGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGAGGCTACCATCCCTTTCATCGGATTGGTTTTCCGCTTGcc        >  1:119634/1‑251 (MQ=255)                                                                                                                                                                                                                          gCTACAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtt     >  1:59008/1‑251 (MQ=255)                                                                                                                                                                                                                             aCAATGGATTCGA‑T‑TCC‑C‑CCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTAAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGAGTGGTTTTCAGCTTGCCGtttgc  >  2:163842/1‑249 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTC  >  NC_000913/3815646‑3816114

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