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.8%	150.5 / 2.2	48	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (30/15);  total (31/17)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCT  >  NC_000913/3815638‑3816125                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    cgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTc‑cccctcgg                                                                                                                                                                                                                                                 >  1:313934/1‑246 (MQ=255)     cGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCa                                                                                                                                                                                                                                             <  1:148831/251‑1 (MQ=255)       gTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCaga                                                                                                                                                                                                                                           <  2:314554/251‑1 (MQ=255)          ttAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGcc                                                                                                                                                                                                                                        >  1:267027/1‑251 (MQ=255)                 gCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCc‑aaca                                                                                                                                                                                                                                   >  1:136007/1‑248 (MQ=255)                    cGGCGATGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATg                                                                                                                                                                                                                              >  2:232337/1‑250 (MQ=255)                           gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGctc                                                                                                                                                                                                                      >  1:239858/1‑251 (MQ=255)                                aaTAGGGGTTTTTGCGCCCGGCTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCACCTTCCTTTTCTACAGACAAAAAAAAGCAGACCCATCAGTCGCCTTAAAAATCATTTTCCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCA‑AGCCCAACAAggcgatttgcctttctt                                                                                                                                                                                                                 >  2:464018/1‑235 (MQ=255)                                             gcgcCCGGATTTCAGCGTAAACTCTCCAAATGTTAACACCTGCTTGCTAAGCGCAAAATCAATAATCTGGCGCTGATATGGTGTCATGCCTTCGCTCCACATCTGACGTGGCTACAGACAAAAAAAAGGCGACTCATCAGTCGACGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCAATGGATGCGA‑TTC‑CCCCTAGGGCCAGAGCC‑GACAAGATGAGTCGCTCTTCATGGGTGAACGGCt                                                                                                                                                                                                     <  2:386636/250‑1 (MQ=255)                                              cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                                                                                                                                   >  2:34882/1‑251 (MQ=255)                                                 ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                                                                                                                                                                                                >  1:228377/1‑251 (MQ=255)                                                  cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAAGAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                               <  2:425610/251‑1 (MQ=255)                                                          aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                                                                                                                                                       >  1:386207/1‑251 (MQ=255)                                                                    cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                             <  2:136007/251‑1 (MQ=255)                                                                     ggcAAACTTGAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGTCACTAATATAGTTTCATTCCTTCGCTCCTCATCATACTTTTGTACAGACAAAAAAAAAGCCACTCATCAGTCGGCGTAAAAATCAGGTTGCCAACGCCGCCATCTGCGTCGCTTCAATGGATTCGA‑TTC‑CCCCGCGGGCCAGAGCC‑AACAAGATGATGAGCTCTGGATGTGTGAACGGCTCGCATTCTGACAACCCCCGCACTTc                                                                                                                                                                            <  2:22279/249‑1 (MQ=255)                                                                            tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                     >  2:154033/1‑251 (MQ=255)                                                                                            tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcc                                                                                                                                                     >  2:420742/1‑250 (MQ=255)                                                                                                         aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                        >  2:123417/1‑251 (MQ=255)                                                                                                                      tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                           >  1:240226/1‑251 (MQ=255)                                                                                                                               ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                  >  1:145148/1‑251 (MQ=255)                                                                                                                                  aTGCCTTCGCTCTTAATCTTACTTTTCTACATACCAAAAAAAGGCGTCTCATCAGTCACCTTAAAAATCAGCTTGCCAGCGCGTCATTCTGCGTCGCTACAATGTATTAGAATGC‑CCC‑TCGGGACAGGGCC‑ATCAAGATGAGTAGCTCTTCATGGGTGAACGGCGCGCCTTCTGTCGTCCCCGGCGCTTCAATGATGCGCCCGTCTTCGGTCATCACTGCGTTCACGGCAGTCTCTGCGACAGAGTCTTCa                                                                                                               <  2:265296/251‑1 (MQ=255)                                                                                                                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCg                                                                                                                                                               >  1:456573/1‑203 (MQ=255)                                                                                                                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCg                                                                                                                                                               <  2:456573/203‑1 (MQ=255)                                                                                                                                         cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                                                                        <  2:161748/251‑1 (MQ=255)                                                                                                                                             cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                      >  2:363680/1‑249 (MQ=255)                                                                                                                                                 aTCTTACTTTTCTATAGAAAAAAAATAGGCGACTCATCAGTCGCCTTCAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCACTGCATTCGA‑TTC‑CCCCTCGGCCCAGAGCC‑AACAAGATGAGTAGCTCTTCACGCGTGAACGGCTCGCCTTCTCACGTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTATTcc                                                                                                      >  2:301206/1‑245 (MQ=255)                                                                                                                                                      aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc                                                                                           <  2:340830/251‑1 (MQ=255)                                                                                                                                                       cTTTTCTACAGACAAACAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCCCAAAcc                                                                                          >  2:402885/1‑251 (MQ=255)                                                                                                                                                                        aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAt                                                                         >  2:203943/1‑251 (MQ=255)                                                                                                                                                                         aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtt                                                                        >  1:140986/1‑251 (MQ=255)                                                                                                                                                                                  tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCACCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGa                                                                >  1:391390/1‑250 (MQ=255)                                                                                                                                                                                   catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaaa                                                              <  2:108219/251‑1 (MQ=255)                                                                                                                                                                                    atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGTTCACAATTCCGACAGAAAc                                                             >  1:345325/1‑251 (MQ=255)                                                                                                                                                                                        gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg                                                         <  1:94337/251‑1 (MQ=255)                                                                                                                                                                                           gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                      >  1:149442/1‑251 (MQ=255)                                                                                                                                                                                           gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                      >  1:29483/1‑251 (MQ=255)                                                                                                                                                                                           gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATGCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                      <  2:436718/251‑1 (MQ=255)                                                                                                                                                                                                            gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                  <  2:157673/174‑1 (MQ=255)                                                                                                                                                                                                            gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                  >  1:157673/1‑174 (MQ=255)                                                                                                                                                                                                                   ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                      >  2:284677/1‑243 (MQ=255)                                                                                                                                                                                                                   ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                      <  1:284677/243‑1 (MQ=255)                                                                                                                                                                                                                      ccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACATTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGCTTTc                           >  2:245273/1‑251 (MQ=255)                                                                                                                                                                                                                       cTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                           >  2:397941/1‑250 (MQ=255)                                                                                                                                                                                                                          cTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                       >  1:401140/1‑251 (MQ=255)                                                                                                                                                                                                                               tCGCTACAATGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCg                  >  2:88451/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaTGGATTCGA‑TTC‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcca           <  1:174702/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggATTCGAGTCC‑CCC‑TCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAGCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTGCAGCGGATTGGTTTTCAGCGGGCCGTTTGccacc         <  2:297538/250‑1 (MQ=255)                                                                                                                                                                                                                                           attcga‑ttc‑CCCCTCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGATTGGTTTTCAGCTTGCCGTTTTccacccgc      >  2:334444/1‑248 (MQ=255)                                                                                                                                                                                                                                              cga‑ttcaCCC‑TCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTc   >  1:378078/1‑251 (MQ=255)                                                                                                                                                                                                                                               gattcc‑ccc‑TCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCt  >  1:403741/6‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TCGGGCCAGAGCC‑AACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCT  >  NC_000913/3815638‑3816125

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