Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	327,999	G→A	S245S (TCC→TCT)	betB ←	betaine aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	327,999	0	G	A	98.6%	205.6 / NA	71	S245S (TCC→TCT)	betB	betaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (34/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGC  >  NC_000913/327754‑328238                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                gggCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAGTGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGagaaga                                                                                                                                                                                                                                            <  1:366174/251‑1 (MQ=255) gggCCGAAGTTAGGTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAAGGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACGCCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATAGCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGagaaga                                                                                                                                                                                                                                            <  2:182503/251‑1 (MQ=255)     cGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGgcc                                                                                                                                                                                                                                        <  1:581512/251‑1 (MQ=255)      gAAGTTAGTTGGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTGTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGgccg                                                                                                                                                                                                                                       <  1:299773/251‑1 (MQ=255)         gttagttTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGccgccg                                                                                                                                                                                                                                    >  1:504919/1‑251 (MQ=255)                   ggATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTATGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACCAACATGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGCTCGAGATACGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGGCCCTTCTTTCAGCGAAGAGGCCGCCGAGTTAGCCAt                                                                                                                                                                                                                          >  1:57652/1‑251 (MQ=255)                    gATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGCAGACGCGGGGGCCCTTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGCTATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCCGCGGTGATTTACCACCCCGTTCCCTGGGCACTTCTTTCAGAGAAGGGGCCGCCGCGTTAGCCATc                                                                                                                                                                                                                         >  1:392421/1‑251 (MQ=255)                           aaCGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACttttt                                                                                                                                                                                                                  >  1:222475/1‑251 (MQ=255)                                 gCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTTATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGc                                                                                                                                                                                                            >  2:160249/1‑251 (MQ=255)                                         gAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGaca                                                                                                                                                                                                    <  1:113599/251‑1 (MQ=255)                                         gAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGACATTGGTACACACCTGACCGGAGCTGAAGAAGGTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAACGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGGAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGac                                                                                                                                                                                                     >  1:378943/1‑250 (MQ=255)                                          aaTGCGCTCAACGCGATCCATAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGAAACTAATACGCGGGTGCCATTAGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCTAGATCCGCATCATCGAAAACGATCAGCTGTGATTTACCGCCCAGTTCCATGTTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGacac                                                                                                                                                                                                   <  1:45826/251‑1 (MQ=255)                                          aaTGCGCTCAAAGCGCGACAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCAAAAAGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCAGAAGAAGATTGCCATCATAGAGATATCGGAGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCCTGGTCACTTCTTTCAGAGAAGAGGCCGACGAGTTAGCCATCACTTTTTTGCCGCTGGCGacac                                                                                                                                                                                                   <  1:45816/251‑1 (MQ=255)                                                  cAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCgg                                                                                                                                                                                             <  1:281678/249‑1 (MQ=255)                                                     cgcgcgCCAGAATTTTCTGCTCAAATGCGGCTTTGCAGTTCGCCGGAACGAAGACGCGGGTGCCATGGGTACACGCCTTGCGGGAGCTGAAGAAGTTTGCCATCAGGGCGTTATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGGGGTGATTTGCCGCCCAGTTCGATGTTCTCTGCTTTCAGAGAAGAGGCCGCGGAGTTGGCCATCCCTTTTTGGCCACTGGCGACTCCGCCGGTAAAt                                                                                                                                                                                        <  2:301010/251‑1 (MQ=255)                                                        gcgcCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCGGGCGACACCGCCGGTAAATGac                                                                                                                                                                                     <  1:178342/251‑1 (MQ=255)                                                              gAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTg                                                                                                                                                                               >  1:164653/1‑251 (MQ=255)                                                               aaTTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTgg                                                                                                                                                                              >  1:572862/1‑251 (MQ=255)                                                                aTTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGc                                                                                                                                                                             >  1:66243/1‑251 (MQ=255)                                                                       gCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATg                                                                                                                                                                         <  1:160537/248‑1 (MQ=255)                                                                           aaaTGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCCTTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCCTCGAACACAATCAGCGGTGATTTACCGACCAGTTCCATGGTCCATTCTTTCAAAGAAGAGGGCGCCGAGGTAGCCATCACTTTTTTGCCGCTGGCGGCAACGCCGGTAAATGACACTTTGGCAATGCTCGGa                                                                                                                                                                   >  1:255715/1‑250 (MQ=255)                                                                           aaaTGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCCTCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGCAAATGACACTTTTGCCATGCCCGGAt                                                                                                                                                                  >  2:140701/1‑251 (MQ=255)                                                                            aaTGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGGTTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCCATGCCCGGATg                                                                                                                                                                 >  2:35162/1‑251 (MQ=255)                                                                             aTGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGGAAACGATCAGCGGTGATTTAACGCCCAGTTCCATGGTCACTTCTTTAAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGacgc                                                                                                                                                                >  1:456686/1‑248 (MQ=255)                                                                               gCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCGGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCt                                                                                                                                                               <  1:160249/250‑1 (MQ=255)                                                                                cGGCTGTGCTTTTCGCGGGACCGAAGAGGCGGGGCCCATTGGTACTGCCCTGGCCGGGGGTGAAGAAGTTTGCCATCATGGCGATATCGTCGGCGAGATCAAAAGCCGCACCATCGAAACCGTTCAGCGGTAATTTACCGCCCATTTCCATGGTCCCTCCTTTCAGGGAAGAGGCCGCTGAGTTACCCATCCCTTTTTTGCCGCTGGCGCCACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCg                                                                                                                                                             <  2:319525/251‑1 (MQ=255)                                                                                  gCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                                                                                                                          >  2:38450/1‑236 (MQ=255)                                                                                  gCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGACCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                                                                                                                          <  1:38450/236‑1 (MQ=255)                                                                                        cctttCGCCGGACCGAAGACGCGGGTGCCATTGGTACCCACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCAAAAACGACCAGCGGTGATTTCCCGCCCAGTCCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGata                                                                                                                                                     <  2:145507/249‑1 (MQ=255)                                                                                        cATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGata                                                                                                                                                     <  1:35162/251‑1 (MQ=255)                                                                                              gCCGAAACGAAGACGCGGGGGCCATTGGTACACCCCTGACCGGAGCTGAAGAATTTTGCAATCATGGCGATATCGGCGGCGAGACCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTCCCGCCCAGTTCAATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGccc                                                                                                                                               <  2:270348/251‑1 (MQ=255)                                                                                               ccGGACCGAAGCCGCGGTTGCCATTGTTACACACCTGACCGGAGATGAAGAAGTTTGCCACCAGGGCGATATCGGCGTCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCg                                                                                                                                              <  1:101637/251‑1 (MQ=255)                                                                                                cGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGAGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCgg                                                                                                                                             >  1:553342/1‑251 (MQ=255)                                                                                                 ggAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGt                                                                                                                                            <  1:533708/251‑1 (MQ=255)                                                                                                  gAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATc                                                                                                                                                                                                                         <  1:6890/173‑1 (MQ=255)                                                                                                  gAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATc                                                                                                                                                                                                                         >  2:6890/1‑173 (MQ=255)                                                                                                   aaCGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtct                                                                                                                                          >  1:313198/1‑251 (MQ=255)                                                                                                      gAAGACGCGGGTGCCATTGGTACACCCCTGACCGGAGCGGAAGAAGTTTGCCACCATGGCGATATCGGCGGCGAGAGCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTCCCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTGGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcg                                                                                                                                       <  2:423691/251‑1 (MQ=255)                                                                                                             cgGGTGCCATTGGTACACCCCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATAGCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTCCCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCAc                                                                                                                                <  2:313198/251‑1 (MQ=255)                                                                                                              gggTGCCATTGGTACACGCCGGACCGGAGCTGAAGAAGTTTGCCAACATGGCGTAATCGGGGGAGAGATCGAGAGCCGAATCTTCGAAAACGACCAGCGGTGATTTCCCGCCCAGTTCCATGGTCCCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACg                                                                                                                               <  2:446561/251‑1 (MQ=255)                                                                                                                     aTTGGCACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGcc                                                                                                                        >  1:186826/1‑250 (MQ=255)                                                                                                                             acacCTGACAGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGAACGCGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCCGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCCGGCGACACCGCCGGTAAAAGACACTTTGGCAATGCCCGGAGGCCCGGTCAGATATTGCCCGGTCTACGCGCCCACGCCCGGCaaaccgtt                                                                                                                 >  2:43902/1‑244 (MQ=255)                                                                                                                              ccccTGACCGGGGCTAAAGAAGTTTGCCATCATGGCGATATCGGGGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTaa                                                                                                               <  2:512332/249‑1 (MQ=255)                                                                                                                                     ccGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCCCTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGcc                                                                                                        <  1:199520/251‑1 (MQ=255)                                                                                                                                       ggAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGt                                                                                                      >  1:539858/1‑251 (MQ=255)                                                                                                                                           cTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGCGTTAGCCATCACTTTTTTGCCGCTGGCGGCACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGACGTCCgg                                                                                                  >  1:520668/1‑251 (MQ=255)                                                                                                                                             gaagaaGTTTGCCATCATGGCGATATCGGCGGCCAGATCGAGATCCGCAACATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGAAACCGCCGGTAACTGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGCCCGGcc                                                                                                >  2:124568/1‑250 (MQ=255)                                                                                                                                             aaagaagTTGGCCATCATGGAGATAGCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGGGGTGATTTCCCGCCCATTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGc                                                                                                 <  2:504919/249‑1 (MQ=255)                                                                                                                                                  agTTTGCCATCTGGGCGATATCGGCGGCGAGATCGAGACCCGCATCATCAAAAACGACCCGCGGTGATTTACCCCCCAGTCCCATGGTCCCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTGGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCGCGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGccc                                                                                           <  2:377817/251‑1 (MQ=255)                                                                                                                                                               ggCGATATCGGTGTCGAGATGGAGATCCGAATCATCGAAAACGATCAGGGGTGATTAACCGCCCAGTCCCATGGTCACTTATTTCAGAGAAGAGTCCTCGGAGTTAGCCATCCCTTTTTTCCCTCTGGCGCCACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGGCAGATAGTGCCCGGTCTCCGCGCCCTCGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTaaa                                                                              <  1:535886/251‑1 (MQ=255)                                                                                                                                                                           tgcgAGATCGAGATCCGCATCATCGAAAACTATCAGCGGTTATTTCCCGCCCAGTTCCATGGTCACTTCTTTCAAAGAAGAGGCGCCCGAGTTAGCCACCACTTTTTTGCCGCGGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCg                                                                                                                                                             <  1:246236/159‑1 (MQ=255)                                                                                                                                                                           ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGCTTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCg                                                                                                                                                             >  2:246236/1‑160 (MQ=255)                                                                                                                                                                           ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCt                                                                  >  1:461654/1‑251 (MQ=255)                                                                                                                                                                                gatcgagatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCATTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTGTTTCCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAAc                                                             <  2:572862/251‑1 (MQ=255)                                                                                                                                                                                 cgcgagatcCTCATCATCGAAACCGATCAGCGTTGATTTACCGCCCAGTTCCATGTTCACTTCGTTCAGAGAAGAGGCCGCCGAGTTGGCCATCACTTTGTGCCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTGCGCTGTAAATTTCAGCCAGCTTTAAcg                                                            <  2:365341/249‑1 (MQ=255)                                                                                                                                                                                   cgagatcCGCATCATCAAAACCGATCAGCGGTGATTTACCGCCCAGTTCCTTGGTCACTTCTTTCAGAAAAGGGGCCGCCGAGTTAGCCATCACGTTTTTGCCGCTGGCGACCCCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTGAAcgcg                                                          <  2:186826/251‑1 (MQ=255)                                                                                                                                                                                                  cGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACt                                           >  1:258057/1‑251 (MQ=255)                                                                                                                                                                                                    aaaccgatcgggggtgtttgcccGCCCATTTCCATGGTCACTTCTGCCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAAc                                            <  2:52867/227‑1 (MQ=255)                                                                                                                                                                                                    aaaaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAAc                                            >  1:52867/1‑248 (MQ=255)                                                                                                                                                                                                     aaaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCg                                        >  1:550471/1‑251 (MQ=255)                                                                                                                                                                                                     aaaCGATCAGCGGTGATTTAACGCCCAGTTCCATGGTCCCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACCCCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATAGTGCCCGGTCTCCGCGCCCAAGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGATGTAAATTTCAGCCAGCTTCAACGCGGTAAGCTGGGtcccttcg                                        >  2:279928/1‑244 (MQ=255)                                                                                                                                                                                                     aaaCGAGCAGCGGTGATTTCCCGCCCAGTCCCATGGTCCCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCg                                        <  2:461654/251‑1 (MQ=255)                                                                                                                                                                                                       ccgaccagcggtttgttACCGCCCAGTTCCATGGTCACTTCTGTCAGAGAAGGGGCCGCCGAGGTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTGCGCt                                      <  2:423081/236‑1 (MQ=255)                                                                                                                                                                                                         gATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCg                                    >  1:394052/1‑251 (MQ=255)                                                                                                                                                                                                          aTCAGTGGTGATTTACCGCCCATTTCAAGGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCACCACTTTTTTGCCGGTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCgg                                   <  2:453171/251‑1 (MQ=255)                                                                                                                                                                                                            aaggggTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCAATTGCGCCATCACTTTTTTGCCCCGGGCGACACCGCCGGTAAATGACACTGTGGCAATGCCCGGATGCGCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATGTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGt                                  <  2:378943/246‑1 (MQ=255)                                                                                                                                                                                                             aGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGttt                                >  2:179509/1‑251 (MQ=255)                                                                                                                                                                                                                ggTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGCGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGa                              >  2:571440/1‑250 (MQ=255)                                                                                                                                                                                                                      ttACCGCCCAGTTCCATGTTCTCTCCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTCCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCa                        <  2:422264/250‑1 (MQ=255)                                                                                                                                                                                                                              cAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTtgc               >  1:268737/1‑251 (MQ=255)                                                                                                                                                                                                                                gTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTtgcct             >  2:317838/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTgccgcc       >  1:205103/1‑251 (MQ=255)                                                                                                                                                                                                                                           aCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCAATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTCACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAgcgc  >  2:372331/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                GGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGC  >  NC_000913/327754‑328238

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