Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	851,631	A→G	Y323H (TAC→CAC)	opgE ←	phosphoethanolamine transferase OpgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	851,631	0	A	G	98.5%	249.3 / NA	68	Y323H (TAC→CAC)	opgE	phosphoethanolamine transferase OpgE
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/30);  total (38/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

GCTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGG  >  NC_000913/851388‑851873                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   gCTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCa                                                                                                                                                                                                                                             <  1:52168/251‑1 (MQ=255)  cTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGCGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTCTCTACTTCAAACTTACCTTCTTTCAGGAAGTGCACTTCAt                                                                                                                                                                                                                                            >  1:270247/1‑251 (MQ=255)    tGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCAt                                                                                                                                                                                                                                            <  1:13476/249‑1 (MQ=255)        cAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATCGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGAGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTCTCTGCTTCACAATTACCTTCTTTCAGGCAGTGCACTTCCTCTGCTc                                                                                                                                                                                                                                      >  2:165741/1‑251 (MQ=255)             aaTCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTg                                                                                                                                                                                                                                 <  2:111262/251‑1 (MQ=255)                       ttGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTg                                                                                                                                                                                                                        <  2:479159/250‑1 (MQ=255)                        tGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGc                                                                                                                                                                                                                      >  1:416458/1‑251 (MQ=255)                              atgttataGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGc                                                                                                                                                                                                                <  2:79114/247‑1 (MQ=255)                                     agagATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCg                                                                                                                                                                                                         >  1:2582/1‑251 (MQ=255)                                           aGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCCTATTc                                                                                                                                                                                                    <  2:481407/250‑1 (MQ=255)                                                 acgTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAtt                                                                                                                                                                                             <  1:409856/251‑1 (MQ=255)                                                  cgTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAttt                                                                                                                                                                                            <  1:411466/251‑1 (MQ=255)                                                   gTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCGGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTg                                                                                                                                                                                           <  1:255248/251‑1 (MQ=255)                                                            gATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGtt                                                                                                                                                                                  <  2:452368/251‑1 (MQ=255)                                                              ttGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTgg                                                                                                                                                                                >  2:132624/1‑251 (MQ=255)                                                                 cacaAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCGGGTTGGaaa                                                                                                                                                                             <  1:132624/251‑1 (MQ=255)                                                                  acaAAGGCTTCGTATTTTCCTTGTGTCCTGACGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGATGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGaaa                                                                                                                                                                             >  1:77956/1‑250 (MQ=255)                                                                   caAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAac                                                                                                                                                                           >  2:373126/1‑251 (MQ=255)                                                                    aaaGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATGCGCCGATTTGACCCTGGTTGGAAAacc                                                                                                                                                                          <  2:406344/251‑1 (MQ=255)                                                                     aaGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAacc                                                                                                                                                                          >  1:565089/1‑250 (MQ=255)                                                                             gTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCt                                                                                                                                                                 <  1:519582/251‑1 (MQ=255)                                                                               aTTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTgg                                                                                                                                                               >  2:284604/1‑251 (MQ=255)                                                                                 tttCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGa                                                                                                                                                              >  1:532666/1‑250 (MQ=255)                                                                                  ttCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGCGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGCTTTGACCCTGGTTGGAAAACCACCAGGTCTGGa                                                                                                                                                              >  1:591427/1‑249 (MQ=255)                                                                                   tCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGCTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGc                                                                                                                                                           >  2:178949/1‑251 (MQ=255)                                                                                       tgtgTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCcgc                                                                                                                                                       >  2:493431/1‑251 (MQ=255)                                                                                              tGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAttt                                                                                                                                                >  1:332279/1‑251 (MQ=255)                                                                                               gTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAtttc                                                                                                                                               >  1:450597/1‑250 (MQ=255)                                                                                                   cAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGACATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGAAAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGcc                                                                                                                                             >  2:508169/1‑249 (MQ=255)                                                                                                             gATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAg                                                                                                                                 <  1:541452/251‑1 (MQ=255)                                                                                                             gATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTGGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGAGACAAAg                                                                                                                                 <  2:92003/251‑1 (MQ=255)                                                                                                                 tgAGCCCATCAGATGTAGCACAATCAGCTGCGGTTGCGCGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTAt                                                                                                                             >  2:181829/1‑251 (MQ=255)                                                                                                                  gAGCCCAACAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGGGGGAAAACCACCAGGTCTGGAATGCCGCGCGATTTGCCAGGGTGACAAAGTTAt                                                                                                                             <  2:93679/250‑1 (MQ=255)                                                                                                                   aGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAttt                                                                                                                                                                                            >  2:60520/1‑186 (MQ=255)                                                                                                                   aGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAttt                                                                                                                                                                                            <  1:60520/186‑1 (MQ=255)                                                                                                                                 aGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTAGCCTGAAACTGTGGtt                                                                                                              <  2:216009/250‑1 (MQ=255)                                                                                                                                       aTCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATc                                                                                                       >  2:136208/1‑251 (MQ=255)                                                                                                                                          aGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGAc                                                                                                    <  2:337398/251‑1 (MQ=255)                                                                                                                                              gCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACaa                                                                                                  <  2:137586/249‑1 (MQ=255)                                                                                                                                                ggTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAAcgcg                                                                                              >  2:426957/1‑251 (MQ=255)                                                                                                                                                 gTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAAcgcg                                                                                              >  2:270738/1‑250 (MQ=255)                                                                                                                                                    gCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTTCCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGTTTTGCCATCGACAACGCGATTg                                                                                          <  1:555067/251‑1 (MQ=255)                                                                                                                                                      gAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGCCCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTg                                                                                          >  1:556288/1‑249 (MQ=255)                                                                                                                                                      gAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTg                                                                                          >  1:145291/1‑249 (MQ=255)                                                                                                                                                         tGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGt                                                                                     >  2:137248/1‑251 (MQ=255)                                                                                                                                                                                ggTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTgg                                                              <  1:323046/251‑1 (MQ=255)                                                                                                                                                                                 gTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTgg                                                              >  1:255762/1‑250 (MQ=255)                                                                                                                                                                                     tatCCAGTAACGCTTGGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATCGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTCCTGCGTGGAGcc                                                          >  2:105226/1‑250 (MQ=255)                                                                                                                                                                                        ccAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAGATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTg                                                      >  1:256789/1‑251 (MQ=255)                                                                                                                                                                                         cAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGGCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTgg                                                     <  2:576294/251‑1 (MQ=255)                                                                                                                                                                                            tAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGcc                                                          >  1:120319/1‑243 (MQ=255)                                                                                                                                                                                            tAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGcc                                                          <  2:120319/243‑1 (MQ=255)                                                                                                                                                                                             aaCGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTggcgg                                                  <  2:37509/250‑1 (MQ=255)                                                                                                                                                                                             aaCGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGAGTTCTGCGTGGAGCCACTggcgg                                                  <  1:489346/250‑1 (MQ=255)                                                                                                                                                                                               cGCTTCCTCTCTGCTGTATTTATCTGATCCAAAATTACGTACTTCCCGGAAGAACACTTCATCTGCTCGTTTGGCGATGCTGTCGATAGCAGCATCGTATTCGCCGATTTGACCCTGGTTGGAAACCCACTAGCACTGGAAGCCAGCGCGATCTGCAAGGGTCACAAAGTTATCGTGAACCTGTGGTTTGCCATCGACAACGCGTTTGAGCGTTAAGCCAAGCGATTTCTGCGTTGAGCCACTggcggaa                                                >  2:278569/1‑248 (MQ=255)                                                                                                                                                                                                     gTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGt                                         >  1:87865/1‑251 (MQ=255)                                                                                                                                                                                                      tCTTTGGTGTTTTTATCTGCTTCAACCTTACCTTCTTTCAGGAAGTGCACTTAATCTGCTCGTTTGGCGATGCTGGCGCTAGCGGTATCGTATTCGCCGATTTGCCCCTGGTTGGACAACCACCAGGTCTGGCAGCACGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCCTCGCCAACGCGCTTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCCCTGTAGTc                                        >  2:277469/1‑251 (MQ=255)                                                                                                                                                                                                       cTTTGGTGATTTTCTCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCCTCTGATCGTTTGGCGATGCTGGCGATCGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGCATCCACTGGCGGCAATGTAgcc                                        >  2:599620/1‑248 (MQ=255)                                                                                                                                                                                                         ttGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAg                                      <  1:270738/250‑1 (MQ=255)                                                                                                                                                                                                                ttttATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGCTTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAatacc                              >  2:368561/1‑249 (MQ=255)                                                                                                                                                                                                                  ttATCTGCTTCAAACTTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGtt                                                                                    >  1:471822/1‑195 (MQ=255)                                                                                                                                                                                                                  ttATCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGtt                                                                                    <  2:471822/195‑1 (MQ=255)                                                                                                                                                                                                                    aTCTGCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGGCAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAAc                          >  2:126947/1‑251 (MQ=255)                                                                                                                                                                                                                        gCTTCAAAATTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGtt                      <  1:493431/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaaTTACCTTCTTTCAGGCAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTCTCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACgc                 >  1:217241/1‑251 (MQ=255)                                                                                                                                                                                                                              aaaTTACCTTCTTTCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACgct                >  1:298665/1‑251 (MQ=255)                                                                                                                                                                                                                                           tCAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGCAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCCAACggg   >  2:330722/1‑251 (MQ=255)                                                                                                                                                                                                                                            cAGGAAGTGCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAGACgggg  <  1:136208/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GCTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGG  >  NC_000913/851388‑851873

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