Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	891,909	C→T	intergenic (+3/‑58)	nfsA → / → rimK	NADPH‑dependent nitro/quinone reductase NfsA/ribosomal protein S6 modification protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	891,909	0	C	T	98.4%	195.0 / NA	63	intergenic (+3/‑58)	nfsA/rimK	NADPH‑dependent nitro/quinone reductase NfsA/ribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/29);  minor base G (0/1);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGATAATCCGGATCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAACCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTGACGCAGTGA  >  NC_000913/891664‑892150                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  gCGGATAATCCGGATCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCa                                                                                                                                                                                                                                              <  1:193119/251‑1 (MQ=255) gCGGATAATCCGGATCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGAAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGcc                                                                                                                                                                                                                                               >  2:399870/1‑250 (MQ=255)     aTAATCCGGATCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTc                                                                                                                                                                                                                                          <  1:331725/251‑1 (MQ=255)      tAATCCGGATCTTAAGCCGCGCTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCg                                                                                                                                                                                                                                         >  2:434348/1‑251 (MQ=255)            gtaTCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTAt                                                                                                                                                                                                                                   <  2:176904/249‑1 (MQ=255)            ggctCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGAATGAAAACAGCTATCAACCGCAGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATgcat                                                                                                                                                                                                                                   >  2:131351/4‑248 (MQ=255)               tCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGAt                                                                                                                                                                                                                                >  2:178965/1‑251 (MQ=255)               tCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGATGGGCGACGCGCTAAAATCGCCACGTCGATGTATGAt                                                                                                                                                                                                                                >  1:465170/1‑251 (MQ=255)                  tAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAGCCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTAGTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCGCGTCGATGTa                                                                                                                                                                                                                                    <  2:226805/244‑1 (MQ=255)                  tAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGa                                                                                                                                                                                                                                 <  2:15046/247‑1 (MQ=255)                           tttACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCtttt                                                                                                                                                                                                                    >  2:342891/1‑251 (MQ=255)                              cccGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATGTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGAc                                                                                                                                                                                                                 <  2:563022/250‑1 (MQ=255)                              aCCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGa                                                                                                                                                                                                                  <  2:468964/250‑1 (MQ=255)                                              gagCATGACAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGACTATTACCTCACACGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCTACGTCGATGTATGACAAGCGGGCTTTTGACCAGGTCTGACagag                                                                                                                                                                                                   >  2:26103/3‑249 (MQ=255)                                                    gAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCACCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGc                                                                                                                                                                                             >  1:6847/1‑249 (MQ=255)                                                          aGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGa                                                                                                                                                                                      >  1:456681/1‑250 (MQ=255)                                                             tATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAAt                                                                                                                                                                                  >  2:135411/1‑251 (MQ=255)                                                                cAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAAt                                                                                                                                                                                  <  2:181297/248‑1 (MQ=255)                                                                cAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAAt                                                                                                                                                                                  >  1:181297/1‑248 (MQ=255)                                                                    cGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCat                                                                                                                                                                            >  2:632250/1‑250 (MQ=255)                                                                        ggATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTg                                                                                                                                                                        <  2:383996/250‑1 (MQ=255)                                                                         gATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTc                                                                                                                                                                      <  1:180965/251‑1 (MQ=255)                                                                                           cAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATGGCCATATTGTCCCGGGAGGGAACGCTCTa                                                                                                                                                    <  2:79324/251‑1 (MQ=255)                                                                                             gTANGACGAGCAACTGGCGGAATATGACCTCACCCGTTGCATAAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCAGTAAAGAAAGCCGCCCATGGATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGGCTGACAGAGATGTGCAGGTTGAAAATTGCCATATAGTCCCGGTTTGGAACGctctgtt                                                                                                                                                  <  2:594415/251‑4 (MQ=255)                                                                                              tATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCAGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTc                                                                                                                                                 >  1:637259/1‑251 (MQ=255)                                                                                                             gCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCg                                                                                                                                  <  2:7616/251‑1 (MQ=255)                                                                                                              cGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGt                                                                                                                                 >  1:338703/1‑251 (MQ=255)                                                                                                                gAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGa                                                                                                                               <  2:456681/251‑1 (MQ=255)                                                                                                                gAATATTACCTCACCAGTGGCAGCAATAACCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATTTATGATACGCGGTCTTTTGACCAGGTCTGACAGAGAGGTGCATGGTGAAAATTGCCATATTGGCCCGGGATGGAACGCTCTAGTCGTTTAAGCGGCTGCGTGa                                                                                                                               <  2:661477/251‑1 (MQ=255)                                                                                                                                     agcaATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTc                                                                                                           >  2:23021/1‑250 (MQ=255)                                                                                                                                      gcaATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCAc                                                                                                         <  2:480352/251‑1 (MQ=255)                                                                                                                                         ataatCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCTCCTg                                                                                                      <  1:178965/251‑1 (MQ=255)                                                                                                                                         ataatCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGATGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGAGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTg                                                                                                      >  2:106079/1‑251 (MQ=255)                                                                                                                                              cGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGAGCGGTCACCTGGTTGa                                                                                                 <  2:129602/251‑1 (MQ=255)                                                                                                                                                     aTACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCtt                                                                                          <  1:249687/251‑1 (MQ=255)                                                                                                                                                                         cgccgAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACAt                                                                      <  1:575897/251‑1 (MQ=255)                                                                                                                                                                         cgccgAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACAt                                                                      >  1:20160/1‑251 (MQ=255)                                                                                                                                                                               aCAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACCGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACAt                                                                >  2:174203/1‑251 (MQ=255)                                                                                                                                                                                 aaTCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGAGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGATCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATa                                                               >  1:627519/1‑250 (MQ=255)                                                                                                                                                                                          agaaagCCGCCCATTTAGTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTgcgg                                                     <  1:146116/251‑1 (MQ=255)                                                                                                                                                                                               gccgccCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGtct                                                >  1:209279/1‑251 (MQ=255)                                                                                                                                                                                                ccgccCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGtctt                                               <  2:615288/251‑1 (MQ=255)                                                                                                                                                                                                ccgccCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGtctt                                               <  1:174203/251‑1 (MQ=255)                                                                                                                                                                                                      ctttATTCTGGATTATTTGCACAAACAGGGTTGCGCGCCGCGCTACAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTc                                         >  1:217699/2‑251 (MQ=255)                                                                                                                                                                                                          aTTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTa                                     <  2:150186/251‑1 (MQ=255)                                                                                                                                                                                                            tCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTAc                                    >  1:397191/1‑250 (MQ=255)                                                                                                                                                                                                            tCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGACACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATACTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTAc                                    >  1:295635/1‑250 (MQ=255)                                                                                                                                                                                                              tGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACaa                                  >  2:345959/1‑250 (MQ=255)                                                                                                                                                                                                               ggATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACaa                                  >  2:445532/1‑249 (MQ=255)                                                                                                                                                                                                                      ttGCACAATCAGGGTTGTGCGACGCGCTAAAAGCGCCCCGGAGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTTCAGGGTGAAAATTGCCATTTTGTCCCGGGATGGAGCGGTCTATTCGTGTAAGTGTCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTTAGCCTCTTTCTTTCTGTATGAACATAAATCCTCCGGCGTCGGCTATTctttatataggcctca                          <  2:449561/250‑16 (MQ=255)                                                                                                                                                                                                                        gCACAAACAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGt                       >  2:567420/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaCAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTAcc                   >  2:478414/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaCAGGGGTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTAcc                   <  2:397191/251‑1 (MQ=255)                                                                                                                                                                                                                              aCAGGGTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCTTCTTCTATTCACTACAAAGTCCGCAAGTTAcccc                 <  1:581213/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTTGGGCGACGCGCTAAACTCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCAtttt            >  1:198633/1‑251 (MQ=255)                                                                                                                                                                                                                                   gTTGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCAtttt            >  1:517508/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGGGCGACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTg           >  1:516177/1‑250 (MQ=255)                                                                                                                                                                                                                                          gACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTGACGc       >  2:176411/1‑249 (MQ=255)                                                                                                                                                                                                                                          gACGCGCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTGACGCAg     <  1:305783/251‑1 (MQ=255)                                                                                                                                                                                                                                            cgcgCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGt                                                                                                            <  1:513460/146‑1 (MQ=255)                                                                                                                                                                                                                                            cgcgCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGt                                                                                                            >  2:513460/1‑146 (MQ=255)                                                                                                                                                                                                                                            cgcgCTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTGACGCAGTg   >  2:96291/1‑251 (MQ=255)                                                                                                                                                                                                                                             gcgcTAAAATCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACACCATTTTGACGCAGTGa  <  2:141721/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GCGGATAATCCGGATCTTAAGCCGCGTTTACCGGCCTCCATTTTGGTGCATGAAAACAGCTATCAACCGCTGGATAAAGGCGCACTGGCGCAGTATGACGAGCAACTGGCGGAATATTACCTCACCCGTGGCAGCAATAATCGCCGGGATACCTGGAGCGATCATATCCGCCGAACAATCATTAAAGAAAGCCGCCCATTTATTCTGGATTATTTGCACAAACAGGGTTGGGCGACGCGCTAAAACCGCCACGTCGATGTATGATACGCGGGCTTTTGACCAGGTCTGACAGAGAGGTGCAGGGTGAAAATTGCCATATTGTCCCGGGATGGAACGCTCTATTCGTGTAAGCGGCTGCGTGAAGCCGCTATACAGCGCGGTCACCTGGTTGAAATTCTTGATCCGCTTTCTTGCTACATGAACATAAATCCTGCGGCGTCTTCTATTCACTACAAAGGCCGCAAGTTACCCCATTTTGACGCAGTGA  >  NC_000913/891664‑892150

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