Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	353,610	C→G	18.7%	A302G (GCG→GGG)	prpE →	propionate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	353,610	0	C	G	18.7%	142.7 / 9.3	71	A302G (GCG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base G (12/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAC  >  NC_000913/353373‑353835                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     gggTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGc                                                                                                                                                                                                                                      <  2:418146/238‑1 (MQ=255)     tccGGGGGCGTGGCTGGAATCCAACGAAACCGCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCGGCTGGCGGGGAT‑GGCGACt                                                                                                                                                                                                                                  <  2:1230757/237‑1 (MQ=255)          tggcgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGt                                                                                                                                                                                                                           <  2:591367/238‑1 (MQ=255)           ggcgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGG‑GGCGGCCC‑T‑CGttt                                                                                                                                                                                                                         >  2:1171926/1‑239 (MQ=255)             cgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGAATA‑T‑CGGTTAc                                                                                                                                                                                                                       >  1:1217163/1‑239 (MQ=255)             cgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGGCCA‑T‑CGGTTAc                                                                                                                                                                                                                       >  1:127673/1‑239 (MQ=255)                          cAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGa                                                                                                                                                                                                          <  1:28679/239‑1 (MQ=255)                                aaCCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                     <  1:1035599/238‑1 (MQ=255)                                      ctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGGCCATT‑CGTTTGC‑AAGGGTGTCCGACCTGGCCgggcg                                                                                                                                                                                              >  2:268894/1‑236 (MQ=255)                                           ttCTCTCCACTCCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                         <  2:1030778/239‑1 (MQ=255)                                             ctctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACCAAGGATTGCCGACCTGGCCGGACTGCGGCgt                                                                                                                                                                                       >  1:250877/1‑239 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGACTA‑T‑CGTTTAAGCAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      >  1:1053242/1‑239 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGACTA‑T‑CGTTTTCGAAGGATTGCCGACCTGGCCGGACTgcggcggg                                                                                                                                                                                      >  1:727930/1‑237 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑TGCGACTA‑T‑CGTTTTCGGAGGATTGCCGACCTGGCCGGGCTgcggcggg                                                                                                                                                                                      >  1:525130/1‑237 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      <  1:1170925/239‑1 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      <  1:263084/239‑1 (MQ=255)                                                tACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGAc                                                                                                                                                                                                                                   <  2:43658/194‑1 (MQ=255)                                                tACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGAc                                                                                                                                                                                                                                   >  1:43658/1‑194 (MQ=255)                                                      tCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGGCTA‑T‑CGTTTCCGAAGGGTTGGCGGCCTGGCCGGACGTCGGCGTGTGGGGGAc                                                                                                                                                                              >  1:1126449/1‑239 (MQ=255)                                                       ccGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGACTA‑T‑CGTTTACGAAGGGTTTCCCACCCGGCCCGACTGGGGCGTGTGGTGGGCa                                                                                                                                                                             >  1:311622/1‑239 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGGACTA‑‑‑CGGTTTCGGAGGGATGCCCACccggccgggcgggcgggggtggggggaaa                                                                                                                                                                            >  2:263084/1‑218 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGGGCT‑T‑CGTTTACTAAGGATTGCGAACCCGGGCCGACTgcggcggttgtggcggac                                                                                                                                                                            >  2:791545/1‑227 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGGCTA‑T‑CGTTTACGCAGGATTGCCGACCCGGCCGGGCTGCGGCGTGTGGTGGTCaa                                                                                                                                                                            >  1:975555/1‑239 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGG‑GGGGCCTA‑A‑CGGTTACCAAGGAGTGTCGCGCTTGGCGGGCTgcggcggggtgggggaaa                                                                                                                                                                            >  2:329412/1‑227 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGGTTACGAAGGATTGCCGACCTGGGCGGACTGCGGCGTGTGGTGGGCaa                                                                                                                                                                            >  1:942808/1‑239 (MQ=255)                                                          gCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGGGTTCTTTTGCGCATCGGATATCGGCTGGGCGGTGGGGCATTTGTGTACCGTTTAGGCGCCGCTACTGGCGGGGGT‑GGGGGCTT‑TATGTTTGTGATGG‑TTGCCGACCTGGCCGGGAGgcggcggggggggggacatc                                                                                                                                                                          >  2:757896/1‑225 (MQ=255)                                                                ccGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCga                                                                                                                                                                     <  2:489527/238‑1 (MQ=255)                                                                 cGGCAAACCTAAAGGCGGGCAGCGGGACGTCCGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCgaga                                                                                                                                                                   <  2:311622/239‑1 (MQ=255)                                                                         cTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAAtcgt                                                                                                                                                                        <  1:211812/226‑1 (MQ=255)                                                                          tAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAgg                                                                                                                                                          <  2:457148/239‑1 (MQ=255)                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGGGGCTA‑T‑CGTTTACGAAGGGTTGCCGACCTGGCCGGGCTGCGGCGTGTTGGGGGCAATCGTCGAGAAATATCCAGTTAGCTgg                                                                                                                                                  >  1:160641/1‑239 (MQ=255)                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGGCAATCGTCGAGAAATATCAGGTTAGcccg                                                                                                                                                  >  1:101088/1‑237 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGGCTA‑T‑CGTTTTCGTAGGATTGCCGCCCCCGCCGGGCTGGGGGGGGTGGGGGGCAAatgtcggggaaaaatctgtttcttat                                                                                                                                                  >  2:211812/1‑218 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGG‑GGCGGCTA‑A‑CGGTTTCGAAGGATTGCCGACCCGGCCGGGGTTCGGCGGGGGGGGGGCAATCCGCGGGAAATAGCAGGttgaccgcg                                                                                                                                                 >  2:829248/1‑232 (MQ=255)                                                                                         tGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGGGACTA‑T‑CGTTTACGAAGGATTGCCGGCCTGGCCGGGCTGCGGGGGGGGGGGGGCAATTGTCGGGAAATAACAGGTTAGCCGGGTGGtct                                                                                                                                           >  1:45609/1‑239 (MQ=255)                                                                                                               gcgcTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGGCCGCCATTcgctg                                                                                                                     >  1:1040967/1‑237 (MQ=255)                                                                                                               gcgcTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACCA‑T‑CGTTTACGCAGGGTTGCCGCCCTGGCCCGACTGCGGGGTGTGGGGGAGAATCGTGGGGAAATATCTGGTTAGGCGGATGTTCTCAGCGCCGACCGcccattgcgt                                                                                                                     >  2:1142402/1‑231 (MQ=255)                                                                                                                         ccTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGGATA‑T‑CGTTTACGGAGGGTTGCCCGCCTGGCCGGACTGCGGCGGGGGGTGGGCAATCTGCGGGAAATATCAGGGTTGCCCGGTGGTCTCATCGCCAGCCGCCCCTCGCGGGCTGagaaaa                                                                                                           >  2:983912/1‑239 (MQ=255)                                                                                                                               tGGACACCATTTTGGGCGGCAAAGCGGGCAGCGTGTTCTTTGGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                      <  1:599680/238‑1 (MQ=255)                                                                                                                                              gcggcAAGGCGGGCGGCGTGTTCTTTTGCGCATCGGATAGCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                       <  2:727930/238‑1 (MQ=255)                                                                                                                                               cggcAAAGCGGGCGGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGc                                                                                      <  1:2449/238‑1 (MQ=255)                                                                                                                                                ggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTTCGGCGTGTGGGGGACAATCGTCGAGAAATATCAGGGTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGCAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    >  1:901869/1‑239 (MQ=255)                                                                                                                                                ggaaaaGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCa                                                                                     <  1:268894/235‑1 (MQ=255)                                                                                                                                                            aGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGCGACTA‑T‑CGTTTACGAAGGATTTCCGAACTGGGCGGGCTGCGGCGCGTGGGGGGCAATCGTCGGGAAATATCAGGTTAGCCGGGTGTTTTCTGCGCCGGCCGCCATTCGCGTGCTGAAAAAAATCCCTACCGCTGTAATTCGCAAACACGAtctct                                                                         >  2:448096/1‑238 (MQ=255)                                                                                                                                                               gtgtTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                        <  1:993544/236‑1 (MQ=255)                                                                                                                                                                 gtTCTTTTGCGCATCGGATGTCGGCGGGGTGGTGGTGCATTCGTATACCGTTTCCGCGCCGCTGCTGGCGGAGAG‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCGGGCCGGACGGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGc                                                                   <  2:1042108/239‑1 (MQ=255)                                                                                                                                                                      tttGCGCATCGGTTATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                              <  2:1245928/239‑1 (MQ=255)                                                                                                                                                                         gcgcATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTg                                                          <  1:1142402/239‑1 (MQ=255)                                                                                                                                                                                  atagCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATct                                                  <  1:1171926/238‑1 (MQ=255)                                                                                                                                                                                  atagCGGCTGGGTGGGGGGGCATTGGTATATCGTTTACGCGCCGCGGCGGGCGGGGAG‑GGCGACTA‑T‑CGTTTACGAAGGATGGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTg                                                          <  1:866357/230‑1 (MQ=255)                                                                                                                                                                                    atCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATctg                                                 >  1:942857/1‑237 (MQ=255)                                                                                                                                                                                          tGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCGGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGCCCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGgag                                          <  1:580851/238‑1 (MQ=255)                                                                                                                                                                                           gggTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGGGGCGTGTGGGGGACAAACGTCGAGAAATATTAGGGTAACCGGGTGGTCTCAGCGCCGACCGCCATTCGCGTGCCGAAAAAATTTCCCAACCCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGGTGC‑CTATCTGGATGGGGaa                                        >  2:1248934/1‑239 (MQ=255)                                                                                                                                                                                            ggtggtGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGGGGGGTGGGCAAACGTCGAGAAATATCAGGGTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGGGGTGAAAAAATTCCCTACCCCTGAAATTTCCAAAAACGAACTCTCGTCGCTGGGA‑GGGGTCCAActggcctgtttta                                       >  2:866357/1‑231 (MQ=255)                                                                                                                                                                                            ggtggtGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGGCCA‑T‑CGTTTACGAAGGATTGCCGACCCGGCCGGACTGCGGCGTGGGGGGGGCAATCGTCGCGGAATATCAGGGTTGCCGGGTGGTCTCAGCGCCGACCCCCCTTTCCGTGCCGGAAAAATTTCCTACCGCTGAAATTCGCAAACACGGTCTCTCCGCGCTGGGC‑GTGCTCTCTCTGGGTGGACaa                                        >  2:813514/1‑238 (MQ=255)                                                                                                                                                                                                 tggGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGGCCATCCGCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCCGAAATTCGCAAACACGATCTCTCGCCCCGGGGA‑GTGCTCTATCTGGCTGGGGAACCGccg                                  >  2:19622/1‑237 (MQ=255)                                                                                                                                                                                                   gggCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTgg                                 >  1:11297/1‑238 (MQ=255)                                                                                                                                                                                                    ggCATTCGTATATCGTTTACGCGCCGCTGTGGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGa                                <  1:754165/238‑1 (MQ=255)                                                                                                                                                                                                      cATTCGTTTATCGTTTACGCGCCGCTGCTGGGGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACg                              <  1:312168/238‑1 (MQ=255)                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGT‑GGGGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGTGGGCAATCGTCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCGCCGACCGCCCTTTGCGTGCTGAAAAAATTCCCCACCGCTGAAATTCGCCAACACGAACTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAg                            >  2:877305/1‑238 (MQ=255)                                                                                                                                                                                                           gTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTAt                                                    <  1:931127/211‑1 (MQ=255)                                                                                                                                                                                                             atacCGTTGCCGCGCCGCTGCGGGGGGGTGG‑GGCGACGT‑T‑GGTTTACGAAGGATGGCCGCCCGGGCGGGCCGGCGGCGGGTGGGGGACAATCGCCGAGAAATACCGGGTTAGCCGGATGTTCCCGGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTCCCGCTGAAATTCGCAAACCCGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGAc                       <  2:45609/238‑1 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    <  2:882167/172‑1 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    >  1:882167/1‑172 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAAGATTGCCGACCTGGCCGGACCGCGGCGCGTGGTGGGCAATCGTCGGGAAAAATCAGGTTAGCCGGATGTTTTCAGCGCCGGCCGCCATTCCCGTGCTGGAAAAATTCCCCACCGCTGAAATTTGCCAACACCATCTCTCGTCCCTGGAA‑GTGCCCTTTATGGTTGGAGAACCGCGGGAGGAGGCGACCGCCCGt                >  2:596298/1‑239 (MQ=255)                                                                                                                                                                                                                    ttACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGCCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGt                <  1:1016906/238‑1 (MQ=255)                                                                                                                                                                                                                      tcggggcgcTGGTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAGGTATGCCCGGCCTGGCCGGACGGCGGCGTGTGGGGGACAATCGGCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg             <  2:901869/233‑1 (MQ=255)                                                                                                                                                                                                                      aCGCGCCGCTGTGGGGGGGGAG‑GGGGACTT‑T‑CGTTTACGAAGGGTTGCCGACCGGGGCGGCCGGGGGCGTGGGGGGGACAATCGCCGAAAAAGAGCAGGTGAGCGGGATGTCCTCAGCGCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg             <  2:425855/239‑1 (MQ=255)                                                                                                                                                                                                                       cgcgCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGGCGGGAAATATCCGGTTAGCCGGCTGTTCTCAGCGCCCACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAACCACCATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCCGGGGAACCCCCGGGCGGGCCGGCcccccggttgg            >  2:34401/1‑232 (MQ=255)                                                                                                                                                                                                                             gcggctgGCGGGGAG‑GGCGACTA‑G‑CGTTTACGAAGGATTGCCGACCGGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCa      <  1:676130/236‑1 (MQ=255)                                                                                                                                                                                                                                  tgGCGGGGAG‑GGGGACTA‑C‑CGTTTACGAAGGTTTGCCGACTGGGCCGGATTGGGGCGGGTGGTGGACAACCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAc  <  2:942857/238‑1 (MQ=255)                                                                                                                                                                                                                                  tgGCGGGGAG‑GGCGACTA‑T‑CTTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAc  <  1:1116787/238‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAT‑GGCGACTA‑T‑CGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAC  >  NC_000913/353373‑353835

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