Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	353,610	C→G	13.0%	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	13.0%	145.4 / 3.6	66	A302G (GCG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base G (7/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAC  >  NC_000913/353373‑353835                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    gggTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGc                                                                                                                                                                                                                                     <  2:418146/238‑1 (MQ=255)     tccGGGGGCGTGGCTGGAATCCAACGAAACCGCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCGGCTGGCGGGGATGGCGACt                                                                                                                                                                                                                                 <  2:1230757/237‑1 (MQ=255)          tggcgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGt                                                                                                                                                                                                                           <  2:591367/238‑1 (MQ=255)           ggcgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGCGGCCC‑TCGttt                                                                                                                                                                                                                         >  2:1171926/1‑239 (MQ=255)             cgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGAATA‑TCGGTTAc                                                                                                                                                                                                                       >  1:1217163/1‑239 (MQ=255)             cgtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCCA‑TCGGTTAc                                                                                                                                                                                                                       >  1:127673/1‑239 (MQ=255)                          cAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGa                                                                                                                                                                                                          <  1:28679/239‑1 (MQ=255)                                aaCCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                     <  1:1035599/238‑1 (MQ=255)                                      ctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGGCCATTCGTTTGC‑AAGGGTGTCCGACCTGGCCgggcg                                                                                                                                                                                              >  2:268894/1‑236 (MQ=255)                                           ttCTCTCCACTCCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                         <  2:1030778/239‑1 (MQ=255)                                             ctctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACCAAGGATTGCCGACCTGGCCGGACTGCGGCgt                                                                                                                                                                                       >  1:250877/1‑239 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTA‑TCGTTTAAGCAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      >  1:1053242/1‑239 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTA‑TCGTTTTCGAAGGATTGCCGACCTGGCCGGACTgcggcggg                                                                                                                                                                                      >  1:727930/1‑237 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATTGCGACTA‑TCGTTTTCGGAGGATTGCCGACCTGGCCGGGCTgcggcggg                                                                                                                                                                                      >  1:525130/1‑237 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      <  1:1170925/239‑1 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                      <  1:263084/239‑1 (MQ=255)                                                tACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGAc                                                                                                                                                                                                                                  <  2:43658/194‑1 (MQ=255)                                                tACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGAc                                                                                                                                                                                                                                  >  1:43658/1‑194 (MQ=255)                                                      tCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGGCTA‑TCGTTTCCGAAGGGTTGGCGGCCTGGCCGGACGTCGGCGTGTGGGGGAc                                                                                                                                                                              >  1:1126449/1‑239 (MQ=255)                                                       ccGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTA‑TCGTTTACGAAGGGTTTCCCACCCGGCCCGACTGGGGCGTGTGGTGGGCa                                                                                                                                                                             >  1:311622/1‑239 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCTA‑TCGTTTACGCAGGATTGCCGACCCGGCCGGGCTGCGGCGTGTGGTGGTCaa                                                                                                                                                                            >  1:975555/1‑239 (MQ=255)                                                        cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGGTTACGAAGGATTGCCGACCTGGGCGGACTGCGGCGTGTGGTGGGCaa                                                                                                                                                                            >  1:942808/1‑239 (MQ=255)                                                                ccGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCga                                                                                                                                                                     <  2:489527/238‑1 (MQ=255)                                                                 cGGCAAACCTAAAGGCGGGCAGCGGGACGTCCGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCgaga                                                                                                                                                                   <  2:311622/239‑1 (MQ=255)                                                                         cTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAAtcgt                                                                                                                                                                        <  1:211812/226‑1 (MQ=255)                                                                          tAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAgg                                                                                                                                                          <  2:457148/239‑1 (MQ=255)                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGGCTA‑TCGTTTACGAAGGGTTGCCGACCTGGCCGGGCTGCGGCGTGTTGGGGGCAATCGTCGAGAAATATCCAGTTAGCTgg                                                                                                                                                  >  1:160641/1‑239 (MQ=255)                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGGCAATCGTCGAGAAATATCAGGTTAGcccg                                                                                                                                                  >  1:101088/1‑237 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGCGGCTA‑ACGGTTTCGAAGGATTGCCGACCCGGCCGGGGTTCGGCGGGGGGGGGGCAATCCGCGGGAAATAGCAGGttgaccgcg                                                                                                                                                 >  2:829248/1‑232 (MQ=255)                                                                                         tGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTA‑TCGTTTACGAAGGATTGCCGGCCTGGCCGGGCTGCGGGGGGGGGGGGGCAATTGTCGGGAAATAACAGGTTAGCCGGGTGGtct                                                                                                                                           >  1:45609/1‑239 (MQ=255)                                                                                                               gcgcTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGGCCGCCATTcgctg                                                                                                                     >  1:1040967/1‑237 (MQ=255)                                                                                                               gcgcTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACCA‑TCGTTTACGCAGGGTTGCCGCCCTGGCCCGACTGCGGGGTGTGGGGGAGAATCGTGGGGAAATATCTGGTTAGGCGGATGTTCTCAGCGCCGACCGcccattgcgt                                                                                                                     >  2:1142402/1‑231 (MQ=255)                                                                                                                         ccTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGATA‑TCGTTTACGGAGGGTTGCCCGCCTGGCCGGACTGCGGCGGGGGGTGGGCAATCTGCGGGAAATATCAGGGTTGCCCGGTGGTCTCATCGCCAGCCGCCCCTCGCGGGCTGagaaaa                                                                                                           >  2:983912/1‑239 (MQ=255)                                                                                                                               tGGACACCATTTTGGGCGGCAAAGCGGGCAGCGTGTTCTTTGGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                      <  1:599680/238‑1 (MQ=255)                                                                                                                                              gcggcAAGGCGGGCGGCGTGTTCTTTTGCGCATCGGATAGCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                       <  2:727930/238‑1 (MQ=255)                                                                                                                                               cggcAAAGCGGGCGGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGc                                                                                      <  1:2449/238‑1 (MQ=255)                                                                                                                                                ggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTTCGGCGTGTGGGGGACAATCGTCGAGAAATATCAGGGTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGCAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    >  1:901869/1‑239 (MQ=255)                                                                                                                                                ggaaaaGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCa                                                                                     <  1:268894/235‑1 (MQ=255)                                                                                                                                                            aGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTA‑TCGTTTACGAAGGATTTCCGAACTGGGCGGGCTGCGGCGCGTGGGGGGCAATCGTCGGGAAATATCAGGTTAGCCGGGTGTTTTCTGCGCCGGCCGCCATTCGCGTGCTGAAAAAAATCCCTACCGCTGTAATTCGCAAACACGAtctct                                                                         >  2:448096/1‑238 (MQ=255)                                                                                                                                                               gtgtTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                        <  1:993544/236‑1 (MQ=255)                                                                                                                                                                 gtTCTTTTGCGCATCGGATGTCGGCGGGGTGGTGGTGCATTCGTATACCGTTTCCGCGCCGCTGCTGGCGGAGAGGGCGACTA‑TCGTTTACGAAGGATTGCCGACCGGGCCGGACGGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGc                                                                   <  2:1042108/239‑1 (MQ=255)                                                                                                                                                                      tttGCGCATCGGTTATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                              <  2:1245928/239‑1 (MQ=255)                                                                                                                                                                         gcgcATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTg                                                          <  1:1142402/239‑1 (MQ=255)                                                                                                                                                                                  atagCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATct                                                  <  1:1171926/238‑1 (MQ=255)                                                                                                                                                                                  atagCGGCTGGGTGGGGGGGCATTGGTATATCGTTTACGCGCCGCGGCGGGCGGGGAGGGCGACTA‑TCGTTTACGAAGGATGGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTg                                                          <  1:866357/230‑1 (MQ=255)                                                                                                                                                                                    atCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATctg                                                 >  1:942857/1‑237 (MQ=255)                                                                                                                                                                                          tGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCTGCGGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGCCCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGgag                                          <  1:580851/238‑1 (MQ=255)                                                                                                                                                                                           gggTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGGGGCGTGTGGGGGACAAACGTCGAGAAATATTAGGGTAACCGGGTGGTCTCAGCGCCGACCGCCATTCGCGTGCCGAAAAAATTTCCCAACCCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGGTGC‑CTATCTGGATGGGGaa                                        >  2:1248934/1‑239 (MQ=255)                                                                                                                                                                                            ggtggtGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGGGGGGTGGGCAAACGTCGAGAAATATCAGGGTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGGGGTGAAAAAATTCCCTACCCCTGAAATTTCCAAAAACGAACTCTCGTCGCTGGGA‑GGGGTCCAActggcctgtttta                                       >  2:866357/1‑231 (MQ=255)                                                                                                                                                                                            ggtggtGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGGCCA‑TCGTTTACGAAGGATTGCCGACCCGGCCGGACTGCGGCGTGGGGGGGGCAATCGTCGCGGAATATCAGGGTTGCCGGGTGGTCTCAGCGCCGACCCCCCTTTCCGTGCCGGAAAAATTTCCTACCGCTGAAATTCGCAAACACGGTCTCTCCGCGCTGGGC‑GTGCTCTCTCTGGGTGGACaa                                        >  2:813514/1‑238 (MQ=255)                                                                                                                                                                                                 tggGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGGCCATCCGCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCCGAAATTCGCAAACACGATCTCTCGCCCCGGGGA‑GTGCTCTATCTGGCTGGGGAACCGccg                                  >  2:19622/1‑237 (MQ=255)                                                                                                                                                                                                   gggCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTgg                                 >  1:11297/1‑238 (MQ=255)                                                                                                                                                                                                    ggCATTCGTATATCGTTTACGCGCCGCTGTGGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGa                                <  1:754165/238‑1 (MQ=255)                                                                                                                                                                                                      cATTCGTTTATCGTTTACGCGCCGCTGCTGGGGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACg                              <  1:312168/238‑1 (MQ=255)                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGTGGGCAATCGTCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCGCCGACCGCCCTTTGCGTGCTGAAAAAATTCCCCACCGCTGAAATTCGCCAACACGAACTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAg                            >  2:877305/1‑238 (MQ=255)                                                                                                                                                                                                           gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTAt                                                    <  1:931127/211‑1 (MQ=255)                                                                                                                                                                                                             atacCGTTGCCGCGCCGCTGCGGGGGGGTGGGGCGACGT‑TGGTTTACGAAGGATGGCCGCCCGGGCGGGCCGGCGGCGGGTGGGGGACAATCGCCGAGAAATACCGGGTTAGCCGGATGTTCCCGGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTCCCGCTGAAATTCGCAAACCCGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGAc                       <  2:45609/238‑1 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    >  1:882167/1‑172 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCaa                                                                                    <  2:882167/172‑1 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAAGATTGCCGACCTGGCCGGACCGCGGCGCGTGGTGGGCAATCGTCGGGAAAAATCAGGTTAGCCGGATGTTTTCAGCGCCGGCCGCCATTCCCGTGCTGGAAAAATTCCCCACCGCTGAAATTTGCCAACACCATCTCTCGTCCCTGGAA‑GTGCCCTTTATGGTTGGAGAACCGCGGGAGGAGGCGACCGCCCGt                >  2:596298/1‑239 (MQ=255)                                                                                                                                                                                                                    ttACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGCCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGt                <  1:1016906/238‑1 (MQ=255)                                                                                                                                                                                                                      tcggggcgcTGGTGGCGGGGATGGCGACTA‑TCGTTTACGAGGTATGCCCGGCCTGGCCGGACGGCGGCGTGTGGGGGACAATCGGCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg             <  2:901869/233‑1 (MQ=255)                                                                                                                                                                                                                      aCGCGCCGCTGTGGGGGGGGAGGGGGACTT‑TCGTTTACGAAGGGTTGCCGACCGGGGCGGCCGGGGGCGTGGGGGGGACAATCGCCGAAAAAGAGCAGGTGAGCGGGATGTCCTCAGCGCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg             <  2:425855/239‑1 (MQ=255)                                                                                                                                                                                                                       cgcgCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGGCGGGAAATATCCGGTTAGCCGGCTGTTCTCAGCGCCCACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAACCACCATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCCGGGGAACCCCCGGGCGGGCCGGCcccccggttgg            >  2:34401/1‑232 (MQ=255)                                                                                                                                                                                                                             gcggctgGCGGGGAGGGCGACTA‑GCGTTTACGAAGGATTGCCGACCGGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCa      <  1:676130/236‑1 (MQ=255)                                                                                                                                                                                                                                  tgGCGGGGAGGGGGACTA‑CCGTTTACGAAGGTTTGCCGACTGGGCCGGATTGGGGCGGGTGGTGGACAACCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAc  <  2:942857/238‑1 (MQ=255)                                                                                                                                                                                                                                  tgGCGGGGAGGGCGACTA‑TCTTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑GTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATAc  <  1:1116787/238‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTA‑TCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAA‑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 < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑