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.5%	129.2 / 3.9	46	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (32/11);  total (35/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

CCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815640‑3816132                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           cccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGc                                                                                                                                                                                                                                                     >  1:53243/1‑251 (MQ=255)    gggaTTAAACAGCCCGGCGTTGAAGAAATAGGTGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCag                                                                                                                                                                                                                                                  <  2:181356/248‑1 (MQ=255)               gcacGGCGTTGAAGAACTAGGAGCTTTTGCGACCGGATTTCAACGTAACCGCGCCAAACTTTAACACCTGCTTGCTAAGCGCCAATTAAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACACAAAAAAGGCGACTCATCAGTCGCCTTAATAATCAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACACGGGATTCGATTCAC‑CC‑TC‑GTGCCAGAGCcaacaac                                                                                                                                                                                                                                        >  2:118861/4‑249 (MQ=255)                   ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                                     >  1:35972/1‑249 (MQ=255)                                tAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                                                                                                                                                      >  1:59094/1‑251 (MQ=255)                                           gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTATCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                                           <  2:59094/251‑1 (MQ=255)                                            cgaccGGATTTCAGCGTACACTCGCCAAACTTTAACACCTGCTTGCTAAGAGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAACTGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCg                                                                                                                                                                                                          >  2:225575/4‑251 (MQ=255)                                                                cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCAc                                                                                                                                                                                      >  1:283239/1‑251 (MQ=255)                                                                       aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                                                                                                                                                                               >  1:193172/1‑251 (MQ=255)                                                                          tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                            >  1:235445/1‑251 (MQ=255)                                                                             aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                                                                                                          <  1:263921/250‑1 (MQ=255)                                                                                 ccTGCTTGCTAAGCGCAAATTCAGTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                                                                                                                     >  2:124728/1‑251 (MQ=255)                                                                                           aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTAAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                           >  2:325581/1‑251 (MQ=255)                                                                                                        aTAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGTATCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCCACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTTCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                                                                                                              >  1:118675/1‑251 (MQ=255)                                                                                                         tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                                                                             >  1:91488/1‑251 (MQ=255)                                                                                                          aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTGCATGt                                                                                                                                             <  2:20226/250‑1 (MQ=255)                                                                                                           aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACCGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                           >  1:198586/1‑251 (MQ=255)                                                                                                             cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGACAAAAAACAGGCGACTCATCATTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CACC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgtt                                                                                                                                         >  2:173399/1‑249 (MQ=255)                                                                                                                     gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                                                                                 >  1:267891/1‑251 (MQ=255)                                                                                                                        atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                                                                               <  1:325581/250‑1 (MQ=255)                                                                                                                                aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                                                                       >  2:276266/1‑250 (MQ=255)                                                                                                                                    cTTCGCTCCTCATCTTACTTTTCTTCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                                  >  2:80345/1‑251 (MQ=255)                                                                                                                                         ctcctcATCTTAATTTTCTACAGACAAAAAAAAGGCGCCTAATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGTCTCTTAGACAGAGTCTTCAACGTATTcc                                                                                                             >  2:310930/1‑251 (MQ=255)                                                                                                                                                            aCAGACAAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGCTTCGATTC‑C‑CCCTC‑GGGACAGAGCCAACAAGATGAGTAGATCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCCCCCGTCTTAGGTCATCACTACGTTCATGTCGGTCTCTGAGGCAGAGTCTTCAACGTATTCTAGATCCCAAACCGCTTCGc                                                                                          >  1:218010/1‑251 (MQ=255)                                                                                                                                                              aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                                                                        >  1:108340/1‑251 (MQ=255)                                                                                                                                                                aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCGC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                                                                      <  1:247371/251‑1 (MQ=255)                                                                                                                                                                aCAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGtt                                                                                      <  2:283239/251‑1 (MQ=255)                                                                                                                                                                  aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CC‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcc                                                                                    >  2:89823/1‑250 (MQ=255)                                                                                                                                                                    aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                                   >  1:20421/1‑250 (MQ=255)                                                                                                                                                                     aaaaaGGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACCAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGttaccc                                                                                  >  2:152997/1‑246 (MQ=255)                                                                                                                                                                        aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                                                                              >  1:267609/1‑251 (MQ=255)                                                                                                                                                                          ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                                                                            >  1:22207/1‑251 (MQ=255)                                                                                                                                                                              aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  1:278367/1‑249 (MQ=255)                                                                                                                                                                                         gCCTTAAAAATCAGTGTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCATATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                             <  1:128373/251‑1 (MQ=255)                                                                                                                                                                                           cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAAATGCGGCTAc                                                           >  2:229883/1‑251 (MQ=255)                                                                                                                                                                                                         tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                             >  1:245280/1‑251 (MQ=255)                                                                                                                                                                                                         tGCCAGCGCCGCCTTCCGCGTCGCTACAATGGATGCGATTC‑C‑CCCGC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcccctgcatc                                              <  2:53865/250‑8 (MQ=255)                                                                                                                                                                                                            cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGAt                                          >  1:162571/1‑251 (MQ=255)                                                                                                                                                                                                               cgccgcCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                       >  1:166161/1‑251 (MQ=255)                                                                                                                                                                                                                  cgccTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                                    <  2:59037/251‑1 (MQ=255)                                                                                                                                                                                                                  cgccCTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                                    >  2:266160/1‑251 (MQ=255)                                                                                                                                                                                                                   gccTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                   <  1:195993/251‑1 (MQ=255)                                                                                                                                                                                                                        cTGCGTCGCTACAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                              >  1:221855/1‑251 (MQ=255)                                                                                                                                                                                                                                   cAATGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcc                   >  1:13427/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaTGGATTCGATTC‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATAACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcc                   >  1:182452/1‑250 (MQ=255)                                                                                                                                                                                                                                          ttcgattc‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCt            >  2:153568/1‑251 (MQ=255)                                                                                                                                                                                                                                                   c‑cc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcg   <  2:236077/251‑1 (MQ=255)                                                                                                                                                                                                                                                     cc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                                   >  2:282162/1‑186 (MQ=255)                                                                                                                                                                                                                                                     cc‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                                   <  1:282162/186‑1 (MQ=255)                                                                                                                                                                                                                                                      c‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGATTCTGTAgcgc  >  2:149937/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           CCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑C‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815640‑3816132

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