Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,863,777	G→C	7.6%	P900P (CCG→CCC)	PP_3413 →	two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,863,777	0	G	C	7.6%	209.0 / 5.8	66	P900P (CCG→CCC)	PP_3413	two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base C (3/2);  total (33/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 = 1.00e+00

CGACGACTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGC  >  NC_002947/3863630‑3863921                                                                                                                                                    |                                                                                                                                                   cgacgaCTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCgg                                                                                                                                                   >  2:275047/1‑149 (MQ=255) cgacgaCTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCgg                                                                                                                                                   >  2:230006/1‑149 (MQ=255)          cccGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGAc                                                                                                                                           <  6:367080/148‑1 (MQ=255)           ccGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACgg                                                                                                                                         >  6:74129/1‑149 (MQ=255)               ccAACCTGCCTTTGCTGGAACGGCAACTGCAAACCCTGCGGCCCCATGTGACACTTGCGGTGAATGGCGAGATCGGCCTTGCACGCTGGCAGGAAGAGCGCTTCGACCTGGTGAGCACCGATTGCAGCATGCCCGTCCATGACCGGCAt                                                                                                                                     <  6:367588/149‑1 (MQ=255)                cAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAt                                                                                                                                     >  5:315464/1‑148 (MQ=255)                    cTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgctg                                                                                                                                 >  8:115071/1‑145 (MQ=255)                    cTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTGCGCCCTGGTGAGCACCGATTGCAGCATGCCGGTCATGGACGGGCATGATCt                                                                                                                                >  4:15138/1‑149 (MQ=255)                     tGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTc                                                                                                                               >  2:209142/1‑149 (MQ=255)                         tttGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCAccc                                                                                                                           >  6:202665/1‑149 (MQ=255)                            gCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                  <  2:216125/139‑3 (MQ=255)                            gCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                  >  1:216125/1‑137 (MQ=255)                             cTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  6:165924/1‑149 (MQ=255)                              tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  5:299623/1‑148 (MQ=255)                              tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  7:86986/1‑148 (MQ=255)                              tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  5:268386/1‑148 (MQ=255)                                     gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCg                                                                                                                 >  3:628/1‑147 (MQ=255)                                     gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCa                                                                                                               <  6:268386/149‑1 (MQ=255)                                      cAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAg                                                                                                              >  2:202994/1‑149 (MQ=255)                                       aaCTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGc                                                                                                             <  5:74129/149‑1 (MQ=255)                                            cAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑g                                                                                                         <  4:628/147‑1 (MQ=255)                                             aaaCCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAg                                                                                                      <  7:301005/149‑1 (MQ=255)                                                 ccTGGGGCACCATGAGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCg                                                                                                  <  7:135553/149‑1 (MQ=255)                                                       gCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGcgcg                                                                                            >  2:96379/1‑149 (MQ=255)                                                       gCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGcgcg                                                                                            <  3:174860/149‑1 (MQ=255)                                                               tgACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGc                                                                                      >  8:288503/1‑147 (MQ=255)                                                                   acTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCcagcag                                                                                >  8:191835/1‑149 (MQ=255)                                                                     tGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGaa                                                                                                                                  >  2:341941/1‑98 (MQ=255)                                                                       gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGc                                                                              <  2:283875/147‑1 (MQ=255)                                                                       gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCg                                                                            >  7:21518/1‑149 (MQ=255)                                                                       gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGGCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTCGGAGGGCGAGCGCGGCCTGCCAGCAGGcc                                                                             <  1:209142/149‑1 (MQ=255)                                                                        cGGAGAATGGCGAGATCGCCCTTGCATGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGTCATGAACTCACCCGGCGCGTTCGCAGCCG‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCg                                                                            >  4:235052/1‑148 (MQ=255)                                                                          gagaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAt                                                                         <  7:191835/149‑1 (MQ=255)                                                                            gaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTc                                                                       <  3:14743/149‑1 (MQ=255)                                                                               tGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTg                                                                    <  6:299623/149‑1 (MQ=255)                                                                               tGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTg                                                                    >  3:108144/1‑149 (MQ=255)                                                                                 gCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑c                                                                 <  8:86986/149‑1 (MQ=255)                                                                                         gCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGc                                                         >  6:69536/1‑149 (MQ=255)                                                                                          cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGc                                                         >  6:153765/1‑148 (MQ=255)                                                                                               gCACGCTGTCAGGAAGCGCGCTGCGACCTGGTGATCACCGATTGCAGCATGCCTGTCCTGGCCGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATg                                                    <  8:330556/148‑1 (MQ=255)                                                                                                  cGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCa                                                <  5:202665/149‑1 (MQ=255)                                                                                                   gCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTc                                                                       <  2:309145/126‑1 (MQ=255)                                                                                                        cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAt                                                                         >  1:74282/1‑119 (MQ=255)                                                                                                        cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTAt                                                                         <  2:74282/119‑1 (MQ=255)                                                                                                        cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCga                                          <  6:315464/149‑1 (MQ=255)                                                                                                        cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCga                                          <  5:69536/149‑1 (MQ=255)                                                                                                         aGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCgag                                         <  7:288503/149‑1 (MQ=255)                                                                                                          ggAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCtt                                                                     <  4:42559/121‑1 (MQ=255)                                                                                                          ggAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCtt                                                                     >  3:42559/1‑121 (MQ=255)                                                                                                             aGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGc                                     <  5:153765/149‑1 (MQ=255)                                                                                                                  gcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACt                                <  1:275047/149‑1 (MQ=255)                                                                                                                  gcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACt                                >  1:246784/1‑149 (MQ=255)                                                                                                                   cTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACt                                <  1:230006/149‑1 (MQ=255)                                                                                                                   cTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGt                                                               >  7:332033/1‑117 (MQ=255)                                                                                                                   cTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGt                                                               <  8:332033/117‑1 (MQ=255)                                                                                                                    ttCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCCGGCCGAGGAGCGTACTCg                              >  2:116438/1‑149 (MQ=255)                                                                                                                       gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCt                            <  2:246784/148‑1 (MQ=255)                                                                                                                        aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                          <  4:108144/149‑1 (MQ=255)                                                                                                                                  cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTc                                                                       >  4:48430/1‑95 (MQ=255)                                                                                                                                  cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTc                                                                       <  3:48430/95‑1 (MQ=255)                                                                                                                                  cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCg                >  3:79513/1‑149 (MQ=255)                                                                                                                                   aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCgg               >  4:129498/1‑149 (MQ=255)                                                                                                                                       aTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATg           >  6:290577/1‑149 (MQ=255)                                                                                                                                           cagcaTGCCGGTCATGGACGGGCATGAACTGACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGAc                                                            >  1:30856/1‑96 (MQ=255)                                                                                                                                           cagcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGAc                                                            <  2:30856/96‑1 (MQ=255)                                                                                                                                           cagcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCTTACTCGCTGCCTGGCCAGCGGCATGGAc        <  1:202994/148‑1 (MQ=255)                                                                                                                                            ttcaTGCCCGTCCATGACCGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACg       >  2:2131/3‑148 (MQ=255)                                                                                                                                                tGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc  <  5:290577/149‑1 (MQ=255)                                                                                                                                                    ggTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGc                                     <  6:226689/110‑1 (MQ=255)                                                                                                                                                    ggTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGc                                     >  5:226689/1‑110 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGACGACTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGC  >  NC_002947/3863630‑3863921

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