Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,935	C→T	A56V (GCC→GTC)	yahH →	putative uncharacterized protein YahH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,935	0	C	T	97.4%	101.1 / ‑3.7	39	A56V (GCC→GTC)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/17);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑CCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTC  >  NC_000913/339654‑340229                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       gTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGNGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATgc                                                                                                                                                                                                                                                                                        >  2:630921/1‑298 (MQ=255)    ggCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGt                                                                                                                                                                                                                                                                                   >  2:694053/1‑300 (MQ=255)                gcgGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                                                                                                                                       <  1:762840/300‑1 (MQ=255)                  gggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                                                                                                                                      <  1:27912/299‑1 (MQ=255)                   ggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                                                                                                                                      >  1:504175/1‑298 (MQ=255)                   ggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑tcacc                                                                                                                                                                                                                                                                                                    <  1:204277/268‑1 (MQ=255)                   ggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑tcacc                                                                                                                                                                                                                                                                                                    >  2:204277/1‑268 (MQ=255)                    gCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                                                                                                                                   >  1:562602/1‑300 (MQ=255)                    gCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGAACGCAGCGTGACCGCCCTATCCGGc                                                                                                                                                                                                                                                                   >  1:563294/1‑300 (MQ=18)                            cATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGCTGCGGCGTGACCGCCTTATCCGGCCTACGAGt                                                                                                                                                                                                                                                           >  2:519038/1‑300 (MQ=255)                                     tttGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCgaga                                                                                                                                                                                                                                                  >  1:17057/1‑300 (MQ=255)                                            gCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAttgta                                                                                                                                                                                                                                            >  1:524209/1‑295 (MQ=255)                                                           ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACTATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTCGGCCTGATAGGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGCTGCGGCGTGAACGCCTTATCCGGCATAAGAATGACGCGAGAAT‑TGTAGGCCTGATAAGCGt                                                                                                                                                                                                                             >  2:380007/1‑298 (MQ=17)                                                             cTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATNGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTTATACGTGTCGCGCATCCGGCATCTG‑TCAACCTTGCCAGATGCGGCGTGAACGCCTTATCCGGCATACGATTGGCGCGAGAAT‑TGTAGGCCTGATACACGTAgc                                                                                                                                                                                                                          >  2:593472/1‑299 (MQ=14)                                                              tgtgCTGGTGCGAACGTTACTGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTATGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                        <  2:301242/297‑1 (MQ=255)                                                                   tGGTGCGAACGCTACGGCGGCTGACCTACATGGTAATGCCGGATGCGGCGTGGACGTCTTCTCCTGCCTGCGAACTGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATTTGTAGGCCTGGTAAGCGTAGCGCATCAGGCATTGG‑TCACCATTGCCGGATGCGGCGTGAACGCCCTATCCGGCCTAGGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGcagcag                                                                                                                                                                                                                   <  2:372271/300‑1 (MQ=11)                                                                    ggTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCGAGGATCTGTAGGACTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAGTCTGCAGGCCTGCTAAGCGTAGCGCATCAGGCATTTGTCCAC‑ATTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACGCGTGGCGCGAGAAT‑TGTAGGCTTGATAAgtttagc                                                                                                                                                                                                                          >  2:165916/1‑286 (MQ=14)                                                                    ggTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTTATAAGCGTAGCGCATCa                                                                                                                                                                                                                    >  2:654147/1‑298 (MQ=255)                                                                               tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                       <  1:635750/300‑10 (MQ=255)                                                                                    cGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCATGCagtctggcgttggtc                                                                                                                                                                                                  >  1:297457/1‑286 (MQ=11)                                                                                                                                   ggCCTACGAATGGCGCAGGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGTGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTa                                                                                                                                                                                                                                                                >  2:304059/1‑192 (MQ=11)                                                                                                                                                                                      gacctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATANGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                                                                 <  1:169452/286‑1 (MQ=11)                                                                                                                                                                                       acctacctgttcattccGGATGCGGCGTGAACGTCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                <  2:725919/285‑1 (MQ=11)                                                                                                                                                                                          tacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                                               <  1:632706/288‑1 (MQ=11)                                                                                                                                                                                          tacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                <  1:457802/287‑1 (MQ=11)                                                                                                                                                                                              tgttcaggCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCGGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                                         <  1:630921/293‑1 (MQ=11)                                                                                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                                  <  1:115802/300‑1 (MQ=255)                                                                                                                                                                                                          aTGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAACCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                             <  2:146846/300‑1 (MQ=255)                                                                                                                                                                                                          aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                             <  2:315158/300‑1 (MQ=14)                                                                                                                                                                                                                 gggAACGCCGTATCCGTCCTACGGATGGCGCAAGAATCTGTAGGCCTGAGAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCTTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                <  2:2418/272‑1 (MQ=11)                                                                                                                                                                                                                        ccTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCATCTGTAGCGGACATAAAAACATGa                                                               >  1:140177/1‑300 (MQ=255)                                                                                                                                                                                                                            aTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTTTAGCGGACATAACAACATGAAag                                                            <  2:763253/299‑1 (MQ=255)                                                                                                                                                                                                                                 gacTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGAATCAGGCAGTTGTCGCCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTAAGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGAGTTGGTCGTGACACCATCAACCTCTGCAGAGGACGTAACAAAATGAAagcgatt                                                       >  2:733067/3‑294 (MQ=11)                                                                                                                                                                                                                                       gAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACNCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGAGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTc                                                 <  1:176104/299‑1 (MQ=255)                                                                                                                                                                                                                                            gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGcc                                           >  1:499865/1‑300 (MQ=255)                                                                                                                                                                                                                                              gcacgAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTT‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTAAGAGTGGCGCGAGAATTGTAGGCCTGATAAGTGTAGCGTATCAGGCAGTCTGGCGTTGCTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGACAGAGCTTTCGGTCGTTGCCAt                                         >  2:333809/5‑300 (MQ=17)                                                                                                                                                                                                                                              gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGAATGGCGCGAGAATTGTAGGCCTGATCAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAACCTCTGTAGCGGACACAACAACATGAAAGAGATTGTGGTCGTTGCCAt                                         >  2:333819/1‑300 (MQ=17)                                                                                                                                                                                                                                                                     gCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACNCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCaaa                   >  1:580455/1‑299 (MQ=255)                                                                                                                                                                                                                                                                     gCGTAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACNCCATCACCCTCTGTAGCGGACATAACAAAATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCaaa                   >  1:580556/1‑299 (MQ=255)                                                                                                                                                                                                                                                                        tAGCGCATCAGGCATTTG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTNTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGAt                <  1:160517/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                      tttG‑TCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTc  >  1:183395/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       GTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTG‑CCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTC  >  NC_000913/339654‑340229

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