Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,941	C→T	T58I (ACT→ATT)	yahH →	putative uncharacterized protein YahH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,941	0	C	T	100.0%	103.9 / NA	38	T58I (ACT→ATT)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/27);  total (11/27)

GCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGAT  >  NC_000913/339680‑340232                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     gCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTATGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAg                                                                                                                                                                                                                                                                <  1:682084/300‑1 (MQ=255)      cactttttGCCTGCTTCTAGGAAGCCATTCTTGGCTAGTGCGAACGTTACGTCGTCTGACCTACATGTTCATGCCGGATGCGGCTTGAACGCCTTATCCTGCCTACGAATTGCGCAAGAATCTGTAGGCCTGATAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCTGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGAGTGgcg                                                                                                                                                                                                                                                           <  2:297603/296‑1 (MQ=17)      cACCTTTTGCCTGCTTCGAAAAAGCCATTCTTGGCTGGTTCAAACGTTACGACTTCCGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGATTGgcg                                                                                                                                                                                                                                                           <  2:322055/300‑1 (MQ=17)                  gCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAg                                                                                                                                                                                                                                              <  1:228063/300‑1 (MQ=255)                          aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATa                                                                                                                                                                                                                                      <  1:304860/300‑1 (MQ=255)                            aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCGGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAg                                                                                                                                                                                                                                    <  2:729740/300‑1 (MQ=255)                                    ttGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                            <  1:741902/300‑1 (MQ=255)                                    ttGGATGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                            <  1:369507/300‑1 (MQ=255)                                     tGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCa                                                                                                                                                                                                                           >  1:688229/1‑300 (MQ=255)                                          ggTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                                                                                                                                                      >  1:637458/1‑300 (MQ=255)                                           ttGCGAACGTTACGTCGTCGGACCTACATGTTCATGCCTGATGCGCCGTGAACGCCTTATCCTTCCTACAAACGGCACAAGAATCTTTATGCCTGATAAGCGTAGCGCATCAGGCAGTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGACGCAAGANTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGANT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                                                                                                                                                     <  2:589111/299‑1 (MQ=14)                                           gTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                                                                                                                                                     <  1:123634/300‑1 (MQ=255)                                                 aCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTATAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGCGTGGCGCGAGAAT‑TGTAAGTCTGATAAGCGTAGCGCATTAGGCagtctg                                                                                                                                                                                                               >  2:296772/1‑295 (MQ=14)                                                     tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                           <  1:379179/300‑10 (MQ=255)                                                     tACGGCGTCTGACCTACATGTTCAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTATGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                           <  1:313317/300‑10 (MQ=11)                                                     tACGGCGTCTGACCTACATGTGCAGGCCGGATGCGGCATGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGTCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                           <  2:139551/300‑10 (MQ=11)                                                     tACGGCGTCTGACCTACAGGATCATGCTGGATGAGGCGTGAAAGCCTTATCCGGNCTGCGAATGGCGCAAGAATCTGTAGGCGTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGGTGCGGCGTGAACTCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                           <  2:601370/300‑10 (MQ=11)                                                     tACGACGTCTGACCGGCATGTTCATGCCGGATGCGGAGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCCGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                           <  2:121455/300‑10 (MQ=11)                                                         gCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgttggt                                                                                                                                                                                                       >  1:699850/1‑287 (MQ=11)                                                                                                                                                            gacctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                                                                     <  1:240260/286‑1 (MQ=11)                                                                                                                                                                tacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCa                                                                                                 >  1:294072/11‑300 (MQ=11)                                                                                                                                                                    tgttcacgCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAACGGCGCAAGAATTTGTAGGCTGGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTNCGAATTGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                                             <  2:200145/293‑1 (MQ=11)                                                                                                                                                                    tgctcatGCCGGATGCGGCGTGGACGCCTTATCCGGCCTATGAATGGCGCAAGGATCCGTAGGCCTGATAAGCGTAGCGCATCAGGAATTTGTCACCATTGCCGGAGGCGGCGTGAACGCCTTATCTTGCCTACGAATGGCGCAAGAATCTGTAGGCCTAATGAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                                             <  2:332965/294‑1 (MQ=18)                                                                                                                                                                      ttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAccc                                                                                           <  1:382351/296‑1 (MQ=255)                                                                                                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                                      <  1:671204/300‑1 (MQ=255)                                                                                                                                                                                agacggcgTGAACGCCTTATCCTGCCTACTAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                 <  2:169336/297‑1 (MQ=14)                                                                                                                                                                                aTGCGGCGTGAACGCTCTATTCGACCTACGAATGGCGCAAGAAGCTTTAGGCTTGATAAGCGTGGCGCAGCAGGCGTTGGTCACTATTACCGGATGCGGCGTGAACGCCTTGTCTGGCCGACGAACGGCGCAAAAATCTGTAGGCCTGATAAGCTTAGCGCATCAGGCATTTGTCNCCATTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTAATAAGCGTAGCGCATCGGGCAGTCTTGCNTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                 <  2:581446/300‑1 (MQ=14)                                                                                                                                                                                aTGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                 <  1:510309/300‑1 (MQ=255)                                                                                                                                                                                aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                 <  1:66826/300‑1 (MQ=255)                                                                                                                                                                                         gAACGCCTTATCCGGCCTACGAATGGCGCAATAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑Taacaa                                                                       <  1:608247/300‑1 (MQ=255)                                                                                                                                                                                             ggattATCCGGACTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGAGCACCAGGCATTTGTCAACGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGGATCTGCAGGCCTGATAGGCGTAGCGCATCAGGCATTTGTCGCCATTACCGGACGCGGCGTGACCGCCTTCTCAGGCCTCCGAGTGGCGCAAGGATTGTAGGCCTGATGAGCGTGGCGCATCCGGCAGTCTGGCGTTGGTCATAACCCAATCACCTTCTGTAggtgaaccgagcaaatg                                                                    >  2:157507/4‑284 (MQ=9)                                                                                                                                                                                                 tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAAAATGaacg                                                               >  1:90103/1‑298 (MQ=255)                                                                                                                                                                                                       gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTtg                                                         >  1:452816/1‑300 (MQ=255)                                                                                                                                                                                                           aCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTTTAGCGGACA‑TAACAACATGAAAGAGCTTGTGGt                                                     <  1:810226/300‑1 (MQ=255)                                                                                                                                                                                                                    gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGTGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAACTTGTGGTCGTTGCCAt                                            >  2:461573/1‑300 (MQ=255)                                                                                                                                                                                                                                   ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCAt                             <  1:612201/300‑1 (MQ=255)                                                                                                                                                                                                                                      gATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTCGTGGTCGTTGCCATTGGTGGCAACAGCCTTAt                          >  1:46536/1‑300 (MQ=255)                                                                                                                                                                                                                                                  gcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCCTAACCCCATCACCCTCTGTAGCGG‑CACTAACCACATGCAAGAACTTGTGATCGTTGCCATTGGTGGCAACAGCATTATCACAGATAACgc              >  2:295485/1‑300 (MQ=17)                                                                                                                                                                                                                                                            tttatcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGAAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATANGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCCCCCTCTGTAGCGGACA‑TAACAACACGAAAGAGCTTATGGTCGTCGCCGTTGGTAGCAACAGCATTACCAAAGATCACGACAGCCAGtag    >  2:206373/6‑298 (MQ=18)                                                                                                                                                                                                                                                              tgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGAt  <  2:330832/297‑1 (MQ=255)                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     GCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACA‑TAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGAT  >  NC_000913/339680‑340232

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: ‑