Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,449,269	G→T	10.0%	L87L (CTG→CTT)	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,269	0	G	T	10.0%	170.0 / 8.3	60	L87L (CTG→CTT)	nicR	HTH‑type transcriptional repressor NicR
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base T (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 = 4.35e-01

CCAGCAGCAGGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATGGGCAGCCTG  >  NC_002947/4449122‑4449413                                                                                                                                                     |                                                                                                                                                 ccAGCAGCAGGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTg                                                                                                                                                  >  5:161794/1‑148 (MQ=255)       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                                                                                              <  7:68292/145‑1 (MQ=255)                                                        gcCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAgg                                                                                              >  8:68292/1‑145 (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                                                 <  2:454/148‑1 (MQ=255)                                                                                                  aaGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCaa                                                 <  3:157014/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                                                          <  4:107712/116‑1 (MQ=255)                                                                                                                         aaCCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGc                                                          >  3:107712/1‑116 (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)                                                                                                                                                     |                                                                                                                                                 CCAGCAGCAGGCCTGCGACCCGCAGCTTACCTCGATCCAGTTCGTCACCCTGTGCGCCCTGCGCGATCACGGCCCCAGCTCCCAGGCTGAACTGATCAAGGCCACCGCGGTGGACCAGGCAACCATCCGCGGCAT‑CGTCGAGCGCCTGAAGGCGCGCGAGCTGGTGCAGCTTTCGCCAGACCCGGGCGACCGGCGCAAGGTCATCGTCGAACTCACCGAAAGCGGCGCGGCCTTGCTCGATGCAATGATCCCCTGCGCCCGGCAGATCAGCGAACTGAGCATGGGCAGCCTG  >  NC_002947/4449122‑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: ‑