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.6%	141.0 / 5.3	48	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (27/17);  total (31/17)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-02

ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815646‑3816132                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          aTTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTGCGCTCCTCATCTTACGTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTACAATGGATTCGATT‑c‑c‑cccgcggtccagag                                                                                                                                                                                                                                                <  2:163817/250‑13 (MQ=255)        aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCcaaca                                                                                                                                                                                                                                          >  2:293924/1‑249 (MQ=255)         gCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGa                                                                                                                                                                                                                                       <  1:20982/251‑1 (MQ=255)                   gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                                                                                                                                                                              <  2:299647/250‑1 (MQ=255)                        aaTAGGTGCTTTTGAGCCCGGATTTCAGCGTAAACTCGCCAAAATTTAAAACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCAATGGATTCGATT‑C‑C‑CCCTCGGGCGAAAGCCAACAAGATGAGTACATCTTCAt                                                                                                                                                                                                                        >  2:371318/1‑251 (MQ=255)                                            gttgtCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGATCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                    <  2:239605/247‑1 (MQ=255)                                                tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 <  1:177912/250‑1 (MQ=255)                                                  agcggaaaccagCCAAACTATAACACCTGCTTGCTAAGCGCAAATACAATAAACTGGCGCTGATATGGTCTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCTCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACTGCTCGCCTTCTTCCGTc                                                                                                                                                                                              <  2:243315/240‑1 (MQ=255)                                                                    tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                            >  1:281532/1‑251 (MQ=255)                                                                        acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                                                                        >  1:55887/1‑251 (MQ=255)                                                                          acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                                                                                                                                                         >  1:131745/1‑216 (MQ=255)                                                                          acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                                                                                                                                                         <  2:131745/216‑1 (MQ=255)                                                                             tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCt                                                                                                                                                                   >  1:121061/1‑251 (MQ=255)                                                                                 tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCGCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                                                                                <  2:55887/250‑1 (MQ=255)                                                                                       gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTACTCATCTTACTTTTCTACAGACAAAAAAAAGACGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCAATTCAATGATGCGCCCGTCTTCGGtcatcc                                                                                                                                                         >  1:59088/1‑250 (MQ=255)                                                                                        cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACATAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                        >  1:260808/1‑251 (MQ=255)                                                                                           aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCAAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTAc                                                                                                                                                     >  2:28160/1‑251 (MQ=255)                                                                                            aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGTGCCCGTCTTCGGTCCTCAGTACg                                                                                                                                                    >  2:182479/1‑251 (MQ=255)                                                                                                         ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACAGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                                                                       >  1:50793/1‑251 (MQ=255)                                                                                                                      tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                                                                          >  1:254778/1‑251 (MQ=255)                                                                                                                                 cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                                                                               <  1:206970/251‑1 (MQ=255)                                                                                                                                   ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                             >  1:137583/1‑251 (MQ=255)                                                                                                                                    tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CAC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTcc                                                                                                             >  2:72669/1‑250 (MQ=255)                                                                                                                                      ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCACCGCCTTCTGCGTCGCTACAATGGATTCGATTAC‑C‑CC‑TCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                             >  1:381701/1‑248 (MQ=255)                                                                                                                                            ttACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCATCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa                                                                                                     >  2:250023/1‑250 (MQ=255)                                                                                                                                                      aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGc                                                                                          >  1:179014/1‑251 (MQ=255)                                                                                                                                                                     gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                                                                            >  1:336232/1‑250 (MQ=255)                                                                                                                                                                     gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                                                                           <  2:96566/251‑1 (MQ=255)                                                                                                                                                                       gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACa                                                                         <  1:157107/251‑1 (MQ=255)                                                                                                                                                                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAg                                                                        <  2:281532/251‑1 (MQ=255)                                                                                                                                                                                           aaTCAGGTTGCCTGCGCCGCCGTCTGCGTCGCTACAATGGATCCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCGTCAACGTATTCCAGATCGCAAACCGCTGCGCCGTTCACAATTCCGACAAAAACTGCGGCTACCATccc                                                     <  2:87657/251‑1 (MQ=255)                                                                                                                                                                                            aTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCGCAACCCGCTTCGCCGTTCACCATTCCGACAGAAACTGCGGCTACCATCCCt                                                    >  2:18839/1‑251 (MQ=255)                                                                                                                                                                                               tgTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCACCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGAGGCAGAGTCTTCACCGTATTCCAGATCGCAAACCGATTCGCCGTTCACAATTCCGACAGAAACCGCAGCTACCCTCCCtttt                                                 >  2:83514/2‑250 (MQ=255)                                                                                                                                                                                                 tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                             >  1:273218/1‑189 (MQ=255)                                                                                                                                                                                                 tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCAGCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                             <  2:273218/189‑1 (MQ=255)                                                                                                                                                                                                   tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                             >  1:150366/1‑251 (MQ=255)                                                                                                                                                                                                     ccAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTCCCATCCCTTTCATCggg                                           >  2:374891/1‑250 (MQ=255)                                                                                                                                                                                                          gccgccTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                                      >  1:200437/1‑251 (MQ=255)                                                                                                                                                                                                            cgccTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGttt                                    >  1:104249/1‑251 (MQ=255)                                                                                                                                                                                                                 tCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAg                                <  2:101090/250‑1 (MQ=255)                                                                                                                                                                                                                   tGCGTCGCTCCAATGGATTCGCTT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtt                             >  2:74524/1‑251 (MQ=255)                                                                                                                                                                                                                    gCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGTCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGc                               <  2:114586/248‑1 (MQ=255)                                                                                                                                                                                                                     cGTCGCTACAATGGATTCGATT‑C‑CTCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                           >  1:99933/1‑251 (MQ=255)                                                                                                                                                                                                                     cGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                           >  1:48037/1‑251 (MQ=255)                                                                                                                                                                                                                      gTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                          >  1:105728/1‑251 (MQ=255)                                                                                                                                                                                                                         gCTACAATGGATGCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtt                       <  1:153575/251‑1 (MQ=255)                                                                                                                                                                                                                              aaTGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccc                  >  1:160596/1‑250 (MQ=255)                                                                                                                                                                                                                                   attcgagtcc‑c‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTGTTCATGGGTGAACGGTTCGCCTTCTGCCGTCCCCTTCACTTCAATGAGGCGCCCGTCTTCTGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGTCGTGTACAATTCCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGATTGGTTTTTTTCTTGCCGTTTCCCACCAg              <  2:50793/241‑1 (MQ=255)                                                                                                                                                                                                                                    ttcgatt‑c‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg              <  2:396712/249‑1 (MQ=255)                                                                                                                                                                                                                                    ttcgatt‑c‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg              >  1:396712/1‑249 (MQ=255)                                                                                                                                                                                                                                      cgactcc‑c‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTCCGTTCATGTCTGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGAGTGGTTTTCAGCTTGCCGTTTTCCACCAGCtt           <  1:112505/244‑1 (MQ=255)                                                                                                                                                                                                                                       gattcc‑c‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCt         >  2:353645/6‑251 (MQ=255)                                                                                                                                                                                                                                       gattcc‑c‑cc‑GCGGGCCATAGCCAACAAGATGAGTAGCGCTTCATGGTTGACCGGCTCGCCTTCTTCCGTCCCCTTCACTTCAATGATGCGCCCGTCTCCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTTGTCTACGTAGGTCAGATCGCAAACCGCTTCGCCTTTCACAATTCCGACAGAAACGGCGGCTCCCAGCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCttat         <  2:337759/246‑3 (MQ=255)                                                                                                                                                                                                                                                cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  1:13416/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815646‑3816132

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