Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	457,253	G→A	A193T (GCG→ACG)	clpP →	ATP‑dependent Clp protease proteolytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	457,253	0	G	A	87.8%	113.3 / 8.4	49	A193T (GCG→ACG)	clpP	ATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (26/17);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

GATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGCCAGCAT  >  NC_000913/457009‑457494                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  gATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTgg                                                                                                                                                                                                                                             >  1:222917/1‑251 (MQ=255)      gcgcTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATAc                                                                                                                                                                                                                                        <  1:546507/251‑1 (MQ=255)       agcTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGTGCGGCTACCCGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTACAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGCAACGGTGGacta                                                                                                                                                                                                                                         >  2:191/2‑246 (MQ=255)           gtCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATGCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCAGTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCGGAAACGGGGGAATACGGt                                                                                                                                                                                                                                     <  2:246625/248‑1 (MQ=255)                        cAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTcgattcgatt                                                                                                                                                                                                                      <  1:243649/251‑1 (MQ=255)                            ggCAAAAGGTAAACGTTTTTGCCTGCAGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTg                                                                                                                                                                                                                   >  2:92588/1‑250 (MQ=255)                                aaaGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCa                                                                                                                                                                                                              <  1:188603/251‑1 (MQ=255)                                     tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGt                                                                                                                                                                                                          >  1:205585/1‑250 (MQ=255)                                      aaaCGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTACAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTa                                                                                                                                                                                                         >  2:213628/1‑250 (MQ=255)                                        aCGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATATCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAAtt                                                                                                                                                                                                      <  1:375528/251‑1 (MQ=255)                                         cGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTg                                                                                                                                                                                                     <  1:352484/251‑1 (MQ=255)                                                     ccGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGgcg                                                                                                                                                                                         >  1:113811/1‑251 (MQ=255)                                                        aaTTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCATGGCCAGGCGCCCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGAGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCa                                                                                                                                                                                       >  2:180329/1‑250 (MQ=255)                                                        aaTTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGCAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGCTGCCAGAGGCGCa                                                                                                                                                                                       >  2:364950/1‑250 (MQ=255)                                                          ttCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACt                                                                                                                                                                                    >  1:529889/1‑251 (MQ=255)                                                              cgcgTGATGATTCACCAACCGTTGGGCGGCTACAAGGGCCCGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGCTGCCAGCGGCGCAACTGTGc                                                                                                                                                                                >  2:214026/1‑251 (MQ=255)                                                               gcgTGATTATTCACCAACCGTTGGGCGGCTGCCAGTGCCAGGCGACCGATATCGAAAGTCATGCCCGTGAAATGCTGAAAGTGAAATGGCTCATGAATGAACTTATGGCGCTTCATACGTGTCAATCATTAGAACAGATTGAACGTGATATCGAGCGCGATCGCTTCCCTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACTCAGCGTAAGTGATGCCAGAGGCGCAGCTGTGcc                                                                                                                                                                               <  2:3787/251‑1 (MQ=255)                                                                             cAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCa                                                                                                                                                                 <  2:20270/251‑1 (MQ=255)                                                                                              ccAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGt                                                                                                                                                >  1:590223/1‑251 (MQ=255)                                                                                                          gACCGTTATCGATATTCAAGCCCGTGAACTTCTCAAAGTTACAGGGCGCATGACTGAACTTATGGCGCTTCATACGGGTCAATCATTAGTACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCAGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATGGATGCCAGCGGCGCAACTGTGCCGCTATCCTTATCCAGGGCGGCACAACGCTGAAAGCGGCTTg                                                                                                                                      >  2:444051/1‑249 (MQ=255)                                                                                                            ccGATAGCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGGTTGAACGTGATTCCGAGCGCCATCGCTTCCTGTCCGCCCCTGAAACGGTGGAAGACGGTCTGGTGGATTCGATTCTGACCTTTCTTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGcc                                                                                                                                  <  1:602420/251‑1 (MQ=255)                                                                                                                         tCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCAtt                                                                                                                     >  1:591494/1‑251 (MQ=255)                                                                                                                          cATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCCTGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCAttt                                                                                                                    >  2:617487/1‑251 (MQ=255)                                                                                                                                   gAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGcgtcgtcg                                                                                                           <  2:123387/251‑1 (MQ=255)                                                                                                                                           gAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGc                                                                                                      <  1:19211/248‑1 (MQ=255)                                                                                                                                           gAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCGGAAACGGTGGAATACGGTCTGGTCGATGCGAGTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGc                                                                                                      <  1:31323/248‑1 (MQ=255)                                                                                                                                           gAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAAGCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTAGACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATTGCATTTGCGTCGTCGTGTGCgg                                                                                                    <  2:64126/250‑1 (MQ=255)                                                                                                                                                ttAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGCATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCCGGGCGACACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGCCGTGTGCggaacca                                                                                               >  1:307540/1‑245 (MQ=255)                                                                                                                                                        cgcATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTAGCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCacaaagaacaaag                                                                                      <  1:294838/251‑1 (MQ=255)                                                                                                                                                          cATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCacaaagaacaaaga                                                                                     >  1:117606/1‑250 (MQ=255)                                                                                                                                                              aatgaaCTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGCGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAgg                                                                                <  2:243438/251‑1 (MQ=255)                                                                                                                                                                  aaCTTATGGCGCTTCATACGGGGCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGtttt                                                                            <  2:590223/251‑1 (MQ=255)                                                                                                                                                                          gcgcTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATg                                                                    <  2:24517/251‑1 (MQ=255)                                                                                                                                                                          gcgcTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATg                                                                    <  1:10065/251‑1 (MQ=255)                                                                                                                                                                              ttCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGAc                                                                  >  1:483002/1‑249 (MQ=255)                                                                                                                                                                                         cAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCaa                                                      >  2:34267/1‑250 (MQ=255)                                                                                                                                                                                                           gAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAAtt                                    <  2:104559/250‑1 (MQ=255)                                                                                                                                                                                                              cGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTg                                >  1:366080/1‑251 (MQ=255)                                                                                                                                                                                                                       gAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGctct                       >  1:122342/1‑251 (MQ=255)                                                                                                                                                                                                                       gAGCGCGATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGctct                       >  1:435639/1‑251 (MQ=255)                                                                                                                                                                                                                           gcgATCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCttttt                   >  1:212297/1‑251 (MQ=255)                                                                                                                                                                                                                              aTCGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGCCCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCCCAACGCTGTAAGCGGCTTGCGCCTGAGACTGGCCTTTACGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGTGTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATAGCTGAATTGCTCtttttacg                >  1:331260/1‑248 (MQ=255)                                                                                                                                                                                                                               tCGCTTCCTTTCCGCCCCTGCAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCgg               >  1:439134/1‑251 (MQ=255)                                                                                                                                                                                                                                cGCTTCCTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGc              >  2:210713/1‑251 (MQ=255)                                                                                                                                                                                                                                     ccTTTCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCaaaa          <  2:352875/250‑1 (MQ=255)                                                                                                                                                                                                                                        ttCCGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGAGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGccc      >  1:77697/1‑250 (MQ=255)                                                                                                                                                                                                                                          ccGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGCTTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGCCAGc    >  2:236776/1‑251 (MQ=255)                                                                                                                                                                                                                                          ccGCCCCTGAAACGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGCCAGc    >  1:36810/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCCCCTGAAGCGGTGGCATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGCCAGCAt  >  1:435500/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  GATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGATATCGAAATTCATGCCCGTGAAATTCTGAAAGTTAAAGGGCGCATGAATGAACTTATGGCGCTTCATACGGGTCAATCATTAGAACAGATTGAACGTGATACCGAGCGCGATCGCTTCCTTTCCGCCCCTGAAGCGGTGGAATACGGTCTGGTCGATTCGATTCTGACCCATCGTAATTGATGCCAGAGGCGCAACTGTGCCGCTATACTTATCCAGGGCGGCACAACGCTGTAAGCGGCTTGCGCCTGAGAATGGCATTTGCGTCGTCGTGTGCGGCACAAAGAACAAAGAAGAGGTTTTGACCCATGACAGATAAACGCAAAGATGGCTCAGGCAAATTGCTGTATTGCTCTTTTTGCGGCAAAAGCCAGCAT  >  NC_000913/457009‑457494

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