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	91.2%	197.8 / 10.6	69	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base C (36/26);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCT  >  NC_000913/3815637‑3816122                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 gcgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGTTTTTGCGCCCTGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTGCAATGGATTCGATT‑c‑c‑ccctcg                                                                                                                                                                                                                                               <  1:133365/251‑5 (MQ=255)  cgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑c‑c‑ccctcgg                                                                                                                                                                                                                                              >  2:38625/1‑246 (MQ=255)     ccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCTCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAAGCCGCCTTAAAAATCAGTTTCCCAGCGCATCCTTCTTCGTCGCTACAATGAATTCGATT‑C‑CACC‑TC‑gcgcg                                                                                                                                                                                                                                          >  2:316438/1‑247 (MQ=255)        gTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCaga                                                                                                                                                                                                                                       <  1:286515/251‑1 (MQ=255)                 aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCcaaca                                                                                                                                                                                                                                >  1:22091/1‑249 (MQ=255)                      ggCGTTGCAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                           >  1:421491/1‑249 (MQ=255)                                    aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                                                                                                                                            <  1:406749/250‑1 (MQ=255)                                                      aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                         <  1:351756/251‑1 (MQ=255)                                                                aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTACAAGGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCt                                                                                                                                                                                <  1:270291/250‑1 (MQ=255)                                                                     cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCCATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                          >  1:398593/1‑251 (MQ=255)                                                                       ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCCATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                         >  1:289141/1‑250 (MQ=255)                                                                         aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                      >  2:362067/1‑251 (MQ=255)                                                                         aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                      >  2:206636/1‑251 (MQ=255)                                                                          aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                                                                                                                                                                     >  2:91377/1‑251 (MQ=255)                                                                                 acacCTGCTTGCTAAGCGCAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATTATGCGcc                                                                                                                                                              <  2:158242/251‑1 (MQ=255)                                                                                  cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                                                             >  1:449042/1‑251 (MQ=255)                                                                                     cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                                                          >  1:353479/1‑251 (MQ=255)                                                                                        cttgctAAGCGCAAATTCAATAAACTGGCGCTGCTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                       >  1:205782/1‑251 (MQ=255)                                                                                         ttgctATGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                                                                      >  1:129138/1‑251 (MQ=255)                                                                                          tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                       <  1:38625/249‑1 (MQ=255)                                                                                                  gcAAAATCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                             >  1:366494/1‑251 (MQ=255)                                                                                                    aaaTTCAATAACCGGGCGCTGATATGGTTTCATGCCTCCGCTCCTCATCTTACTTGTCTACAGAAAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑CCCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAGGGGTGAACGGCTCGCTTTCTTCCGGCCCCTGCTCTGCAATGATGCGCCCGTATTCGGTCATtactac                                                                                                                                           <  2:431793/251‑1 (MQ=255)                                                                                                      aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                                                                         >  1:335690/1‑251 (MQ=255)                                                                                                      aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCATCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGTGCCCGTCTTCGGTCATCACTACGt                                                                                                                                         >  2:422806/1‑251 (MQ=255)                                                                                                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGATCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATTCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                                                                   <  1:462138/250‑1 (MQ=255)                                                                                                                                 tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGTTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                                                              >  2:106719/1‑251 (MQ=255)                                                                                                                                           gatcctcATCTTACTTTTCTACCGACAAAAAAAAGACGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCATCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                                                                    >  2:308689/3‑251 (MQ=255)                                                                                                                                              cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                   >  1:178792/1‑249 (MQ=255)                                                                                                                                                tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                   >  1:169549/1‑247 (MQ=255)                                                                                                                                                       aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGACTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CAC‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc                                                                                        >  1:42645/1‑251 (MQ=255)                                                                                                                                                          tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                                                                       <  1:136685/249‑1 (MQ=255)                                                                                                                                                          tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                                                                       <  2:129652/249‑1 (MQ=255)                                                                                                                                                          tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGc                                                                                     >  2:124712/1‑251 (MQ=255)                                                                                                                                                             cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                                                                  >  1:190064/1‑251 (MQ=255)                                                                                                                                                               aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                                                                  <  2:398593/249‑1 (MQ=255)                                                                                                                                                                 aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                                                              >  2:140877/1‑251 (MQ=255)                                                                                                                                                                 aGACAAAAAAAAGGCGACGCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCTCTACAATGGATTCGATT‑C‑C‑CCCGC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCTGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGTTCGCAAACTGCTTCGCCg                                                                              <  2:421491/251‑1 (MQ=255)                                                                                                                                                                      aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                         >  2:380246/1‑251 (MQ=255)                                                                                                                                                                        aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcaca                                                                        >  1:180609/1‑250 (MQ=255)                                                                                                                                                                            aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                                                                >  2:22198/1‑222 (MQ=255)                                                                                                                                                                            aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                                                                <  1:22198/222‑1 (MQ=255)                                                                                                                                                                            aGGCGACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCTGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCTCGTCTTCGGTCATCACTACGTTCATTTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCCTCGCCGTTCACAATTcc                                                                   <  2:120641/251‑1 (MQ=255)                                                                                                                                                                                gACTCATCAGCCGCCTTAACAATCAGTTTGACAGAGCCGAATTATGCGTCGCTACAATGGATTCGATT‑C‑CTCC‑TC‑GGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                >  2:47391/1‑250 (MQ=255)                                                                                                                                                                                   tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGa                                                             >  1:459759/1‑250 (MQ=255)                                                                                                                                                                                         gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg                                                      <  2:190064/251‑1 (MQ=255)                                                                                                                                                                                           cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                     <  1:23821/250‑1 (MQ=255)                                                                                                                                                                                           cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGc                                                    >  2:157847/1‑251 (MQ=255)                                                                                                                                                                                            gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                   >  2:354991/1‑251 (MQ=255)                                                                                                                                                                                             ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                                                                                                <  2:232137/173‑1 (MQ=255)                                                                                                                                                                                             ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                                                                                                >  1:232137/1‑173 (MQ=255)                                                                                                                                                                                             ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                   >  1:377893/1‑250 (MQ=255)                                                                                                                                                                                                       cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCttt                                        >  2:157710/1‑251 (MQ=255)                                                                                                                                                                                                       cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGACTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCttt                                        >  2:328897/1‑251 (MQ=255)                                                                                                                                                                                                        aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTc                                       <  2:311394/251‑1 (MQ=255)                                                                                                                                                                                                         gTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGCTTAC‑C‑AC‑TC‑GGGCCAGAGCCAACAAGATGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGACTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTc                                       >  1:413141/1‑250 (MQ=255)                                                                                                                                                                                                           tttccAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGATT‑C‑C‑CCATC‑GGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCATCACCTGCACTTCAATGCTACGCCCGTCTTCGGACATCACTAAGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTGCACAATTCCGACAGAAACTGCGGCGACCATCCCttttctc                                    >  2:355800/4‑247 (MQ=255)                                                                                                                                                                                                            tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAGCGCAAACCGCTTCGCCGTTCACAAGTCCGACAGAAACTGCGGCTACCATCCCTTTCATc                                    <  1:43780/250‑1 (MQ=255)                                                                                                                                                                                                              ccAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCgg                                  >  2:343146/1‑250 (MQ=255)                                                                                                                                                                                                                  cgccgcCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                             <  2:22091/251‑1 (MQ=255)                                                                                                                                                                                                                   gccgccTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                            >  2:351493/1‑251 (MQ=255)                                                                                                                                                                                                                   gacgccTTCAGAGTCGCTACAATGGATTAGCTC‑C‑C‑CCCTC‑GGGCCAGGGCCAACAAGATGAGTAGCTCTTCATAGGTGAACGGCTCGCCTTAAGCCGTACAATGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGACGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCCCACACAGCTTCGCCGTTACCAATTCCGACCGAAACTGCGGCTCACATCACTTTCATCGGATTGGt                            >  2:448373/3‑251 (MQ=255)                                                                                                                                                                                                                    ccgccTTCTGCGTCGCTACAATGGATTCGATT‑CAC‑CC‑TC‑GAGCAAGAGGCAACAAGAGGAGTAGCTCTTCATGGGTGAAAGGCTCGCATTCTGCCGTCCCCTGCACTTCAATGATGCCCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGCATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTACGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtt                           >  1:220944/1‑251 (MQ=255)                                                                                                                                                                                                                             gCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGc                     <  2:342619/248‑1 (MQ=255)                                                                                                                                                                                                                             gCGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTg                  <  2:353479/251‑1 (MQ=255)                                                                                                                                                                                                                              cGTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCGTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAGGCCGGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTTCGGCTGCCAGCCCTTTCATCGGATTGGTTTTCAGCTTGc                 <  1:28351/251‑1 (MQ=255)                                                                                                                                                                                                                               gTCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGCTTGGTTTTCAGCTTGcc                >  2:276162/1‑251 (MQ=255)                                                                                                                                                                                                                                tCGCTACAATGGATTCGATT‑C‑C‑CCCTC‑GTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCTTCTTCGGTCATCACTACGTTCGTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                <  1:356600/250‑1 (MQ=255)                                                                                                                                                                                                                                     aCAATGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTc          <  1:180073/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTGGATTCGATT‑C‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccac       <  1:276162/251‑1 (MQ=255)                                                                                                                                                                                                                                           gattcgatt‑c‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACTTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg    <  1:210641/251‑1 (MQ=255)                                                                                                                                                                                                                                             ttcgatt‑c‑C‑CCCTC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCt  >  2:357290/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 GCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑C‑CC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCT  >  NC_000913/3815637‑3816122

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