Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	492,768	C→T	A226V (GCC→GTC)  A226V (GCC→GTC)	dnaX → dnaX →	DNA polymerase III subunit tau DNA polymerase III subunit tau

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	492,768	0	C	T	97.8%	136.2 / NA	46	A226V (GCC→GTC)  A226V (GCC→GTC)	dnaX dnaX	DNA polymerase III subunit tau DNA polymerase III subunit tau
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/20);  minor base G (0/1);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GCCGCCGGAGCACGTTAAGTTTCTGCTGGCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGCCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACTTTCGCCTGCTGCACTTGGC  >  NC_000913/492526‑493000                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         gcctccgGAGCACGTTAAGTTTCTGCTGGTGACGACCGATCCACAGAAATTGCCGGTGACTATGTTGTCACGCTGTCTGCAATTGCAGCTCAAGTCGCTGGATGCCTAGCAAATGCGCCATCTGCTTGAGCTCATCCTCGACGTAGAACCGAGCGCTCTCGTGCCGCGGGCGCGGCAATTGCTGGCCCGCGCCTCTGATGGCAGCCTGCGAGATGCCTTGAGCCTGACCGACAAGGCGATTTGCAGCGGTg                                                                                                                                                                                                                                  <  2:155676/247‑1 (MQ=255)    gccgggacACGTTAAGTTTATGCTGGCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCAGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGATGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGCTGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACg                                                                                                                                                                                                                               >  2:111573/8‑251 (MQ=255)        gAGCAAGTTAAGTTTCTGCTGGCGAAGACAGCTCCACAGACATTGCCGGCGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGCGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCCGCGGTGACGGcc                                                                                                                                                                                                                            >  2:15742/1‑250 (MQ=255)                tAAGTTTCTGCTGGCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCa                                                                                                                                                                                                                    >  2:182322/1‑250 (MQ=255)                 aaGTTTCTGCTGGCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAAc                                                                                                                                                                                                                  >  1:196157/1‑251 (MQ=255)                     ttcggcgggCGACGACTGATCCACAGAAATTGCCGGTGACGATTTTGTCTCGCTGGCTGCAATTTCCTCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACAGCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCTATTGTCAGCGGTGACTGCCAGGTTTCAACCCa                                                                                                                                                                                                               <  2:6581/243‑1 (MQ=255)                            ggCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCAGCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCTGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCa                                                                                                                                                                                                       <  1:232353/251‑1 (MQ=255)                                acgacAGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGCGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGCCGGCCAGGTTTCAACCCAGCCGGGCAGTgg                                                                                                                                                                                                   >  2:248399/1‑250 (MQ=255)                                  gacCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCGGTCGTCAAGTTCATCTCAAGGCGCTGGATGTCGAGCAAATGCGCCAGCCGCTTGAGCACATCCTCAACGAAGAACATAGCCCTCACGAGCCGCGGGCGCTGCAATTGCGTGCACGCGCCGCTGAAGGCAGCCTGCGATATGCCTTAAGTCGGACCGACCAGGCGATTGTCAGCTGTGACGGCAAGGTTTCAACCCATTCGGGCAGTGCg                                                                                                                                                                                                  <  2:163125/250‑1 (MQ=255)                                    ccGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCggtg                                                                                                                                                                                               >  2:86652/1‑248 (MQ=255)                                    ccGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATg                                                                                                                                                                                               >  1:222504/1‑251 (MQ=255)                                     cGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGc                                                                                                                                                                                              >  1:230088/1‑251 (MQ=255)                                                      cGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTgacgac                                                                                                                                                                             >  1:49588/1‑251 (MQ=255)                                                                     cACGCTGTCTGCAATTTCATATCAAGGCGCTGGATGTCGAGAAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGATGCAATTGCTGGCACGCGCCGCTGAAGGCAGCATGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCCGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTggcg                                                                                                                                                              >  2:254094/1‑248 (MQ=255)                                                                             cTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGCCGCTGGTTg                                                                                                                                                       >  2:58324/1‑250 (MQ=255)                                                                               gCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAg                                                                                                                                                    <  1:54203/251‑1 (MQ=255)                                                                                 aaTTCCATCGCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCg                                                                                                                                                  <  2:74036/251‑1 (MQ=255)                                                                                  aTTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCg                                                                                                                                                  >  2:65675/1‑250 (MQ=255)                                                                                    ttCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGCTTGTCAGCGGTGACGGCCAGGGTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATg                                                                                                                                               >  2:198917/1‑251 (MQ=255)                                                                                       aTCTCAAGGCGCTGGATGTCGAGCAACTTCGCCATCAGTTTGAGCACATCCTCAACGAAGAACATACCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGACGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCCGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGGCCGTGCGATGCTGGGGACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGGAGCGATGGTc                                                                                                                                            >  2:30144/1‑251 (MQ=255)                                                                                                      aTGTCGAGCAAATTCGCAATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGGCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCg                                                                                                                               >  1:61018/1‑249 (MQ=255)                                                                                                             gCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTa                                                                                                                       >  1:237644/1‑250 (MQ=255)                                                                                                              cAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGCTCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTac                                                                                                                      >  2:81093/1‑249 (MQ=255)                                                                                                              cAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAAt                                                                                                                     <  1:211378/251‑1 (MQ=255)                                                                                                                 aTTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGCTGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGCGGCCAACGGCGAGCGCGTAATGgc                                                                                                                  >  2:132481/1‑251 (MQ=255)                                                                                                                                  gCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGGCGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCg                                                                                                 <  1:232788/251‑1 (MQ=255)                                                                                                                                    acaTCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCt                                                                                               <  1:130996/251‑1 (MQ=255)                                                                                                                                     caTCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTTCGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCTCGTAATGGCGCTGATTAATGAAGCCGCTg                                                                                              <  1:204489/251‑1 (MQ=255)                                                                                                                                                           gCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGc                                                                        >  1:141454/1‑251 (MQ=255)                                                                                                                                                             tCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGt                                                                      >  1:250004/1‑251 (MQ=255)                                                                                                                                                                cGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGc                                                                   >  1:52929/1‑251 (MQ=255)                                                                                                                                                                  aGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCtg                                                                 <  1:86652/251‑1 (MQ=255)                                                                                                                                                                      gcgGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCtggtgg                                                             <  1:55713/251‑1 (MQ=255)                                                                                                                                                                          gcgctgctagtcccggCACGCGCCGCTGAAGGCGGACTGCGAGATGGCGTAGGGCTGACCGACCAGGTGAGTGTTAGAGGTGACGGCCAGGGTTCAACCGAGGCGGTCGGTGGGATGCTGGGTCCGCTTGACGTCGATCAGGTGCTATCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGGGCGCGTAATGGCGCTGAGTAATGAAGCCGCTGCCCGGGGTATCGATTGGGAAGCGTTGCtggttgacat                                                         <  2:231083/237‑7 (MQ=255)                                                                                                                                                                            gcTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCATGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGc                                                       <  2:141541/251‑1 (MQ=255)                                                                                                                                                                                     gCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGt                                              <  1:36102/251‑1 (MQ=255)                                                                                                                                                                                     gCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCATTGCTGGTGGAAATGCTCGGCCTGt                                              >  1:96786/1‑251 (MQ=255)                                                                                                                                                                                        ggCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGc                                           <  2:223682/251‑1 (MQ=255)                                                                                                                                                                                                  gcTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGc                                 <  1:198917/251‑1 (MQ=255)                                                                                                                                                                                                             ccTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCGGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAAc                      <  2:262400/251‑1 (MQ=255)                                                                                                                                                                                                                gCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCTGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACttt                   <  1:263165/251‑1 (MQ=255)                                                                                                                                                                                                                gCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGGGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTGGAAGCGATGGTCGAGGCCAACGGCGAGCGCTTAATTGCGCTGATTAATGAAGCCGCTGCCCGGGGCAGCGAGTGGGAAGCGTTGCGGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAAtttt                   <  1:113440/251‑1 (MQ=255)                                                                                                                                                                                                                   agaTGCCTTAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACTTTCGc                >  1:217333/1‑251 (MQ=255)                                                                                                                                                                                                                         cgtAAGTCTGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACTTTCGCctgctg          <  1:89835/249‑1 (MQ=255)                                                                                                                                                                                                                                 tGACCGACCAGGCGATTGTCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACTTTCGCCTGCTGCACTTGGc  >  1:104189/1‑251 (MQ=255)                                                                                                                                                                                                                                  gtccggccAGTCGATTGACAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCGGGGTACGCTGTACGACGATCAGGCTCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGGGGTATCGAGTGGGAAGCGTTGCTGGTGTAAAGTCTCGGCCTGTTGCATCTTATGGCGATGTTACAACTTGCCCCTGGTGCACTTGGc  <  2:170034/244‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         GCCGCCGGAGCACGTTAAGTTTCTGCTGGCGACGACCGATCCACAGAAATTGCCGGTGACGATTTTGTCACGCTGTCTGCAATTTCATCTCAAGGCGCTGGATGTCGAGCAAATTCGCCATCAGCTTGAGCACATCCTCAACGAAGAACATATCGCTCACGAGCCGCGGGCGCTGCAATTGCTGGCACGCGCCGCTGAAGGCAGCCTGCGAGATGCCTTAAGTCTGACCGACCAGGCGATTGCCAGCGGTGACGGCCAGGTTTCAACCCAGGCGGTCAGTGCGATGCTGGGTACGCTTGACGACGATCAGGCGCTGTCGCTGGTTGAAGCGATGGTCGAGGCCAACGGCGAGCGCGTAATGGCGCTGATTAATGAAGCCGCTGCCCGTGGTATCGAGTGGGAAGCGTTGCTGGTGGAAATGCTCGGCCTGTTGCATCGTATTGCGATGGTACAACTTTCGCCTGCTGCACTTGGC  >  NC_000913/492526‑493000

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