Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,570,469	A→C	7.9%	intergenic (+35/‑153)	pcaF‑I → / → pcaT	beta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,570,469	0	A	C	7.9%	192.3 / 7.0	63	intergenic (+35/‑153)	pcaF‑I/pcaT	beta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (16/41);  new base C (3/2);  total (19/44)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

ACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCC  >  NC_002947/1570320‑1570618                                                                                                                                                      |                                                                                                                                                       aCGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCACCa                                                                                                                                                        >  1:1585944/1‑150 (MQ=255) aCGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTAGCGGGGAAGCCCCCCCCa                                                                                                                                                        >  1:1738719/1‑150 (MQ=255)    ccTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACa                                                                                                                                                     <  2:4762817/150‑1 (MQ=255)     cTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAg                                                                                                                                                    <  2:482239/150‑1 (MQ=255)     cTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAg                                                                                                                                                    <  2:482177/150‑1 (MQ=255)     cTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAg                                                                                                                                                    <  1:2187226/150‑1 (MQ=255)      tGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAgg                                                                                                                                                   <  1:1732907/150‑1 (MQ=255)      tGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCAGGTGTGTGGGTGTCGGCCAAGGTCTGGCGGTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAgg                                                                                                                                                   >  2:3704474/1‑150 (MQ=255)       ggTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAggt                                                                                                                                                  <  1:238073/150‑2 (MQ=255)       ggTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAggt                                                                                                                                                  <  1:238151/150‑2 (MQ=255)           ctgactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATcc                                                                                                                                             <  1:4150695/150‑1 (MQ=255)           ctgactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATcc                                                                                                                                             <  2:3285572/150‑1 (MQ=255)             gactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCAc                                                                                                                                           <  2:2657819/150‑1 (MQ=255)                tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCa                                                                                                                                        >  1:5010891/1‑150 (MQ=255)                tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCa                                                                                                                                        <  2:218499/150‑1 (MQ=255)                   gTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATc                                                                                                                                     >  2:1677458/1‑150 (MQ=255)                       cACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCt                                                                                                                                 <  2:3159846/150‑1 (MQ=255)                           aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTgg                                                                                                                             <  2:572234/150‑1 (MQ=255)                                 agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCgg                                                                                                                       >  2:3719426/1‑150 (MQ=255)                                 agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCgg                                                                                                                       <  1:1081734/150‑1 (MQ=255)                                 agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCgg                                                                                                                       >  2:3719434/1‑150 (MQ=255)                                 agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCgg                                                                                                                       <  2:219912/150‑1 (MQ=255)                                  gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCggg                                                                                                                      <  1:683166/150‑1 (MQ=255)                                  gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCggg                                                                                                                      <  2:379705/150‑1 (MQ=255)                                  gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCGGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCggg                                                                                                                      <  1:683196/150‑1 (MQ=255)                                   aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                     <  2:1339842/150‑1 (MQ=255)                                   aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                     <  2:2417799/150‑1 (MQ=255)                                          gCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCACGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCcg                                                                                                              <  2:4408187/150‑1 (MQ=255)                                          gCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGACAAGTTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCTCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCcg                                                                                                              <  2:4408445/150‑1 (MQ=255)                                                gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                                        <  1:1677458/150‑1 (MQ=255)                                                gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCACCAACAGC‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                                        <  2:2468021/150‑1 (MQ=255)                                                  aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAgg                                                                                                      <  2:5084452/150‑1 (MQ=255)                                                         tGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGATCCCCCGCCAACAGG‑GGTCCACGCCCACACCGGGGGGGGCGGGCGGCCCCCCGCAGCGGccgcgcc                                                                                               >  1:3121550/1‑145 (MQ=255)                                                             gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACa                                                                                                                                                     >  2:4624346/1‑93 (MQ=255)                                                             gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACa                                                                                                                                                     <  1:4624346/93‑1 (MQ=255)                                                               cctggggtgtgGGTGTCGGCAATGGTGGGGGGTGGGCCATCGCGGGGGTTTGACCTGTCCCGGCCTCTTCGCGGGGAAGCCCGCGCCAACTGG‑GATCCACGAAATCACCGGTGGGAGCGGGCTCCCCCGCGAAGAGGCCTCTCCTGTgag                                                                                         <  1:4157683/144‑1 (MQ=255)                                                               ccATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTgag                                                                                         <  2:5012552/150‑1 (MQ=255)                                                                      tgtgGGTGTCGGCCACGGTGTGGGGTTGGCCATCGCGCGGGTTTGACCTGTACCGGACTCTTCGGGGGGACGCCCGCTCCACCAGG‑GATCCACGCAATCACCTGGGGGCGGGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAg                                                                                  <  2:2293711/150‑1 (MQ=255)                                                                             gtCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCAtt                                                                           <  2:1984153/150‑1 (MQ=255)                                                                                ggCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGcc                                                                        >  2:3816128/1‑150 (MQ=255)                                                                                   cAATGTATGGCGTTGGCCATCGAGCTGGTTAGACCAGTAACGGCCTCTTCGCGGGGAAGCCCTCACCAAAAGG‑GATACACGCAATAACCGGTGGGAGCGGGCTTCCCCGCGAAGAGGCCGCTCCTGGGAGACAAAAGGCCCATTGCCAGc                                                                     >  2:1181951/1‑150 (MQ=255)                                                                                   cAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGc                                                                     >  2:1181953/1‑150 (MQ=255)                                                                                   cAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCACCTGTGAGACATAAGGCCCATTGCCAGc                                                                     >  2:1181892/1‑150 (MQ=255)                                                                                         cgggggTTGGCCTTGGCGGGGGTTTGACCTGTACCGGCCTCGTCGCGGGGAAGCCCGCTCCAACCGG‑GATCCCCGAAATCACGGGTGGGCGCGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATa                                                               <  1:3742646/145‑1 (MQ=255)                                                                                               ttGGCCAGCGCGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑TATCCACGCAATCCCCTGTGGGCGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGaa                                                         <  2:3517709/150‑1 (MQ=39)                                                                                                              gTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAg                                          <  2:2404338/150‑1 (MQ=255)                                                                                                                ttGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCCCCTGTGGGCGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGt                                        <  2:1738719/150‑1 (MQ=255)                                                                                                                    ccTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCTACAGG‑GTTCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGTCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATATCTGAACGAGTGTGTTACTAGGTATGt                                    <  2:4528012/150‑1 (MQ=255)                                                                                                                                tctTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATC‑CCTGTTGGAGCGGGCTTCCCCACGAAGAGGCCTCTCCTGTGAGACATAATGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                        >  2:485337/1‑150 (MQ=18)                                                                                                                                tctTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATC‑CCTGTTGGAGCGGGCTTCCCCACGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                        >  2:485466/1‑150 (MQ=18)                                                                                                                                tctTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATC‑CCTGTTGGAGCGGGCTTCCCCACGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                        >  2:485582/1‑150 (MQ=18)                                                                                                                                tctTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                        >  2:4070817/1‑150 (MQ=255)                                                                                                                                 ttgggcgGGGACGCCCGCACCAACCGG‑GACCCCCGCAACCCCGGGTGGGCGCGGGCGGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa                       <  1:3549236/146‑1 (MQ=255)                                                                                                                                 ctTCGCGGGGAAGCCCGCTCCAACCGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAGTTCTAGACATACCTa                       <  2:3219532/150‑1 (MQ=255)                                                                                                                                 ctTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa                       <  2:2447336/150‑1 (MQ=255)                                                                                                                                           aaGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGTGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCa             <  1:799517/150‑1 (MQ=255)                                                                                                                                            aGCCCGCCCCAACAGG‑GATCCACGCACTCACCTGTGGGTGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg            <  2:3121550/150‑1 (MQ=255)                                                                                                                                             gCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAgg           >  1:2921389/1‑150 (MQ=255)                                                                                                                                              cccGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa          <  2:389887/150‑1 (MQ=255)                                                                                                                                              cccGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa          <  1:2139901/150‑1 (MQ=255)                                                                                                                                              cccGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa          >  1:1815153/1‑150 (MQ=255)                                                                                                                                              cccGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa          <  2:3268721/150‑1 (MQ=255)                                                                                                                                              cccGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa          >  2:2139901/1‑150 (MQ=255)                                                                                                                                                  cTCCAACCGG‑GAGCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTaa      <  1:4564843/150‑1 (MQ=255)                                                                                                                                                   tCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCAAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTaaa     >  1:4771099/1‑150 (MQ=255)                                                                                                                                                      aaCAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAccc  >  2:603458/1‑150 (MQ=255)                                                                                                                                                      aaCAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAcac  >  1:4294785/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGG‑GATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCC  >  NC_002947/1570320‑1570618

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