Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,037,773	T→G	5.7%	intergenic (+15/+101)	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,773	0	T	G	5.7%	194.5 / 4.6	70	intergenic (+15/+101)	PP_0896/PP_0897	carbon‑nitrogen hydrolase family protein/class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

AGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGG  >  NC_002947/1037631‑1037922                                                                                                                                               |                                                                                                                                                       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)                                                           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                                                                                   >  1:383973/1‑110 (MQ=255)                                                                                                      ccGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCtt                                                                                   <  2:383973/110‑1 (MQ=255)                                                                                                             aGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGAc                                    <  2:748037/150‑1 (MQ=255)                                                                                                               ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAgg                                                   >  2:113458/1‑133 (MQ=255)                                                                                                               ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAgg                                                   <  1:113458/133‑1 (MQ=255)                                                                                                               ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                                  <  2:663413/150‑1 (MQ=255)                                                                                                               ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                                  >  1:1005472/1‑150 (MQ=255)                                                                                                               ttCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACgg                                  >  1:113870/1‑150 (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                                                                                              <  1:13763‑M1/71‑1 (MQ=255)                                                                                                                                 taACTGCGGCGGCTGCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTgg                                                                                              >  2:13763‑M1/3‑73 (MQ=255)                                                                                                                                   ACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTg              >  1:141015‑M1/1‑150 (MQ=255)                                                                                                                                   ACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTg              <  1:366479‑M1/150‑1 (MQ=255)                                                                                                                                     TGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGct            <  2:706593‑M1/150‑1 (MQ=255)                                                                                                                                      GTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGctc           >  1:1151324‑M1/1‑150 (MQ=255)                                                                                                                                         gcggcCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGc        <  1:363127/150‑1 (MQ=255)                                                                                                                                           ggcCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTac                                               <  1:293852/109‑1 (MQ=255)                                                                                                                                           ggcCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTac                                               >  2:293852/1‑109 (MQ=255)                                                                                                                                             ccTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCcaca    >  2:708081/1‑150 (MQ=255)                                                                                                                                             ccTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCcaca    >  2:865327/1‑150 (MQ=255)                                                                                                                                               tctTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCAcagg  <  1:255339/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       AGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCCAAATGCGTAGCCTTCTGGCTTGCGCTGTACGTGTGGGAGCGG‑GCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGG  >  NC_002947/1037631‑1037922

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