Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,449,270	A→C	10.0%	K88Q (AAG→CAG)	nicR →	HTH‑type transcriptional repressor NicR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,449,270	0	A	C	10.0%	162.7 / 10.5	60	K88Q (AAG→CAG)	nicR	HTH‑type transcriptional repressor NicR
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base C (3/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

GCAGGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATGGGCAGCCTG  >  NC_002947/4449128‑4449413                                                                                                                                                |                                                                                                                                                gcagGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGgcgc                                                                                                                                           <  5:186077/149‑1 (MQ=255)   agGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGgcgcgc                                                                                                                                         <  3:21940/149‑1 (MQ=255)   agGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGgcgcgc                                                                                                                                         <  3:154014/149‑1 (MQ=255)    ggCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGCGCGCCTTCAGGcgctcg                                                                                                                                        >  3:269909/1‑149 (MQ=255)     gCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGgcgcgc                                                                                                                                         >  1:254951/1‑147 (MQ=255)       cTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGCGCGCCTTCAGGCGCTCGacg                                                                                                                                     <  4:269909/149‑3 (MQ=255)       cTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGc                                                                                                                                     >  6:50899/1‑149 (MQ=255)        tGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCt                                                                                                                                    >  8:192942/1‑149 (MQ=255)          cGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTgg                                                                                                                                  <  5:258491/149‑1 (MQ=255)          cGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTg                                                                                                                                   <  6:161794/148‑1 (MQ=255)           gACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCATCCGTCGCGCGCCTTCAGGCGCTCGacgatg                                                                                                                                  <  1:175223/149‑6 (MQ=255)                gCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGc                                                                                                                            >  4:74336/1‑149 (MQ=255)                       aCCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGc                                                                                                                      <  7:192942/148‑1 (MQ=255)                              tCCAGTTCGTCACCCTGTGCGCACTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGCCGCGCGACCAGGTGAAGCTTTCGCCAGACCCg                                                                                                              >  6:105680/1‑149 (MQ=255)                                  gTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCg                                                                                                          >  5:253932/1‑149 (MQ=255)                                           ccTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGgcgc                                                                                                  >  3:115692/1‑148 (MQ=255)                                               tgCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAgg                                                                                              <  3:74336/148‑1 (MQ=255)                                               tgCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGt                                                                                             >  1:116476/1‑149 (MQ=255)                                                  gcCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAgg                                                                                              >  8:68292/1‑145 (MQ=255)                                                  gcCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAgg                                                                                              <  7:68292/145‑1 (MQ=255)                                                       gcgcGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCAACGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCAtcgtc                                                                                      >  5:145836/1‑148 (MQ=255)                                                               aCGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCAcc                                                                             <  2:116476/149‑1 (MQ=255)                                                                  gCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCCCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGaa                                                                          <  4:1124/149‑1 (MQ=255)                                                                     ccAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAgc                                                                       >  7:212677/1‑149 (MQ=255)                                                                       aGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCAtc                                                                                         >  3:49522/1‑129 (MQ=255)                                                                       aGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCAtc                                                                                         <  4:49522/129‑1 (MQ=255)                                                                          tCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAgcggcgc                                                                  >  8:18465/1‑149 (MQ=255)                                                                           cccAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAgcggcgcg                                                                 >  4:157014/1‑149 (MQ=255)                                                                             cAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAgcggcgcggc                                                               >  6:138032/1‑149 (MQ=255)                                                                                       tGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGAt                                                     >  3:109704/1‑149 (MQ=255)                                                                                       tGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGAGCATGAAGGCGCGCGAGCTGGTGCAGCTTTCGCAAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGAt                                                     >  5:72429/1‑149 (MQ=255)                                                                                        gATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCAtcgtc                                                                                      >  6:35484/1‑115 (MQ=255)                                                                                        gATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCAtcgtc                                                                                      <  5:35484/115‑1 (MQ=255)                                                                                         aTCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGAt                                                     <  2:254951/147‑1 (MQ=255)                                                                                          tCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGc                                                   >  1:28255/1‑148 (MQ=255)                                                                                          tCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGc                                                   >  1:454/1‑148 (MQ=255)                                                                                            aaGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCaa                                                 <  3:157014/148‑1 (MQ=255)                                                                                            aaGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCaa                                                 <  2:454/148‑1 (MQ=255)                                                                                             aGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATg                                               <  6:253932/149‑1 (MQ=255)                                                                                               gCCACCGCGGTGGACCAGGCAACCATCCGCGGCAA‑CGTCGAGCGCCTGAAGGCCCGCGACCTGGTTCAGCATTCGCCAGACACGGGCGACCCGCCCAACGTCATCGTCGAACACACCGAAAGCGCCGCGGCCTTGCTCGATGCAATGa                                              >  2:283621/1‑148 (MQ=255)                                                                                                ccACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCAc                                                                              <  8:164035/115‑1 (MQ=255)                                                                                                ccACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCAc                                                                              >  7:164035/1‑115 (MQ=255)                                                                                                  aCCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGTAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATccc                                          <  2:62227/149‑1 (MQ=255)                                                                                                    cgcgGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCt                                        <  7:18465/149‑1 (MQ=255)                                                                                                          ggACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGccc                                  <  2:28255/149‑1 (MQ=255)                                                                                                           gACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAg                                                                                               <  7:99748/87‑1 (MQ=255)                                                                                                               aGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAg                             <  6:145836/149‑1 (MQ=255)                                                                                                                 gCATCCATCCTCGGGAT‑CGTCGAGCGCCTGAAGTCGCGCGAGCTGGTGCAGCTTTCGCCAGACGCGGGCGACCGGCGCAAGGTCATCGTCGAACTAACCGAAAGCGGAGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGAt                           <  6:72429/149‑1 (MQ=255)                                                                                                                 gCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGAt                           <  5:50899/149‑1 (MQ=255)                                                                                                                   aaCCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGc                                                          >  3:107712/1‑116 (MQ=255)                                                                                                                   aaCCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGc                                                          <  4:107712/116‑1 (MQ=255)                                                                                                                     ccATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGc                       <  4:115692/149‑1 (MQ=255)                                                                                                                       aTCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCGGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCg                      >  1:5688/1‑148 (MQ=255)                                                                                                                           gcgGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTg                 >  7:286372/1‑149 (MQ=255)                                                                                                                            cgGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACt                  >  7:122382/1‑147 (MQ=255)                                                                                                                                   cgtcgAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATgg          >  2:194788/1‑149 (MQ=255)                                                                                                                                    ctcgCGCGCCTTCAGGCGCTCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATggg         >  4:141348/2‑149 (MQ=255)                                                                                                                                         gcgcCTTCAGGCGCTCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGAt                                             <  7:253608/108‑1 (MQ=255)                                                                                                                                         gcgcCTTCAGGCGCTCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGAt                                             >  8:253608/1‑108 (MQ=255)                                                                                                                                            ccTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATGGGCAGCCTg  >  7:45501/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                GCAGGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATGGGCAGCCTG  >  NC_002947/4449128‑4449413

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