Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,863,787	G→C	7.1%	G904R (GGG→CGG)	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,787	0	G	C	7.1%	222.4 / 5.2	70	G904R (GGG→CGG)	PP_3413	two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (3/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGC  >  NC_002947/3863640‑3863921                                                                                                                                                    |                                                                                                                                         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                                                                                                                                  >  1:216125/1‑137 (MQ=255)                  gCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                  <  2:216125/139‑3 (MQ=255)                   cTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  6:165924/1‑149 (MQ=255)                    tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  7:86986/1‑148 (MQ=255)                    tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  5:268386/1‑148 (MQ=255)                    tGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgcg                                                                                                                       >  5:299623/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                                          <  5:69536/149‑1 (MQ=255)                                                                                              cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCga                                          <  6:315464/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)                                                                                                                                             cATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                     >  6:60798/1‑123 (MQ=255)                                                                                                                                             cATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                     <  5:60798/123‑1 (MQ=255)                                                                                                                                                  aCGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                       <  1:96722/116‑1 (MQ=255)                                                                                                                                                  aCGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                       >  2:96722/1‑116 (MQ=255)                                                                                                                                                    |                                                                                                                                         CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCT‑GGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGG‑CGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGC  >  NC_002947/3863640‑3863921

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