Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,164,563	G→A	G157R (GGA→AGA)	nagZ →	beta‑N‑acetylglucosaminidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,164,563	0	G	A	98.4%	206.9 / NA	63	G157R (GGA→AGA)	nagZ	beta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/27);  minor base C (0/1);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCCGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGG  >  NC_000913/1164318‑1164806                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    gtACCCGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGTCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAAt                                                                                                                                                                                                                                                 <  1:439372/249‑1 (MQ=255)   aCCCGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGa                                                                                                                                                                                                                                               >  1:110073/1‑250 (MQ=255)    cccGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGaaa                                                                                                                                                                                                                                             <  1:224550/251‑1 (MQ=255)    cccGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTGCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGaaa                                                                                                                                                                                                                                             <  2:547228/251‑1 (MQ=255)     ccGCTTGCCAGCGGCGTAATCATTAGCTGCGCTGTCAGGAATGGTAGAGGAGGGCATATTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAACTGATCGATATGGATCGTGATAAAAGCTTTGCGCCCGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTAATTGATGGTATGCATGAAGCCAGAATGaaa                                                                                                                                                                                                                                             >  2:612343/1‑250 (MQ=255)          tgcaaGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGaa                                                                                                                                                                                                                                              >  2:146055/5‑244 (MQ=255)          tGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATTTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGaa                                                                                                                                                                                                                                              <  1:146055/244‑1 (MQ=255)          tGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGAcc                                                                                                                                                                                                                                       >  1:396901/1‑251 (MQ=255)          tGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGCCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGAc                                                                                                                                                                                                                                        <  1:111342/250‑1 (MQ=255)                           ttCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAgg                                                                                                                                                                                                                      >  1:220487/1‑251 (MQ=255)                            tCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAgg                                                                                                                                                                                                                      <  2:76087/250‑1 (MQ=255)                             cgctgcgctgTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGac                                                                                                                                                                                                                    <  1:310659/251‑1 (MQ=255)                             cgctgcgctgTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGac                                                                                                                                                                                                                    <  2:648657/251‑1 (MQ=255)                                 gngctgTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGacac                                                                                                                                                                                                                  <  2:606441/249‑1 (MQ=255)                                           gAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACg                                                                                                                                                                                                      >  2:321314/1‑251 (MQ=255)                                             aTGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGc                                                                                                                                                                                                    >  1:656383/1‑251 (MQ=255)                                                     gggTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcac                                                                                                                                                                                            <  1:492338/251‑1 (MQ=255)                                                       gTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcacac                                                                                                                                                                                          >  2:506170/1‑251 (MQ=255)                                                              aCTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGaaa                                                                                                                                                                                   <  2:103257/251‑1 (MQ=255)                                                                    gcAGGAGGAAGGTTGGTTGAGGGCCAGCGAAATGATCGCTATGGATAGTGATATCAGCTTTGCGCCTGTGCGGGATGTCGGGCATAGCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGt                                                                                                                                                                             <  2:179997/251‑1 (MQ=255)                                                                                ttggttgATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTc                                                                                                                                                                 >  1:429677/1‑251 (MQ=255)                                                                                               cGAAATGATCGCTATGGATATTGATATCAGCTTGGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTc                                                                                                                                                  <  2:249117/251‑1 (MQ=255)                                                                                                gAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcg                                                                                                                                                 >  2:259050/1‑251 (MQ=255)                                                                                                  aaTGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCAGTAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACGCAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcgcg                                                                                                                                               <  1:653337/251‑1 (MQ=255)                                                                                                     gATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCt                                                                                                                                             >  1:73528/1‑250 (MQ=255)                                                                                                      aTCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTa                                                                                                                                            >  1:307056/1‑250 (MQ=255)                                                                                                         gCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCTTATCGGCATGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACCAAAAGCACTGGCAATTGCCAGCCGGTTTATTGATGGGATCCATGAAGCCAGAATGAAAACGACCGGGAAAGACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTaa                                                                                                                                           >  1:64630/1‑248 (MQ=255)                                                                                                              ggATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGAATCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGt                                                                                                                                   >  2:355857/1‑251 (MQ=255)                                                                                                                atatTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCg                                                                                                                                 <  1:67311/251‑1 (MQ=255)                                                                                                                   ttGATATCAGCTTTGCGCCTGAGCTGGATGTCGGGCATATCAGCGCGGAGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTCTTGCTGGTATGCATGAAGCCAGCATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACAACGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTc                                                                                                                              >  2:115343/1‑251 (MQ=255)                                                                                                                     gATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCtt                                                                                                                            >  1:50221/1‑251 (MQ=255)                                                                                                                        atCAGCTGTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCGGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCa                                                                                                                          <  1:106316/250‑1 (MQ=255)                                                                                                                          cAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGtt                                                                                                                       >  2:221658/1‑251 (MQ=255)                                                                                                                                 gcgcCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAAt                                                                                                                <  2:214512/251‑1 (MQ=255)                                                                                                                                     cTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGACAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCcgc                                                                                                            >  2:560260/1‑251 (MQ=255)                                                                                                                                       gtgCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCcgcg                                                                                                           >  2:107906/1‑250 (MQ=255)                                                                                                                                            ggATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAATATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAAt                                                                                                      <  2:397706/250‑1 (MQ=255)                                                                                                                                                gTCGGCCATATCAGCGCTGCAATTGGCTAGCGGTCTTATCATGTCAATCCATAAAAAGCCCTGGCAATTGCCTGCCTGTTTATTGATGGTATGCAGGAAGTCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTTCAATAACTGCATACTCACACAAAGAAACCCCTTGCGATTCACGTCCTCAAGCGGAGATTCGCGCTACAGATATGTCGGTCTTCAGTTCCGTAAGCCGCAAAAATAAAc                                                                                                  <  2:539089/250‑1 (MQ=255)                                                                                                                                                   gggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCg                                                                                               >  1:531356/1‑250 (MQ=255)                                                                                                                                                      cATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGc                                                                                           >  2:620397/1‑251 (MQ=255)                                                                                                                                                        tatCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGcc                                                                                          >  2:514301/1‑250 (MQ=255)                                                                                                                                                            aGCGCGGCGTTTGTCGAGCGTTCTTATCAAGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAATCACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTTCGATCCACGTCCACAAGCGGAGATTCGCGCTAAACATATGTCGGTTTTCAGTTCCTTAATCCGCGAAAATAAACTCGTCTCCATTAt                                                                                     <  2:316845/251‑1 (MQ=255)                                                                                                                                                            aGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTAt                                                                                     >  1:275084/1‑251 (MQ=255)                                                                                                                                                              cggggcgGTTGGGGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATGCGCGCGAAAGATATGTCGGTCTTCAGTGCCGTAAACCGCGa                                                                                                          <  2:266929/225‑1 (MQ=255)                                                                                                                                                              cgcgGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCACAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAACGATATGTCGGTCTTCAGTTCCTTAATCCGCGa                                                                                                          >  1:266929/1‑228 (MQ=255)                                                                                                                                                                 ggcgATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTg                                                                                >  2:265598/1‑251 (MQ=255)                                                                                                                                                                  gcgATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGAAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACTCCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACGCGACGCCATTATGCCTgc                                                                               <  2:220487/251‑1 (MQ=255)                                                                                                                                                                    gATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTgcgc                                                                             <  2:366161/251‑1 (MQ=255)                                                                                                                                                                              cGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCt                                                                    >  1:112386/1‑250 (MQ=255)                                                                                                                                                                                 tCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGaaa                                                                                                        >  1:125935/1‑211 (MQ=255)                                                                                                                                                                                     aTCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGAt                                                            >  1:116777/1‑251 (MQ=255)                                                                                                                                                                                       cATGCCGATCCCCAAAAAGCCCTGGCAATTGCCAGCCGGTATATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGAGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGCTATGTCGATCTTCAGCTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGACTGCGCATGGGATCTACAGTGATGt                                                          >  2:243418/1‑251 (MQ=255)                                                                                                                                                                                            cGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATc                                                     <  1:506170/251‑1 (MQ=255)                                                                                                                                                                                                    aaaaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCg                                             >  2:176355/1‑251 (MQ=255)                                                                                                                                                                                                                           cGGTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCt                      <  2:301965/251‑1 (MQ=255)                                                                                                                                                                                                                             gTTTATTGATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTg                     >  2:483068/1‑250 (MQ=255)                                                                                                                                                                                                                                 aTTNATGGTATGCATGAAGCACGCATGAAAACGAACGGGAACCACTTCCCATGACACGGTGCAGTAACGGCAGAATCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGTTCTTCAGTTCCTTAATCCGAGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTCCTGGCTGaaacc                <  2:584876/251‑3 (MQ=255)                                                                                                                                                                                                                                    gATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGCCACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGt             >  1:86339/1‑251 (MQ=255)                                                                                                                                                                                                                                    gATGGTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGt             >  2:232563/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTATGCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGTAGACTCACACAAAGAAACCCCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACGCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCAGCGTCCGGCGAGCGGTTATCCCTACTGGCTGAAAACCGTTTTg         <  1:302848/251‑1 (MQ=255)                                                                                                                                                                                                                                            gCATGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTc     >  1:93388/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTGAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTc     >  1:30215/1‑249 (MQ=255)                                                                                                                                                                                                                                                gAAGCCAGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg  <  2:336448/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    TTACCCGCTTGCCAGCGGCGCAATCATTCGCTGCGCTGTCAGGAATGGAAGAGGGTGGCAAACTGGCGCAGGAGGCAGGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGG  >  NC_000913/1164318‑1164806

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