Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,956,964	A→G	6.7%	Q123R (CAG→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,956,964	0	A	G	6.7%	217.3 / 3.9	74	Q123R (CAG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base G (3/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGA  >  NC_002947/5956816‑5957104                                                                                                                                                      |                                                                                                                                             gcgGGCGTATCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCa                                                                                                                                              <  4:116261/149‑1 (MQ=255) gcgGGCGTATCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCa                                                                                                                                              <  2:51240/149‑1 (MQ=255)         aTCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGa                                                                                                                                        <  4:279924/147‑1 (MQ=255)         aTCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGa                                                                                                                                        >  3:279924/1‑147 (MQ=255)            ggCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGca                                                                                                                                     >  2:336262/1‑147 (MQ=255)                     gCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATg                                                                                                                           >  6:165980/1‑148 (MQ=255)                       cAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGATCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGAc                                                                                                                        <  1:99593/149‑1 (MQ=255)                       cAAGCGCCGAGCTGGGGGCAGATGCCCAACGCCTGGCACGCGCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGAc                                                                                                                        <  8:162143/149‑1 (MQ=255)                          gcgcCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCtg                                                                                                                     >  5:309363/1‑149 (MQ=255)                              cGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCCGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGCCTC‑GCTGGCGCATGCCGAGCACAGCGACATGGACCtgtgcg                                                                                                                 >  5:364056/1‑147 (MQ=255)                              cGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGCTCCACGGCCTGTTCCTGAGGGGCAGCCTGGGCTC‑GCTGGCGCGGGCCGAGCACAGCGACATGGACCTGTggg                                                                                                                 >  3:2918/1‑149 (MQ=255)                                 gCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCGCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGgtgt                                                                                                              <  3:212084/149‑1 (MQ=255)                                      tggCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCAt                                                                                                         >  6:273776/1‑149 (MQ=255)                                          aGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCGGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAc                                                                                                     >  3:234662/1‑149 (MQ=255)                                          aGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCGGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAc                                                                                                     >  1:374938/1‑149 (MQ=255)                                          aGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAc                                                                                                     >  6:198793/1‑149 (MQ=255)                                               ctcggcGCTTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGc                                                                                                <  3:365479/144‑1 (MQ=255)                                               cccAACGCCTGGCACGCTCGTTCACCTACAAAGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGc                                                                                                >  3:219006/1‑149 (MQ=255)                                                ccAACGCCTGGCACGCTCGTTCAACTACAAGGCCAAGCATGGCAACCCGCCACGCCCGATCCACGTCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGAGACATGGACCTGTCGGTGTGCCATGCACCCGGc                                                                                                >  4:297331/1‑148 (MQ=255)                                                 cAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTgg                                                                                                                  >  1:42778/1‑129 (MQ=255)                                                 cAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTgg                                                                                                                  <  2:42778/129‑1 (MQ=255)                                                       cTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCg                                                                                         >  3:29236/1‑148 (MQ=255)                                                          gCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGaa                                                                                       <  4:361458/147‑1 (MQ=255)                                                            aCGCTCGTTCACGTTCAAGGCCGAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAgc                                                                                   <  6:364056/149‑1 (MQ=255)                                                             cGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCGGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAgct                                                                                  <  2:374938/149‑1 (MQ=255)                                                             cGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCGGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAgct                                                                                  <  4:234662/149‑1 (MQ=255)                                                                  gTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                             >  2:73631/1‑149 (MQ=255)                                                                     cACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGt                                                                           <  1:336262/148‑1 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                             <  4:120831/145‑1 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                             >  3:120831/1‑145 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGt                                                                           >  8:201969/1‑147 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGt                                                                           <  2:217429/147‑1 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGt                                                                           >  1:217429/1‑147 (MQ=255)                                                                      aCCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGt                                                                           >  2:276581/1‑147 (MQ=255)                                                                        cTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGc                                                                       <  1:236995/149‑1 (MQ=255)                                                                        cTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGc                                                                       >  7:299803/1‑149 (MQ=255)                                                                         tACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCt                                                                      >  2:4645/1‑149 (MQ=255)                                                                           cAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTgc                                                                    >  8:351584/1‑149 (MQ=255)                                                                           cAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTgc                                                                    <  6:309363/149‑1 (MQ=255)                                                                              ggCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgcc                                                                 <  1:73631/149‑1 (MQ=255)                                                                              ggCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgcc                                                                 >  7:19233/1‑149 (MQ=255)                                                                                cccAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgccg                                                                >  3:40498/1‑148 (MQ=255)                                                                                  cAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgccgc                                                               >  3:7853/1‑147 (MQ=255)                                                                                  cAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgccgc                                                               >  5:218742/1‑147 (MQ=255)                                                                                    gCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGt                                                           <  4:129630/149‑1 (MQ=255)                                                                                      aTGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGc                                                         <  6:326023/149‑1 (MQ=255)                                                                                         gCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAccc                                                                                                   <  1:355818/104‑1 (MQ=255)                                                                                         gCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAccc                                                                                                   >  2:355818/1‑104 (MQ=255)                                                                                         gCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAg                                                      <  8:299803/149‑1 (MQ=255)                                                                                           aaCCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCt                                                                      >  7:354839/1‑131 (MQ=255)                                                                                           aaCCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCt                                                                      <  8:354839/131‑1 (MQ=255)                                                                                               cGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGctgctg                                                <  1:276581/149‑1 (MQ=255)                                                                                                  cACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAg                                             <  5:165980/149‑1 (MQ=255)                                                                                                   aCGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGgtgt                                                                                                              >  1:72808/1‑84 (MQ=255)                                                                                                   aCGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGgtgt                                                                                                              <  2:72808/83‑1 (MQ=255)                                                                                                     gCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGcc                                          >  6:51671/1‑149 (MQ=255)                                                                                                              aCGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGggcggc                                  >  3:219320/1‑148 (MQ=255)                                                                                                               cGGCCTTTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCaa                                <  8:34351/149‑1 (MQ=255)                                                                                                                ggCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAg                               >  1:151028/1‑149 (MQ=255)                                                                                                                     gTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGGGAGCTGCGCCGCAAGGGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCggg                          <  1:4645/149‑1 (MQ=255)                                                                                                                         cTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGgcgc                      >  6:361344/1‑149 (MQ=255)                                                                                                                          tGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCCAGCCTGGGCGcc                     >  6:259788/1‑149 (MQ=255)                                                                                                                            aTGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGcc                     <  4:40498/147‑1 (MQ=255)                                                                                                                                gCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGaa                  >  6:192261/1‑146 (MQ=255)                                                                                                                                gCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGaa                  <  5:192261/146‑1 (MQ=255)                                                                                                                                gCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGca               >  5:261193/1‑149 (MQ=255)                                                                                                                                gCAGCCT‑GGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGca               >  8:69075/1‑148 (MQ=255)                                                                                                                                  aGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGcaca             <  5:273776/149‑1 (MQ=255)                                                                                                                                   gCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGCGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCacac            <  5:198793/149‑1 (MQ=255)                                                                                                                                    ccTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACt           <  7:201969/149‑1 (MQ=255)                                                                                                                                      tGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACttt         >  1:187094/1‑149 (MQ=255)                                                                                                                                         gCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTc      >  4:328162/1‑149 (MQ=255)                                                                                                                                         gCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTc      <  4:219006/149‑1 (MQ=255)                                                                                                                                              gCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGa  <  7:351584/149‑1 (MQ=255)                                                                                                                                                   cgcAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTc      <  2:13352/140‑1 (MQ=255)                                                                                                                                                   cgcAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTc      >  1:13352/1‑140 (MQ=255)                                                                                                                                                      |                                                                                                                                             GCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGCTGGTGGCAGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTC‑GCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGA  >  NC_002947/5956816‑5957104

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