Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,037,765	T→C	8.3%	intergenic (+7/+109)	PP_0896 → / ← PP_0897	carbon‑nitrogen hydrolase family protein/class 1 fumarate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,037,765	0	T	C	8.3%	193.6 / 5.2	68	intergenic (+7/+109)	PP_0896/PP_0897	carbon‑nitrogen hydrolase family protein/class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (3/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTC  >  NC_002947/1037618‑1037913                                                                                                                                                    |                                                                                                                                                      gCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTgg                                                                                                                                                     >  1:205713/1‑150 (MQ=255)              agcTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTcgcg                                                                                                                                        <  2:192613/150‑1 (MQ=255)                     cccGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGc                                                                                                                                 >  1:980503/1‑150 (MQ=255)                     cccGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGCGGCGGCTGCTTCGCGGGCAAGc                                                                                                                                 >  1:231182/1‑150 (MQ=255)                     aaggacatgccggtgaccgtcgccgtggacagcaacggtgaatcggtacacatcactggccctgccctgtggcagagcaagattgccgaaagcctggccgtcgaagtgaagtaacctgatgtaACTGTGGCGGCCTCTTCGCGGGCAAGc                                                                                                                                 <  1:1069556‑M1/27‑1 (MQ=255)                           cTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGCGGCGGCTGCTTCGCGGGCAAGCCCGCTc                                                                                                                           <  2:231182/150‑1 (MQ=255)                              cAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCa                                                                                                                        >  1:663413/1‑150 (MQ=255)                               ggCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCaa                                                                                                                       >  1:994547/2‑150 (MQ=255)                                 ccAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAAt                                                                                                                     >  2:573531/1‑150 (MQ=255)                                 ccAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAAt                                                                                                                     >  2:1186730/1‑150 (MQ=255)                                   aGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGc                                                                                                                   >  2:363127/1‑150 (MQ=255)                                     tGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGt                                                                                                                 <  2:107790/150‑1 (MQ=255)                                      ggCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGAAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCTGGCAAGCCCGCTCCCAAATGCGTa                                                                                                                <  2:619008/150‑1 (MQ=255)                                       gCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGGGCCTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAg                                                                                                               <  1:979176/150‑1 (MQ=255)                                         tGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGcc                                                                                                             >  2:1156246/1‑150 (MQ=255)                                         tGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGCGCTGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGcc                                                                                                             >  2:141015/1‑150 (MQ=255)                                                tgAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGc                                                                                                      >  1:225062/1‑150 (MQ=255)                                                tgAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGc                                                                                                      >  1:985805/1‑150 (MQ=255)                                                 gAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGGGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCt                                                                                                     >  1:728884/1‑150 (MQ=255)                                                 gAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCt                                                                                                     >  1:741196/1‑150 (MQ=255)                                                    gTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTg                                                                                                        >  1:233538/1‑144 (MQ=255)                                                    gTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTg                                                                                                        <  2:233538/144‑1 (MQ=255)                                                       gATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCt                                                                                               <  2:985805/150‑1 (MQ=255)                                                           tGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTAc                                                                                           <  2:3318/150‑1 (MQ=255)                                                               ggaAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACgtgt                                                                                       <  2:1093021/150‑1 (MQ=255)                                                                 aaGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTgg                                                                                     <  1:573531/150‑1 (MQ=255)                                                                       ggcggGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGCCTGCCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCgg                                                                               <  1:1250285/150‑1 (MQ=255)                                                                        gcggGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTggcg                                                                                                                                                   >  1:529203/1‑81 (MQ=255)                                                                        gcggGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTggcg                                                                                                                                                   <  2:529203/81‑1 (MQ=255)                                                                        atcactggccctgccctgtggcagagcaagattgccgaaagcctggccgtcgaagtgaagtaacctgatgtaACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCgg‑g                                                                             <  2:194739‑M1/78‑1 (MQ=255)                                                                         tcactggccctgccctgtggcagagcaagattgccgaaagcctggccgtcgaagtgaagtaacctgatgtaACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑Gc                                                                            >  2:1069556‑M1/72‑150 (MQ=255)                                                                         cggGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑Gc                                                                            <  2:741196/150‑1 (MQ=255)                                                                           actggccctgccctgtggcagagcaagattgccgaaagcctggccgtcgaagtgaagtaacctgatgtaACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                          >  1:194739‑M1/70‑150 (MQ=255)                                                                             tGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGc                                                                        <  1:116539/150‑1 (MQ=255)                                                                                  cGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCAAGCCcgct                                                                   <  1:1156246/150‑2 (MQ=255)                                                                                    ggAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGaa                                                                 <  2:205713/150‑1 (MQ=255)                                                                                         ggCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGCagcagc                                                            >  2:840399/1‑150 (MQ=255)                                                                                         ggCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAagcagc                                                            <  2:735531/150‑1 (MQ=255)                                                                                          gCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAgcc                                                           >  1:1128688/1‑150 (MQ=255)                                                                                           cAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAgccg                                                          >  2:1233080/1‑150 (MQ=255)                                                                                           cAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCGTGCCCGCGAAGCAgccg                                                          >  1:65450/1‑150 (MQ=255)                                                                                                agcaagattgccgaaagcctggccgtcgaagtgaagtaacctgatgtaACTGCGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGCCCGATTGG‑TCTTGCCCGCGAAGCAGCCGCCGCa                                                     >  1:33220‑M1/49‑150 (MQ=255)                                                                                                     ccTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTAc                                                >  1:748037/1‑150 (MQ=255)                                                                                                         gCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCa                                            <  2:225062/150‑1 (MQ=255)                                                                                                            ggACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGt                                         <  2:344750/150‑1 (MQ=255)                                                                                                              aCTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTa                                       <  2:65450/150‑1 (MQ=255)                                                                                                               cTGGCCGGGGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTac                                      <  1:484489/150‑1 (MQ=255)                                                                                                               cTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACATG‑GGGCGTGGTGGGTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTac                                      <  1:1186730/150‑1 (MQ=255)                                                                                                                 ggCCGGAGCAGTTCGTGAGCGTCCTGTGGTCCCACCAGTCGACTGTCAGCGGGGATTCCCGTTCCCAGATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGACAGG‑GCTTGCCCGCGACGCAGCCGCCGCAGTTACATCAGGTTactt                                    <  2:920281/150‑1 (MQ=255)                                                                                                                  gCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                          <  2:627636/111‑1 (MQ=255)                                                                                                                  gCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                          >  1:627636/1‑111 (MQ=255)                                                                                                                   ccGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                          <  2:383973/110‑1 (MQ=255)                                                                                                                   ccGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                          >  1:383973/1‑110 (MQ=255)                                                                                                                          aGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGAc                           <  2:748037/150‑1 (MQ=255)                                                                                                                            ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAgg                                          <  1:113458/133‑1 (MQ=255)                                                                                                                            ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAgg                                          >  2:113458/1‑133 (MQ=255)                                                                                                                            ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                         >  1:113870/1‑150 (MQ=255)                                                                                                                            ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                         >  1:1005472/1‑150 (MQ=255)                                                                                                                            ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                         <  2:663413/150‑1 (MQ=255)                                                                                                                                 gagagTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAgg                    >  2:399339/1‑150 (MQ=255)                                                                                                                                  agagTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGc                   >  2:1151324/1‑150 (MQ=255)                                                                                                                                        cTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTTCTTCACTTCGACGGCCAGGCTTTCgg             >  2:551185/1‑150 (MQ=255)                                                                                                                                              taACTGCGGCGGCTGCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTgg                                                                                     >  2:13763‑M1/3‑73 (MQ=255)                                                                                                                                              taACTGCGGCGGCTGCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTgg                                                                                     <  1:13763‑M1/71‑1 (MQ=255)                                                                                                                                                ACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTg     <  1:366479‑M1/150‑1 (MQ=255)                                                                                                                                                ACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTg     >  1:141015‑M1/1‑150 (MQ=255)                                                                                                                                                  TGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGct   <  2:706593‑M1/150‑1 (MQ=255)                                                                                                                                                   GTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGctc  >  1:1151324‑M1/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTC  >  NC_002947/1037618‑1037913

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