Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,935,858	C→G	8.6%	intergenic (+153/‑109)	PP_2568 → / → PP_2569	ring‑cleaving dioxygenase/MHS family metabolite MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,935,858	0	C	G	8.6%	259.0 / 11.9	81	intergenic (+153/‑109)	PP_2568/PP_2569	ring‑cleaving dioxygenase/MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (3/4);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCC  >  NC_002947/2935712‑2936005                                                                                                                                                     |                                                                                                                                                    ccccGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGccc                                                                                                                                                   >  6:363170/1‑149 (MQ=255)      ccAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtc                                                                                                                                              <  7:247548/149‑1 (MQ=255)         aaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGTGGCCCGCCtcgtc                                                                                                                                           <  4:189310/149‑1 (MQ=255)         aaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                           <  8:134758/149‑1 (MQ=255)           aTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcggatc                                                                                                                                         <  5:2658/149‑5 (MQ=255)               tCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCGCAGGGCCCGGCGGGTGGAAATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGc                                                                                                                                       >  1:291652/1‑149 (MQ=255)                  gCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTcgc                                                                                                                                  <  5:72126/149‑1 (MQ=255)                  gCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTcgc                                                                                                                                  <  6:332566/149‑1 (MQ=255)                   ccAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTcgcg                                                                                                                                 <  3:352740/149‑1 (MQ=255)                     aaCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCgcgc                                                                                                                                >  3:122198/1‑148 (MQ=255)                     aaCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCgcgc                                                                                                                                >  4:330038/1‑148 (MQ=255)                      aCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCaa                                                                                                                              <  4:202095/149‑1 (MQ=255)                       cccGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAAc                                                                                                                             >  4:224547/1‑149 (MQ=255)                          gCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCa                                                                                                                           >  5:282646/1‑148 (MQ=255)                             aTCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaag                                                                                                                        >  8:194001/1‑148 (MQ=255)                                gCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagcaa                                                                                                                     <  7:194001/148‑1 (MQ=255)                                      cACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAg                                                                                                              <  8:114526/149‑1 (MQ=255)                                      cACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAg                                                                                                              <  7:233073/149‑1 (MQ=255)                                         cccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCCGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGc                                                                                                           >  1:186835/1‑149 (MQ=255)                                         cccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGc                                                                                                           >  5:42209/1‑149 (MQ=255)                                          ccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGcc                                                                                                          >  1:327831/1‑149 (MQ=255)                                                  gcGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCg                                                                                                  >  3:37462/1‑149 (MQ=255)                                                   cGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgt                                                                                                                                            >  4:223876/1‑106 (MQ=255)                                                   cGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgt                                                                                                                                            <  3:223876/106‑1 (MQ=255)                                                   cGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgg                                                                                                       <  8:45554/143‑1 (MQ=255)                                                   cGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgg                                                                                                       >  7:45554/1‑143 (MQ=255)                                                     tgtCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACg                                                                                               >  6:141674/1‑149 (MQ=255)                                                      gtcgtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACgg                                                                                              >  8:217993/1‑149 (MQ=255)                                                       tcgtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGt                                                                                             <  4:37462/149‑1 (MQ=255)                                                       tcgtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGt                                                                                             <  5:363170/149‑1 (MQ=255)                                                         gtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagcaagc                                                                                                                   >  6:211906/1‑125 (MQ=255)                                                         gtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagcaagc                                                                                                                   <  5:211906/125‑1 (MQ=255)                                                         gtcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTg                                                                                            <  4:267271/148‑1 (MQ=255)                                                          tcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGaa                                                                                          <  4:162217/149‑1 (MQ=255)                                                          tcCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGaa                                                                                          <  3:224547/149‑1 (MQ=255)                                                           ccAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAg                                                                                         >  3:79775/1‑149 (MQ=255)                                                               gccgccCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGt                                                                                     >  4:103862/1‑149 (MQ=255)                                                               gccgccCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGt                                                                                     <  6:42209/149‑1 (MQ=255)                                                                     ctctCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAg                                                                               <  7:217993/149‑1 (MQ=255)                                                                      tctcGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAgg                                                                              >  5:96113/1‑149 (MQ=255)                                                                       ctcGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAg                                                                               <  3:304597/147‑1 (MQ=255)                                                                       ctcGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAg                                                                               >  4:304597/1‑147 (MQ=255)                                                                          gT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTg                                                                            <  5:141674/147‑1 (MQ=255)                                                                          gT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCt                                                                          >  3:279167/1‑149 (MQ=255)                                                                                     ccGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTg                                                                  >  3:294799/1‑148 (MQ=255)                                                                                      cggcggGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTa                                                                <  2:222730/149‑1 (MQ=255)                                                                                            gTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTAGTGAGGTACTACGGCCTGTATGacaa                                                          >  4:65851/1‑149 (MQ=255)                                                                                             tGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaac                                                         <  4:79775/149‑1 (MQ=255)                                                                                              ggATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaaca                                                        >  8:355/1‑149 (MQ=255)                                                                                                 tatCCAGCCAATCCCGATTATGCTTCATGAATCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGa                                                     >  6:201872/1‑149 (MQ=255)                                                                                                   tCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACg                                                   >  2:48991/1‑149 (MQ=255)                                                                                                     cAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGAGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGATCCGGCCTGTATGACAAAATGACg                                                   >  2:371240/1‑147 (MQ=255)                                                                                                      aGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAg                                                >  8:148340/1‑149 (MQ=255)                                                                                                      aGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAg                                                <  8:316708/149‑1 (MQ=255)                                                                                                       gCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAgg                                               <  2:327831/149‑1 (MQ=255)                                                                                                       gCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAgg                                               <  7:148340/149‑1 (MQ=255)                                                                                                        ccAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAgg                                               >  3:227112/1‑148 (MQ=255)                                                                                                            tCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAcc                                          <  4:279167/149‑1 (MQ=255)                                                                                                            tCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAcc                                          <  4:227112/149‑1 (MQ=255)                                                                                                              ccGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCg                                                                       >  1:262826/1‑118 (MQ=255)                                                                                                              ccGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCg                                                                       <  2:262826/118‑1 (MQ=255)                                                                                                               cGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCAt                                       <  2:186835/149‑1 (MQ=255)                                                                                                                 aTTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATgg                                     <  3:103862/149‑1 (MQ=255)                                                                                                                 aTTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATgg                                     <  3:330038/149‑1 (MQ=255)                                                                                                                     tGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCgg                                  <  2:291652/148‑1 (MQ=255)                                                                                                                           aTGAACCTTGCTGCGATCCGAGGCGGCCCACCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGa                           >  1:27519/1‑149 (MQ=255)                                                                                                                             gAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGa                           >  3:174173/1‑147 (MQ=255)                                                                                                                               aCCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGc                       <  6:282646/149‑1 (MQ=255)                                                                                                                                ccTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGca                      >  8:35551/1‑149 (MQ=255)                                                                                                                                  ttGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGa                                                 <  2:355400/120‑1 (MQ=255)                                                                                                                                  ttGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGa                                                 >  1:355400/1‑120 (MQ=255)                                                                                                                                  ttGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGcaca                    >  3:241257/1‑149 (MQ=255)                                                                                                                                      tgcGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCa                >  6:246405/1‑149 (MQ=255)                                                                                                                                        cGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCAc               <  5:201872/148‑1 (MQ=255)                                                                                                                                          cgacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGa                                                                                <  4:21855/78‑1 (MQ=255)                                                                                                                                          cgacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGa                                                                                >  3:21855/4‑81 (MQ=255)                                                                                                                                           gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaaca                                                        <  2:89392/102‑1 (MQ=255)                                                                                                                                           gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaaca                                                        >  1:89392/3‑104 (MQ=255)                                                                                                                                           gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACa                         >  3:26828/3‑135 (MQ=255)                                                                                                                                           gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACa                         <  4:26828/133‑1 (MQ=255)                                                                                                                                              gAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCg                            <  5:295290/129‑1 (MQ=255)                                                                                                                                              gAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCg                            >  6:295290/1‑129 (MQ=255)                                                                                                                                                    gcCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGccc  <  7:355/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGT‑‑GCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCC  >  NC_002947/2935712‑2936005

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