Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,903	A→G	I216T (ATC→ACC)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,903	0	A	G	98.4%	205.0 / NA	61	I216T (ATC→ACC)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (37/23);  minor base T (0/1);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCAC  >  NC_000913/3815659‑3816147                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                                 <  1:224317/251‑1 (MQ=255)         gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGACAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTc                                                                                                                                                                                                                                         >  1:395653/1‑251 (MQ=255)              aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGCAGCTCTTCATggg                                                                                                                                                                                                                                    >  2:349891/1‑251 (MQ=255)               ggggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATgggg                                                                                                                                                                                                                                   >  1:175350/1‑250 (MQ=255)                 ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                                                 <  1:296712/251‑1 (MQ=255)                    ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                                              >  1:417956/1‑251 (MQ=255)                     tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                               <  2:102553/249‑1 (MQ=255)                          gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                                                                                                                                                                        >  1:384951/1‑251 (MQ=255)                             cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGNAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                                                     >  2:246107/1‑251 (MQ=255)                              ggATTTCAGCGTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                                                                                                                                                                    <  2:262865/251‑1 (MQ=255)                               aaggtccggcgaaaCTCGACAAACTTTAACAGCCGCTCGCTGAGCGGAAAATCAATAATCTTGCGCTGATAGGGTGTCATGCCTTCGCTGCTTATCTTACCTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGGCTTAAAAATCTGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGCGCCAGCGCCAACAAGGTGGTTAGCTCTTCATGGGTGGACGTCGCGCCtact                                                                                                                                                                                                                   <  2:385157/240‑4 (MQ=255)                                aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                                                  >  2:120919/1‑251 (MQ=255)                                aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                                                  <  1:48753/251‑1 (MQ=255)                                                    aaCTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAatg                                                                                                                                                                                              <  2:175350/251‑1 (MQ=255)                                                                    tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                                                <  2:186771/249‑1 (MQ=255)                                                                    tgctAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCAACAGCCAAAAAAAAGGCGACTAATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCTGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCGGCCGTCCCCTGCACTTCGATGATGCGCCCGTCTTCg                                                                                                                                                                               <  2:184425/250‑1 (MQ=255)                                                                        aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                                             <  1:445219/248‑1 (MQ=255)                                                                           cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                                       <  1:349891/251‑1 (MQ=255)                                                                           cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                                       <  1:421522/251‑1 (MQ=255)                                                                            gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                                                      >  1:199131/1‑251 (MQ=255)                                                                              aaaTTCACTAAACAGGCGCTGATATGGTTTCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAACAAGGCGCATCCTCAGTCGCCTTAAAAATAAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAAACAGGTGAGAAGCCCTTAATGGGTGAACGGCCCGCCTTCAGCCGTCCCATGCACTCAAATGATGCCACCGTCTTCGGtcatccccac                                                                                                                                                                    >  2:436232/1‑246 (MQ=255)                                                                                   cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                                                                >  1:407867/1‑250 (MQ=255)                                                                                        aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGCCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCATTCTGCCGTCCCCTGCACTTCAATGAAGCGCCCGTCTTACGTCATCCCTACGTTCATGGCgg                                                                                                                                                         >  2:463087/1‑251 (MQ=255)                                                                                                 tGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                                  <  1:482581/250‑1 (MQ=255)                                                                                                  gATATAGTTTCTTGCCGTCTCTCCACATCTTGTTTTTCTACAGTCTAAAAAAAGGCGACTCCTCAGTCTCCTTAAAAATCCGTTTGCCCGCGCCGCCTGCTGTGTCCCTCCAATTCATTCGA‑GGC‑CCCC‑TCGCGCCAGCGCCAACTAGGTGAGTAGCTCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctcttt                                                                                                                                                  <  2:350526/249‑3 (MQ=255)                                                                                                     atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGACGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTACTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCgtctg                                                                                                                                             >  2:289857/1‑247 (MQ=255)                                                                                                     atGGTTACATGCCTTCGCTTCTCAGCTTCCCTTTCTCCAGACAAAAAAAAGGCGATTCATCAGTCGCTTTAAAAACCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTAAAATGTATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCGTTATGGGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCATGCGGTCATCCCTACTCTCATGTCGGTCTCTTCGGct                                                                                                                                              <  2:381285/250‑2 (MQ=255)                                                                                                               gccctCGCTCCTCATGTTACTTTTCTACAAACAAAAAAAAGGCGACACATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTTCAATGGATTCGC‑T‑T‑CCCCCTCGGGCCAGCGTCAACAAGGTGCTTAGTGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCaat                                                                                                                                   <  2:54483/247‑2 (MQ=255)                                                                                                                     gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTATCAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCCACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATTCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCTTGGTCTTCAACGTCTTc                                                                                                                             >  2:156377/1‑251 (MQ=255)                                                                                                                        cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                            >  1:366602/1‑249 (MQ=255)                                                                                                                           cATCTGACTTCTCTAACGACAAAAAAAAGGCGACACAGCAGTCGCCTTAAAAATCACTTTGCCGGCGCCGCCTTCTGCGTCGCTACAATGGATTCGACT‑C‑CCCC‑GCTGGCCAGCGCCAACCAGTTGTGTAGCTCTTCATTGTTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGTCTGCGTCTTCAACGTATTCCAGagc                                                                                                                       <  2:20196/251‑3 (MQ=255)                                                                                                                                tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaa                                                                                                                   >  2:252826/1‑250 (MQ=255)                                                                                                                                 aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                                  <  1:137367/250‑1 (MQ=255)                                                                                                                                    tttCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaaattgc                                                                                                              >  2:161382/1‑247 (MQ=255)                                                                                                                                               aaataaaaGGCGACTAATCAGTCGCCTTTCACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TTCACA‑TCAGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTCCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCATGTCGGTCTCGGCGGCGGCGTCTGCCACGTATCCCAGATCGCaatttgcttcaacgttt                                                                                                    >  2:28341/5‑235 (MQ=255)                                                                                                                                                  aaaaaGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTAATCACTACGTTCATGTCGGTCTCTGCGCCTGAGTCTTCAACGTATTCCAGCTCGctaatttcttcaccgttttct                                                                                                 >  2:256339/1‑230 (MQ=255)                                                                                                                                                     atggCGATTCATCAGTAACCTTAAAGATCAGTTTGACAGCCACGAATTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCTGAGTCTTCAACGTATTTCAGCTCGCAAATTGCTTCCCCGTTTACAAtac                                                                                             >  2:273046/3‑249 (MQ=255)                                                                                                                                                        gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                           >  1:350878/1‑250 (MQ=255)                                                                                                                                                        gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTAAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCCTCACTACGTTCATTTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                           >  1:245359/1‑250 (MQ=255)                                                                                                                                                              catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGa                                                                                      <  2:377798/249‑1 (MQ=255)                                                                                                                                                                 caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                                 >  1:355915/1‑251 (MQ=255)                                                                                                                                                                        cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTTCGGCTAc                                                                          >  1:413564/1‑251 (MQ=255)                                                                                                                                                                            aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATc                                                                      >  1:284774/1‑251 (MQ=255)                                                                                                                                                                                 cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCCCCGTTTACAATACCCACAGAAACTGCGGCTACCATCCCttt                                                                 >  2:41603/1‑251 (MQ=255)                                                                                                                                                                                       gccacccccgccTTCTGCGTCTCTACAATGGATTCGAGT‑C‑CCCC‑TCGGGCCAGCGCCAACAAGGTGGGGAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGGCCGGCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTGCACCGTTTACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                                            <  1:273046/243‑1 (MQ=255)                                                                                                                                                                                         cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                                   <  1:344277/97‑1 (MQ=255)                                                                                                                                                                                         cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                                   >  2:344277/1‑97 (MQ=255)                                                                                                                                                                                            cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCACGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTgg                                                      >  1:102560/1‑251 (MQ=255)                                                                                                                                                                                             gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                     >  1:451433/1‑251 (MQ=255)                                                                                                                                                                                             gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                     <  1:3097/251‑1 (MQ=255)                                                                                                                                                                                             accgccTTCTGCGTCGATACAATGGATTCGA‑T‑T‑CCCCCTCGGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATTATGCGCCCATCTTCGGTCATCACTACGTTCATGTCGATCTCTTCGGCTGCGCCTTCAACGTATTCCCGCTCGCACATTGCTTCACCTTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCggttttgt                                                     >  2:63/2‑245 (MQ=255)                                                                                                                                                                                                gccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGtttt                                                  >  1:58825/1‑251 (MQ=255)                                                                                                                                                                                                    tcggCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGGGCCAGCGTCAACAAGGTGAGTAGTTCGTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAg                                               <  2:333823/247‑1 (MQ=255)                                                                                                                                                                                                     cTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTAAGCt                                             >  2:484868/1‑251 (MQ=255)                                                                                                                                                                                                       gCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATAGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTg                                           >  2:54780/1‑251 (MQ=255)                                                                                                                                                                                                        cGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACTGCGGCTACCATCCCCTTCATCGGATTGGTGTTAAGCTTGc                                          >  2:441607/1‑251 (MQ=255)                                                                                                                                                                                                         gTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTAAGCtttcc                                         >  2:336014/1‑248 (MQ=255)                                                                                                                                                                                                                                         ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCTTCAACGTATTCCAGATCGCCAATTGCTTCACCGTTTACAATACCGGCCGAAACTGCGGCTACCATCCCTTTCCTCGGGTTGGTTTTTAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg            >  2:3875/1‑250 (MQ=255)                                                                                                                                                                                                                                               cgCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcg      <  2:384951/250‑1 (MQ=255)                                                                                                                                                                                                                                                 ccAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGcc    >  2:343680/1‑251 (MQ=255)                                                                                                                                                                                                                                                   aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc  >  1:359826/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     GCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCAC  >  NC_000913/3815659‑3816147

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