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

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

ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815634‑3816132                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          aTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTcc‑ccc                                                                                                                                                                                                                                                           <  1:327404/251‑2 (MQ=255)       cccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGc                                                                                                                                                                                                                                                     >  2:299720/1‑251 (MQ=255)           gTATTAAACAGCCCGGCGTTGCAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTACGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGTTTA‑GATTCC‑CCCTCGGGcc                                                                                                                                                                                                                                                    >  2:332980/1‑248 (MQ=255)           gTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGCCGCCTTCTGCGTCGCTACAACGGATTC‑GATTCC‑CCCTCGGGCCagg                                                                                                                                                                                                                                                 >  1:314169/1‑250 (MQ=255)             aTTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTTCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGc                                                                                                                                                                                                                                               <  2:314169/251‑1 (MQ=255)               tAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGcc                                                                                                                                                                                                                                              >  2:85671/1‑250 (MQ=255)                  aCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCcaaca                                                                                                                                                                                                                                          <  1:233538/251‑1 (MQ=255)                   cAGCCCGGCGGTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCATAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCcaaca                                                                                                                                                                                                                                          >  2:250192/1‑250 (MQ=255)                    aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCcaaca                                                                                                                                                                                                                                          >  2:239229/1‑249 (MQ=255)                      atcggacttcatgaactAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAACCTTTCACACCTGCTTGCTAAGCGCAAATTAAATCAACTAGCGCTGATATGGTTTCATGACTTCGATCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCAATGGATTC‑GATTCCACC‑TCGGGCCAGAGCcaaca                                                                                                                                                                                                                                          >  2:61903/17‑247 (MQ=255)                       ccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                                     >  1:134651/1‑251 (MQ=255)                        cGGCGTTGAAGAATTAGGGGCTTTTGCGCCCGGGTTTCAGCGTAACCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCCTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                                     >  2:312324/1‑250 (MQ=255)                        cGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCACCAAGATg                                                                                                                                                                                                                                     >  2:410557/1‑250 (MQ=255)                          gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                  <  2:273187/251‑1 (MQ=255)                           cGTTGAAGAACTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGAGCAAATTCAATCAACTGGCGCTGATATGGTTTCATGCCTTAGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCCACAAGATGAga                                                                                                                                                                                                                                  >  2:443419/1‑249 (MQ=255)                           cGTTGAAGAAATGTGGGCTTTTGAGCCCGGATTTCAGCGTAAAATCGCCAAACATGCACACCTGCTTTCTAAGCGCAAATACAATAAACTGGCGCTGATATGGTTTCATGCCATCGCTCCACATCTTACTTATCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCGTCGCTACAATTGAGACTAATCCC‑CC‑TCGGGCCAGAGCCATCAAGaagatta                                                                                                                                                                                                                                 <  2:449680/251‑7 (MQ=255)                           cGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                  >  2:244954/1‑250 (MQ=255)                                 agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                                                                                                                                                           >  1:278103/1‑251 (MQ=255)                                  gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCACTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                                                                                                                                                                                          >  1:37730/1‑251 (MQ=255)                                   aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTACCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                                                                                                                                                                                          >  2:447085/1‑250 (MQ=255)                                              tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                <  1:94184/249‑1 (MQ=255)                                                 gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                                           <  1:244528/251‑1 (MQ=255)                                                  cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                                                                                                                                          <  1:239229/251‑1 (MQ=255)                                                    cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGcc                                                                                                                                                                                                        >  2:367115/1‑251 (MQ=255)                                                      cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                                      <  1:8439/251‑1 (MQ=255)                                                        gATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCGTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                    >  1:171872/1‑251 (MQ=255)                                                            tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 <  1:151408/250‑1 (MQ=255)                                                                   aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCGTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCt                                                                                                                                                                                          <  2:171872/250‑1 (MQ=255)                                                                            aaaCTATAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACACATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                                <  1:62407/251‑1 (MQ=255)                                                                             aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                                                                                                                                                                               <  1:212468/251‑1 (MQ=255)                                                                                 ttAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCTTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                                                                                                             <  2:67872/249‑1 (MQ=255)                                                                                  tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCATGGGAGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                                                                                                          >  1:270892/1‑251 (MQ=255)                                                                                    acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGTTATGGTTTCACTCCTTCGCTCCTCAACGTACGTTTCGACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGGGCCGCCGTCTGCGTCGCTACAATGTATTC‑GATTCC‑CCCGCTGGCCATAGCCAACAAGAGGAGTATCTCTTGATGGGTGGACGGCCCGCCTTTGGCCGTCCTCTGCACTTCAATGATGCGcc                                                                                                                                                                        <  2:218097/251‑1 (MQ=255)                                                                                             tgctTAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                                 <  2:332504/249‑1 (MQ=255)                                                                                                   gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                                                                          >  2:429279/1‑250 (MQ=255)                                                                                                        aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCTACAGACAAAAAACAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCATGCACTTCAATGATGCGCCAGTCTTCGGTCATCACTACg                                                                                                                                                    >  1:76049/1‑251 (MQ=255)                                                                                                           tCAAGAAACTGGCGCTGATATGGTTTCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                                                                  <  2:458341/250‑1 (MQ=255)                                                                                                            cAATAAACTGGCGCTGATATGGTTTAATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTcc                                                                                                                                                >  2:236442/1‑250 (MQ=255)                                                                                                                 aaCTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATgacg                                                                                                                                           <  2:120064/251‑4 (MQ=255)                                                                                                                       cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                                                                     >  1:131307/1‑251 (MQ=255)                                                                                                                         cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                   >  1:179568/1‑251 (MQ=255)                                                                                                                              atGGTTACAGGCCTTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                                                                              <  1:244954/251‑1 (MQ=255)                                                                                                                                 gTTTCATTCCTTCGCGCCCCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCGTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                                                                                           <  2:433568/251‑1 (MQ=255)                                                                                                                                   ttCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCTGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                                                                          <  2:40778/250‑1 (MQ=255)                                                                                                                                           ttCGCTCCTCATCTTACTTTTCTACAGACACAAAAAAGGCGACTCATCAGTCGCATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                                  >  1:347241/1‑250 (MQ=255)                                                                                                                                           ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                                  >  1:48576/1‑250 (MQ=255)                                                                                                                                            tctctcctcATCTTACTTATCTACAGACAAAAAAAAGGCGACTCAACAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATGC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCTTCAACGt                                                                                                                  <  2:96804/246‑1 (MQ=255)                                                                                                                                                 cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATGC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg                                                                                                           <  1:301341/251‑1 (MQ=255)                                                                                                                                                                 tACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                            <  2:288679/122‑1 (MQ=255)                                                                                                                                                                 tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCCACC‑TCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                            >  1:288679/1‑122 (MQ=255)                                                                                                                                                                     gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGt                                                                                       >  2:154547/1‑251 (MQ=255)                                                                                                                                                                       cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCACCGCCGACTTCTGCGTCGCTCCAATGGATTC‑GATTCC‑CCCTCGGGCCAAAGCCAACAAGACGAGCAGCTCTTCATGGGTGCACGGCACGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTATTCGGTCATGCCTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                                                                     >  1:420860/1‑251 (MQ=255)                                                                                                                                                                          aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcaca                                                                                  <  1:90044/251‑1 (MQ=255)                                                                                                                                                                            aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATTCGCCCGTCTTCGGTCATCACTACGTTCAGGTCTGTCTCTGCGGCAGAGTCGTCAACGTATTCCAGATCGCAATCCGCTTCGCCGTTCACAAt                                                                                <  2:270038/251‑1 (MQ=255)                                                                                                                                                                                ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                                                                            >  1:113565/1‑251 (MQ=255)                                                                                                                                                                                             tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                               >  1:28201/1‑251 (MQ=255)                                                                                                                                                                                                     aaaaTCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATc                                                       >  1:131122/1‑251 (MQ=255)                                                                                                                                                                                                       aaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGCCAGAAACTGAGGCTACCATccc                                                     >  2:56041/1‑251 (MQ=255)                                                                                                                                                                                                          cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACCGAAACTGCGGCTACCATCCCttt                                                  >  1:402277/1‑251 (MQ=255)                                                                                                                                                                                                             tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCCGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAt                                               >  2:99472/1‑251 (MQ=255)                                                                                                                                                                                                                   aGCGCCGCCTTCGTCGGCGCTACAATGGATTC‑GATTCC‑CCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATTGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGAtt                                         <  1:7447/251‑1 (MQ=255)                                                                                                                                                                                                                     cgccgcCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                       >  1:132303/1‑251 (MQ=255)                                                                                                                                                                                                                         gccTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                   >  1:421560/1‑251 (MQ=255)                                                                                                                                                                                                                                gCGTCGCTACAATGGATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTg                            <  1:367115/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCTACAATGGATTC‑GATTCC‑CCCTCGGGCAAGAGCCCACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTACCGtt                       >  2:385200/1‑251 (MQ=255)                                                                                                                                                                                                                                             ggATTC‑GATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacc                >  1:121756/1‑250 (MQ=255)                                                                                                                                                                                                                                               attc‑gattcC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGc             >  1:321939/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gattccc‑cc‑TCGGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCt         >  2:331144/6‑251 (MQ=255)                                                                                                                                                                                                                                                    attccc‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTg        >  1:115103/5‑251 (MQ=255)                                                                                                                                                                                                                                                      tccc‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGCTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  2:221248/3‑250 (MQ=255)                                                                                                                                                                                                                                                            c‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGAAgcgc  >  1:283776/1‑250 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          ATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC‑GATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815634‑3816132

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