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	97.7%	136.8 / ‑3.3	44	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (27/16);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTA  >  NC_000913/3815636‑3816128                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCGTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTc‑cccct                                                                                                                                                                                                                                                       <  2:298740/250‑3 (MQ=255)       cGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCa                                                                                                                                                                                                                                                <  1:139101/251‑1 (MQ=255)                              aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                         >  1:376060/1‑251 (MQ=255)                                                  cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGcc                                                                                                                                                                                                     >  2:382294/1‑251 (MQ=255)                                                         ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                               <  1:339056/250‑1 (MQ=255)                                                           catcGTAAACCCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCATGCCTGCGCTCCGCATCTTTCGTTTCTACAGACAAAAAAAAGGCGGCTCTTCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCGAACAAGATGAGTATCTCTTCATGGGTTAACGGCTCGCCTACTGCCGt                                                                                                                                                                                            <  2:535800/248‑1 (MQ=255)                                                                      cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                                 >  1:103701/1‑251 (MQ=255)                                                                       gCCAAACTTTACCTCCGGCGTGCTTAGCGCAAATTCAATAATCGGGCGCTTATATGGTGTCATGTCTTCTCTCCATATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCATGTTGCCAACGCCGCCCTCTGCGTCGCTACAATGGACTCGATTC‑CCCCTCGGGCCAGTGCCAACAAGATGAGTTGCTCTTCATGTGTAAACGTCTCGCCTTCTCCCTTccccgtcccttc                                                                                                                                                                                <  2:207514/251‑9 (MQ=255)                                                                             cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGACAGCGCCAACAAGTTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCACTGCACTTCAaggat                                                                                                                                                                          >  2:26158/1‑251 (MQ=255)                                                                                        gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                                                                 <  1:56016/249‑1 (MQ=255)                                                                                            gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                           >  1:376676/1‑251 (MQ=255)                                                                                                     aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCa                                                                                                                                                                                                                      >  2:452331/1‑183 (MQ=255)                                                                                                     aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                                                                                                                                                                                      <  1:452331/183‑1 (MQ=255)                                                                                                           aaTAGACTGGCGCTGATATGGTTTCATGCCTTCGCTACTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACCAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGAGCCCGTCTTCGGTCATCACTCCGTTctt                                                                                                                                            >  2:433664/1‑249 (MQ=255)                                                                                                            aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGCCCTCACTACGTTcctg                                                                                                                                           >  2:425516/1‑248 (MQ=255)                                                                                                               aaCTGGCGCTGCTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGACTCATCAGCCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGATTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                        <  2:2323/251‑1 (MQ=255)                                                                                                                     cgcTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                                                                  >  2:521320/1‑251 (MQ=255)                                                                                                                          atatGGTTTCATGCCTTCGCTCCACATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGc                                                                                                                             <  2:6081/251‑1 (MQ=255)                                                                                                                                          tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                                                                                                             >  1:70402/1‑251 (MQ=255)                                                                                                                                                    tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGc                                                                                                   >  1:181988/1‑251 (MQ=255)                                                                                                                                                       tacgtttCTACAGACAAAAAAAAGGCGACGCATCAGCCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg                                                                                                    <  2:70402/243‑1 (MQ=255)                                                                                                                                                                        aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTATTCACCGTATTCCAGATCGCAAACCGCTTCGCC‑GTTcac                                                                               >  2:146324/1‑250 (MQ=255)                                                                                                                                                                             aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGATACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAAGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTcc                                                                         >  2:309191/1‑251 (MQ=255)                                                                                                                                                                                 gacccatcatTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGAGGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACa                                                                     <  1:244227/247‑1 (MQ=255)                                                                                                                                                                                          gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCg                                                            >  1:283932/1‑251 (MQ=255)                                                                                                                                                                                            cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCgg                                                           <  1:86620/250‑1 (MQ=255)                                                                                                                                                                                            cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGc                                                          >  1:477118/1‑251 (MQ=255)                                                                                                                                                                                                    aaagccgttcgCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCTGCCCAGAGCCAACAAGACTAGTAGCTTTTCATGGTGGAACGGCTCGCCTTCTGCCTCCCCCTGCACTTCAATGATGCGCCCGTCGTCTGTCATCACTTCGTTCATGTCGGTCTCTGCGGCAGAGTCGTCAACTAATTCCATATCGCAACCCGCTGCGCC‑GTTCACAATTCCGTCAGAAACTGCGGCTCCCATcc                                                  <  2:488563/241‑1 (MQ=255)                                                                                                                                                                                                                 cgcgcCGCCTTCTGCGTCGCTAAAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAACTCTTCATGGGTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTAATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCATAACGCTTCGCC‑GTTCACAATTCCGACAGTAACTGCGGCTACCATCACTTTCCTCggctt                                     >  2:117879/2‑248 (MQ=255)                                                                                                                                                                                                                   cgccgcCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                   <  1:478328/251‑1 (MQ=255)                                                                                                                                                                                                                    gccgccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                                  >  1:78976/1‑251 (MQ=255)                                                                                                                                                                                                                    gccgccTTCTGCGTCGCTACAAAGGATTCGATTCACCC‑TCGGGCCAGAGCCAACAAGATGAGTAACTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTAATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATCCCAGATCGAAAACAGCTTCGCC‑GGTCACAATACCGAAAGAAACTGCGGCTACCATCCCTTTCATCGGATTggg                                  >  2:346996/1‑250 (MQ=255)                                                                                                                                                                                                                     ccgccTTCTGCGTCGCTACAATGGATTAGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtt                                 >  1:431710/1‑251 (MQ=255)                                                                                                                                                                                                                      cgccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCCGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGGCTTCGGTCATTACTCCGTTCATGTCGGTCTCTGCGGCAGCGTCCTCAAAGTATTCCAGATCGCAACCCACCTCACCTGTCCA‑AATTCCGTCCGCAACTGCTGCTACCATCACTTACATTGGCTTGGttt                                >  2:149704/1‑251 (MQ=255)                                                                                                                                                                                                                      cgccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTGCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                                <  1:521320/251‑1 (MQ=255)                                                                                                                                                                                                                        ccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCCGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCATAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                              >  1:418602/1‑251 (MQ=255)                                                                                                                                                                                                                        ccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTACTCGGATTGGttttt                              >  1:109840/1‑250 (MQ=255)                                                                                                                                                                                                                        ccTTCTCCGTCGCTACACTGGTTTCGATTC‑CCCCTCGGACCAGCGCCAACACGATCATTAGCTCTTCATGCGCGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCCCCCGTCTTCGGTCATCACTCCGTTCATGTCGGACTCTGCGGCCGAGCCTTCAACGTCTTACAGATCGCAAACCGCTTCGCC‑GTTCACCATTCCGACAGCAACGGCGGCTACCATCCCTTTCATAGGATTGAtttt                               >  2:348719/1‑250 (MQ=255)                                                                                                                                                                                                                           tCTGCGTCGCTAAAATGGATCCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCa                             <  2:128350/249‑1 (MQ=255)                                                                                                                                                                                                                            cTGCGTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                          >  1:231089/1‑251 (MQ=255)                                                                                                                                                                                                                             tGCGACGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCCCAATTCCGACAGACACTGCGGCTACCATCCCTTTCTTCGGATTGGTTTTCAACtt                         >  2:402667/1‑251 (MQ=255)                                                                                                                                                                                                                                gTCGCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                      >  1:161877/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCTACAATGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCCAACCGCTTCGCC‑GTTCACACTTCCGACAGAACCTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCCGCTTGCCGtt                   >  2:61314/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaTGGATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAAATGCGGCt                                                         >  2:287004/1‑208 (MQ=255)                                                                                                                                                                                                                                           ggATTCGATTC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacca           <  1:473250/251‑1 (MQ=255)                                                                                                                                                                                                                                                gattcc‑ccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGACCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTAACATCCCTTTCATCTGATTGGTTTTCAGCTTGCCGTTTTCCACCCGCtttt     >  2:174759/6‑249 (MQ=255)                                                                                                                                                                                                                                                   tcc‑ccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGCTTGGGTTTCAGCTTGCCGTTTTACACCAGCTTTTGtc  >  2:385133/3‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCC‑GTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTA  >  NC_000913/3815636‑3816128

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