Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,239,343	C→T	T155T (ACC→ACT)	dadX →	alanine racemase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,239,343	0	C	T	98.5%	212.6 / NA	65	T155T (ACC→ACT)	dadX	alanine racemase 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (34/30);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGC  >  NC_000913/1239098‑1239586                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ggACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGAt                                                                                                                                                                                                                                                 >  1:483112/1‑251 (MQ=255) ggACCGATCCTGAGGCGGGAAGGATTTTTCCATGCTCAGGATCTGGTGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGCGGGCGAAATGACTCTGAt                                                                                                                                                                                                                                                 <  2:546353/251‑1 (MQ=255)   aCCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGt                                                                                                                                                                                                                                               <  2:150830/251‑1 (MQ=255)     cGATCCTGATGCTGGAAGGATTTTTCCCTGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCg                                                                                                                                                                                                                                             <  2:140529/251‑1 (MQ=255)       aTCCTGATGCTGGAAGGATTTTTCCATGCTCCGGATCTGGAGATTTATGAC‑CAGCACCGACTGACCAACTGCGTAAACAGCAACTGGCAGCTCAACGCACTGCAACACGCGCGGCTAAAAGCACCGTTGGCTATTTATCTTAAAGTGAAAAGTGGGATGAATCGGTTGGCATTCCAGACCGCTAGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGACAGGACTCTGATGTCGc                                                                                                                                                                                                                                            >  1:526671/1‑250 (MQ=255)            gATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGCTGTCGCATTTTg                                                                                                                                                                                                                                      >  2:150040/1‑251 (MQ=255)                           tttCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATc                                                                                                                                                                                                                       >  2:97849/1‑251 (MQ=255)                             tCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATcc                                                                                                                                                                                                                      <  1:233855/250‑1 (MQ=255)                                      aGGATCTGGAGATTTAT‑ACACAGCACCGCCTGACCCCCTGCGTACCCAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGTACAGAGGGATGAATCGGTTGGGCTTCCAGCCCGATAGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCCTTTTGCCGAAGCGGAACATCCTGATGGaaat                                                                                                                                                                                                            >  2:528295/1‑249 (MQ=255)                                                   ttATGAC‑CAGCACCGCCGGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATg                                                                                                                                                                                                <  2:205680/250‑1 (MQ=255)                                                      tGAC‑CACCACCGCCTGATCACCAGCGTACACCGCAACTGGCAGCCCAAAGCACTGCAAGATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCTATGGCGAATGTTGGCGAAATGACCCTGATGACGCATTTTGCCGAAGCGGAACATCCTGGTGGAATTTCCGTCGCGTTGgccc                                                                                                                                                                                            >  2:501209/1‑249 (MQ=255)                                                               aCCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAg                                                                                                                                                                                    <  2:109551/251‑1 (MQ=255)                                                                  gCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACCGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGATTGTTGGCGAAATGACTCTGATGTCGCATTTTGTCGAAGCGGAACAGCCTGATGGAATTTCCGGTTCGATGGCGCGTATTGAGCAg                                                                                                                                                                                 <  2:290555/251‑1 (MQ=255)                                                                       accaccTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCGGAGGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAggcggc                                                                                                                                                                            <  1:557624/251‑1 (MQ=255)                                                                                  acgcaGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTGCCAGCCCGATCGCGTGCTTACCGGCTGGCAGCAGTGGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTGGCTGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTg                                                                                                                                                                  <  1:521874/247‑1 (MQ=255)                                                                                  acacAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTgg                                                                                                                                                                 >  1:32454/1‑251 (MQ=255)                                                                                  acacAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTg                                                                                                                                                                  <  2:218109/250‑1 (MQ=255)                                                                                     caGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAgt                                                                                                                                                              <  2:307797/251‑1 (MQ=255)                                                                                            tGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGACTGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcg                                                                                                                                                       >  1:242261/1‑251 (MQ=255)                                                                                               cAGCTCAAAGCACCGCTAACTGAGAGGCTATAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCTGCACGCGCGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGGATGTTGGCGAAATGACTCTGATGTCGCCTTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcggtc                                                                                                                                                    >  2:60736/1‑248 (MQ=255)                                                                                                  cTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTTAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcgttcgtt                                                                                                                                                 >  2:512268/1‑251 (MQ=255)                                                                                                  cTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcgttcgtt                                                                                                                                                 >  2:125565/1‑251 (MQ=255)                                                                                                         cACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGTCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTGGTTGTCCAAt                                                                                                                                          >  2:363517/1‑251 (MQ=255)                                                                                                          aCTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAAtt                                                                                                                                         <  1:123773/251‑1 (MQ=255)                                                                                                               aaaaTGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCACGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTcggcg                                                                                                                                    <  1:186744/251‑1 (MQ=255)                                                                                                                     cgcgGCTAAAAGCACCGTTGGATAGTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGAct                                                                                                                              <  1:22189/251‑1 (MQ=255)                                                                                                                            aaaaGCACCGTTGGATATTTATCTTAAAGTGAACCGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGACGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGGGTCGGCGGGCGTTGTGCAATTGGGCGGCGACTCTGTggg                                                                                                                       >  2:112789/1‑250 (MQ=255)                                                                                                                            aaaaGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGc                                                                                                                       >  1:563569/1‑251 (MQ=255)                                                                                                                                       tGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAgcg                                                                                                            >  1:511711/1‑251 (MQ=255)                                                                                                                                            atTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCAtt                                                                                                        <  1:150040/250‑1 (MQ=255)                                                                                                                                                    tAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTCCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTggg                                                                                               >  1:36820/1‑251 (MQ=255)                                                                                                                                                     aaaGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGCCTCTGTGGCACCCGGAAGCGCATTTTGACTGGGt                                                                                              >  2:224089/1‑251 (MQ=255)                                                                                                                                                     aaaGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGt                                                                                              >  1:364579/1‑251 (MQ=255)                                                                                                                                                       aGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTc                                                                                            <  1:129191/251‑1 (MQ=255)                                                                                                                                                          gAACAGTGGGATGAATCGGTTGGGCTTCCCGCCCTATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATGTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGc                                                                                         <  2:39765/251‑1 (MQ=255)                                                                                                                                                                       aaTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTAtt‑tt                                                                           >  1:150631/1‑251 (MQ=255)                                                                                                                                                                        aTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTg                                                                          <  1:97849/251‑1 (MQ=255)                                                                                                                                                                        aTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTg                                                                          >  2:449512/1‑251 (MQ=255)                                                                                                                                                                           ggTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTg                                                                                      <  1:197657/237‑1 (MQ=255)                                                                                                                                                                           ggTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTg                                                                                      >  2:197657/1‑237 (MQ=255)                                                                                                                                                                               gggCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGc                                                                                                                                                                    >  1:346230/1‑155 (MQ=255)                                                                                                                                                                               gggCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGc                                                                                                                                                                    <  2:346230/155‑1 (MQ=255)                                                                                                                                                                               gggCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGgcg                                                                   <  2:324099/251‑1 (MQ=255)                                                                                                                                                                                    tCCAGCCCGATCGCGTGCTTACCGACTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGCGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTCTGGCGCTTCg                                                              >  2:52977/1‑251 (MQ=255)                                                                                                                                                                                        gCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGt                                                          >  1:27834/1‑251 (MQ=255)                                                                                                                                                                                             gccgcgTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGt                                                     <  1:418062/249‑1 (MQ=255)                                                                                                                                                                                             aTCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGt                                                     >  2:469495/1‑251 (MQ=255)                                                                                                                                                                                                    gCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTgg                                               <  1:536071/250‑1 (MQ=255)                                                                                                                                                                                                    gCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTgg                                               <  1:172732/250‑1 (MQ=255)                                                                                                                                                                                                     cTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCg                                             >  1:263515/1‑251 (MQ=255)                                                                                                                                                                                                        aCCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTg                                           <  1:403205/250‑1 (MQ=255)                                                                                                                                                                                                               ggCAGCAGTTGCGGGCAATGTCGAATGTCGGCGAAATGATTCTGATGTCGCATTTTGCCGAAGCGGAACATCCCGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCTGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCTGCGGCGACTCTTTGGCGCCCGGATTCGCTTTTTGACGGGGTTCGGTCTGGCATTATTCTTG‑ATGGCGGTTCTCCGTCCGGTCAGTGGCGTGATATc                                      <  2:64594/248‑1 (MQ=255)                                                                                                                                                                                                               ggCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGcc                                   >  1:33922/1‑251 (MQ=255)                                                                                                                                                                                                                     agTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATAcc                             >  1:299230/1‑251 (MQ=255)                                                                                                                                                                                                                      gTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCg                            <  1:495579/251‑1 (MQ=255)                                                                                                                                                                                                                         gCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGCGGCGTGATATCGCCAATACCGGAt                         >  1:311045/1‑251 (MQ=255)                                                                                                                                                                                                                             gCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACg                     >  1:290684/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCgg                >  2:291961/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACAGCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGGCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCgg                <  1:361651/251‑1 (MQ=255)                                                                                                                                                                                                                                     gAATGTTGGCGAAATGACTCTGCTGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCCGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGAACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGAAGGGCGCTTCGCCGGCCGGTCAGTGGCGTGAGATCGCCACTACCGGATTAGGGCCGGtt              >  2:389484/1‑249 (MQ=255)                                                                                                                                                                                                                                             gCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGGGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTa     <  1:486439/251‑1 (MQ=255)                                                                                                                                                                                                                                               gAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTa     >  1:30980/1‑249 (MQ=255)                                                                                                                                                                                                                                                 aaTGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAgt                                                                                                                                                              >  2:134276/1‑95 (MQ=255)                                                                                                                                                                                                                                                 aaTGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAgt                                                                                                                                                              <  1:134276/95‑1 (MQ=255)                                                                                                                                                                                                                                                 aaTGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGc  >  2:480276/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGAC‑CAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATT‑TTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGC  >  NC_000913/1239098‑1239586

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