Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,159,062	G→C	8.8%	A361P (GCA→CCA)	opuA →	glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,159,062	0	G	C	8.8%	205.3 / 8.7	68	A361P (GCA→CCA)	opuA	glycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base C (3/3);  total (31/37)
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.31e-01

TCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAG  >  NC_002947/3158914‑3159202                                                                                                                                                     |                                                                                                                                             tCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                                                                                                                                              >  4:160408/1‑149 (MQ=255)   gCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                                                                                                                                              >  2:186675/1‑147 (MQ=255)   gCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                                                                                                                                              >  1:200517/1‑147 (MQ=255)    ccAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAg                                                                                                                                           >  7:98299/1‑149 (MQ=255)     cAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGGGTGGCGCGAGGGGGCGGGGGGGCTATCGGCCACAGCACTGGCAGt                                                                                                                                          <  1:218925/149‑1 (MQ=255)     cAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGt                                                                                                                                          >  4:144156/1‑149 (MQ=255)           gtgtTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                    <  5:335299/149‑3 (MQ=255)                 ggCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCaa                                                                                                                              >  6:3376/1‑149 (MQ=255)                  gCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAAt                                                                                                                             <  4:303747/149‑1 (MQ=255)                  gCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                    <  5:37442/142‑3 (MQ=255)                  gCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                    >  6:37442/1‑140 (MQ=255)                   cccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                      <  1:194904/139‑1 (MQ=255)                   cccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                      >  2:194904/1‑139 (MQ=255)                   cccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATc                                                                                                                            >  1:49195/1‑149 (MQ=255)                    ccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCg                                                                                                                           >  8:275185/1‑149 (MQ=255)                    ccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                    <  5:65256/140‑3 (MQ=255)                    ccGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                    >  6:65256/1‑138 (MQ=255)                          gcgcTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTggg                                                                                                                     <  5:231897/149‑1 (MQ=255)                            gctcgctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                   >  7:161821/1‑149 (MQ=255)                                gctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAgtggtgtg                                                                                                                                    >  1:194435/1‑126 (MQ=255)                                 ctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGAGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCt                                                                                                              >  1:249540/1‑149 (MQ=255)                                    ccgccgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           <  8:98299/149‑1 (MQ=255)                                    ccgccgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           <  2:31833/149‑1 (MQ=255)                                    ccgccgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGGGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGGAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           <  4:66889/149‑1 (MQ=255)                                     cgccgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCAAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTggg                                                                                                          >  4:133337/1‑149 (MQ=255)                                          gACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGc                                                                                                     <  2:329039/149‑1 (MQ=255)                                             gAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTg                                                                                                  <  1:186675/149‑1 (MQ=255)                                                  ggttggttGTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           <  7:67044/135‑1 (MQ=255)                                                  ggttggttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           >  8:67044/1‑135 (MQ=255)                                                   gttggttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCtt                                                                                                                        <  4:178846/121‑1 (MQ=255)                                                   gttggttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCtt                                                                                                                        >  3:178846/1‑121 (MQ=255)                                                   gttggttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGAAGTGCTGCTGGTCAATCGCttcg                                                                                                                      >  4:217917/1‑121 (MQ=255)                                                       gttCTGGCTGGCGCGGCGGGGGGAGGGGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggccac                                                                                         <  2:200517/148‑1 (MQ=255)                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTCGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggccacgg                                                                                       >  2:314982/1‑149 (MQ=255)                                                            ggctggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGt                                                                                   >  4:15015/1‑149 (MQ=255)                                                             gctggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGt                                                                                   <  8:161821/148‑1 (MQ=255)                                                              ctggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGa                                                                                 >  7:190287/1‑149 (MQ=255)                                                                ggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGgccgtggccgtgcag                                                                               <  3:270839/149‑4 (MQ=255)                                                                 gcgcTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACt                                                                               >  8:98886/1‑148 (MQ=255)                                                                  cgcTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCa                                                                             <  2:49195/149‑1 (MQ=255)                                                                             gAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggccacgg                                                                                       >  2:129127/1‑128 (MQ=255)                                                                             gAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggccacgg                                                                                       <  1:129127/128‑1 (MQ=255)                                                                              aGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCt                                                                 <  2:249540/149‑1 (MQ=255)                                                                                gTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCtt                                                                >  1:197300/1‑148 (MQ=255)                                                                                  ggcggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTAc                                                             >  1:276602/1‑149 (MQ=255)                                                                                     ggcCCTGGGCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCt                                                          >  8:109904/1‑149 (MQ=255)                                                                                     ggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCt                                                          >  3:11470/1‑149 (MQ=255)                                                                                        ccTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTgg                                                       <  1:314982/149‑1 (MQ=255)                                                                                        ccTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTgg                                                       <  4:20884/149‑1 (MQ=255)                                                                                        ccTGGCCATGCTGGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTgg                                                       <  7:275185/149‑1 (MQ=255)                                                                                          tGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGgcg                                                     >  5:345160/1‑149 (MQ=255)                                                                                            gCCATGGTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCg                                                                   <  6:240970/133‑1 (MQ=255)                                                                                            gCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgct                                                   <  8:190287/149‑1 (MQ=255)                                                                                            gCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgct                                                   <  3:15015/149‑1 (MQ=255)                                                                                            gCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGATGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCg                                                                   >  5:240970/1‑133 (MQ=255)                                                                                              cATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           <  2:9206/91‑1 (MQ=255)                                                                                              cATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                           >  1:9206/1‑91 (MQ=255)                                                                                                  cttgctGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTAc                                              <  2:276602/148‑1 (MQ=255)                                                                                                      ctGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGt                                          <  5:3376/148‑1 (MQ=255)                                                                                                       tggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGa                                        <  7:98886/149‑1 (MQ=255)                                                                                                        ggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctg                                                  <  3:188452/138‑1 (MQ=255)                                                                                                        ggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctg                                                  >  4:188452/1‑138 (MQ=255)                                                                                                           gcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCa                                    <  3:144156/149‑1 (MQ=255)                                                                                                               ggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCa                                >  5:279203/1‑149 (MQ=255)                                                                                                                   caaGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGa                            <  3:133337/147‑1 (MQ=255)                                                                                                                          gCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATg                     >  5:17123/1‑149 (MQ=255)                                                                                                                           cGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATgg                    <  2:326221/149‑1 (MQ=255)                                                                                                                             gggTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctgc                                                 >  2:111733/1‑118 (MQ=255)                                                                                                                             gggTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctgc                                                 <  1:111733/118‑1 (MQ=255)                                                                                                                                      tCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGAc         <  4:11470/149‑1 (MQ=255)                                                                                                                                             caGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCagag  <  2:197300/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             TCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAG  >  NC_002947/3158914‑3159202

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