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

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

CGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815638‑3816132                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          cgcCCGGTATGAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCGGATATGGTTTCATGCCTTCGCTCCTCATCTTTCGGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAGTCGATT‑c‑c‑ccctcgg                                                                                                                                                                                                                                                       <  2:76692/251‑6 (MQ=255)       gTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGcc                                                                                                                                                                                                                                                    >  1:175650/1‑248 (MQ=255)                    cGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGa                                                                                                                                                                                                                                    <  1:216628/251‑1 (MQ=255)                            aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                            >  1:249290/1‑251 (MQ=255)                                   aGGGGCTTTTGCGCCCGGATTTGAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTTGCGCTGATGTGGTTTCATGCCTTCGCGCCTCAGCTGACTTATCTACAGACAAAAAAAAGGCGACTCAGCAGGCGCCTGAAAAATCAGTTGGCCAGCGCCGCCTTCTGCGTCGCGACAATGGATTCGATG‑C‑CGCC‑TCGGGCTAGAACCAGCAAGATGTGTAGCTCTTCATggt                                                                                                                                                                                                                     <  2:244660/251‑2 (MQ=255)                                   aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                                                                                                                                                      <  1:99947/250‑1 (MQ=255)                                       gCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                                  >  2:97954/1‑250 (MQ=255)                                          tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                               <  2:67143/250‑1 (MQ=255)                                           ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACACAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGCGCCAACAAGATGCGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                             >  2:245532/1‑251 (MQ=255)                                              cgcCCGTATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATTTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCTCCGCCTTCTGCGTCGCGACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                                                                                                                                          <  2:227928/251‑1 (MQ=255)                                                 ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATAAGTAGCTCTTCATGGATGAACGGCTCGCCt                                                                                                                                                                                                       >  2:172653/1‑251 (MQ=255)                                                       ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 >  1:40744/1‑251 (MQ=255)                                                             gTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACGGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCGTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                           <  2:101315/251‑1 (MQ=255)                                                                            tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                            >  1:175874/1‑251 (MQ=255)                                                                                    cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                                                                    >  2:4406/1‑251 (MQ=255)                                                                                     tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGCTT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATAAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCAGCCGTCCCCTGCACTTCAatcaagaccccctct                                                                                                                                                                   >  2:215130/1‑238 (MQ=255)                                                                                         tgctAAGCGCAAATTCAATAAACGGGCGCTGATATGGTGTCATGCCTTCGCTCCTCATCCTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACTGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGTCTTc                                                                                                                                                                 <  2:245993/249‑1 (MQ=255)                                                                                            tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                <  2:113347/199‑1 (MQ=255)                                                                                            tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                >  1:113347/1‑199 (MQ=255)                                                                                                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                                       >  1:226131/1‑251 (MQ=255)                                                                                                    aaTTCAATAAACTGTCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGTATTCGATT‑C‑C‑CCCTCGGGCCATAGCCAACAAGATGAGTAGCTCTTCATGGTTGAACGGCGCGTCTTCTGCCGTCCCCTGCACTTCAATGATGTGCCCGTCTTCGGTCATCACTACg                                                                                                                                                    <  2:130568/251‑1 (MQ=255)                                                                                                                 ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                                                                                                                        <  1:231019/250‑1 (MQ=255)                                                                                                                             gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                                                                                           >  1:35265/1‑251 (MQ=255)                                                                                                                                    gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                                                                     <  2:175874/250‑1 (MQ=255)                                                                                                                                                  tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CAC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGc                                                                                                      >  1:173953/1‑251 (MQ=255)                                                                                                                                                      aCTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTTCCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                   <  1:127994/250‑1 (MQ=255)                                                                                                                                                         tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg                                                                                                <  1:49838/250‑1 (MQ=255)                                                                                                                                                               cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGATCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGcc                                                                                         >  1:67330/1‑251 (MQ=255)                                                                                                                                                                  aCAAAAAAAAGGCGCCTCATCAGTCGACTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGGGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCCGAGTCTTCAACGACTTCCATATCGCATACCGCTTCGCCggt                                                                                      >  2:114826/1‑249 (MQ=255)                                                                                                                                                                     aaaaaaaCGCGACTCATCAGTAGCCTTACAAATCGGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTAC‑C‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCCCGCCTTCTGCCGTCCCCGGCGCTTGAATGGTGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCCGAGCCTTCAACGTATTCCATGTCGCCACCAGCTTCGCCGTCcac                                                                                   >  2:50608/1‑251 (MQ=255)                                                                                                                                                                       aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACaa                                                                                 <  1:242462/251‑1 (MQ=255)                                                                                                                                                                        aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAt                                                                                >  1:61587/1‑251 (MQ=255)                                                                                                                                                                         aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtt                                                                               >  1:49867/1‑251 (MQ=255)                                                                                                                                                                            ggCGACTCATCAGTTGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCATCGGGACAGAGACAACAAGATGAGTAGCTCTTCATGGGCGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCCCTCCGTTCATGTAGGTCTCTGCGGCCGAGTCTTCAACGTATTCCCGACCGCACACCGCTTACCCGTTCACAATTccc                                                                            >  2:253554/1‑250 (MQ=255)                                                                                                                                                                              cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  1:252748/1‑251 (MQ=255)                                                                                                                                                                                 cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  1:204991/1‑248 (MQ=255)                                                                                                                                                                                   catcaGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTTCTGAGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTGCATGGGTGAACGGCTCGCCTGCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTTATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGa                                                                       <  1:140822/249‑1 (MQ=255)                                                                                                                                                                                            ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAATCCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTa                                                            <  1:172653/251‑1 (MQ=255)                                                                                                                                                                                                  aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCACGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc                                                      >  2:55086/1‑251 (MQ=255)                                                                                                                                                                                                     tCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt                                                   >  1:234577/1‑251 (MQ=255)                                                                                                                                                                                                                  gccgccTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTCCCATCCCTTTCATCGGATTGGt                                      >  2:31484/1‑251 (MQ=255)                                                                                                                                                                                                                   ccgccTTCTGCGGCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAGTCCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGATTGGtt                                     <  1:157056/251‑1 (MQ=255)                                                                                                                                                                                                                     gccTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCAGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                   >  1:30205/1‑251 (MQ=255)                                                                                                                                                                                                                         tCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGc                               >  1:242577/1‑251 (MQ=255)                                                                                                                                                                                                                          cTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                              >  2:55969/1‑251 (MQ=255)                                                                                                                                                                                                                            gCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTg                            >  1:13321/1‑251 (MQ=255)                                                                                                                                                                                                                             cGTCGCTACAATGGATTCGATT‑C‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                           >  1:5691/1‑251 (MQ=255)                                                                                                                                                                                                                                             tcgatt‑c‑C‑CCCTCGGGCCAGAGCCAACAAGAGGAGTAGCTCGTCATGGGTTAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATACGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAATCCGCTTCGCTGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCGTTCATCGGATTGGTTTTCAGCTTGCCGTTGTCCACCAGCt            <  2:234538/250‑1 (MQ=255)                                                                                                                                                                                                                                                  tcc‑c‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCAGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATCCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  1:197803/3‑250 (MQ=255)                                                                                                                                                                                                                                                        cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  1:275066/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          CGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815638‑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: ‑