Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,417,603	C→G	6.7%	R166G (CGC→GGC)	PP_3905 →	LuxR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,417,603	0	C	G	6.7%	186.2 / 3.8	60	R166G (CGC→GGC)	PP_3905	LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base G (2/2);  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.50e-01

GTTGCCCGGCGCGCCGCGATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAA  >  NC_002947/4417456‑4417741                                                                                                                                                     |                                                                                                                                           gTTGCCCGGCGCGCCGCGATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGcgccg                                                                                                                                           >  2:153946/1‑149 (MQ=255)  ttGCCCGGCGCGCCGCGATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGcgccgc                                                                                                                                          <  8:45912/149‑1 (MQ=255)                 cgATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGCCCCCATGCGCCGCGATCGGCGCGCTCC‑GGAGCGCCGCACGGCCCAAGGTCg                                                                                                                            >  4:247943/1‑148 (MQ=255)                  gATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGaa                                                                                                                          <  8:70847/149‑1 (MQ=255)                         gCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTc                                                                                                                    <  4:108007/148‑1 (MQ=255)                         gCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCg                                                                                                                   >  7:71989/1‑149 (MQ=255)                               gCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGc                                                                                                             >  2:119873/1‑149 (MQ=255)                                 aaTGCGTCGCCACATCATCAGCCCTCTCTCCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAACATCGAATCGTTAGCAAGGcc                                                                                                            >  8:23108/1‑148 (MQ=37)                                   tGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATACGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAg                                                                                                         <  2:96105/149‑1 (MQ=255)                                         gCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTc                                                                                                   <  2:285381/149‑1 (MQ=255)                                              atcatcAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATa                                                                                              >  3:131212/1‑149 (MQ=255)                                              atcatcAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATa                                                                                              >  1:40652/1‑149 (MQ=255)                                               tcatcaGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATa                                                                                              <  8:71989/148‑1 (MQ=255)                                                 atcaGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCa                                                                                           >  2:14391/1‑149 (MQ=255)                                                  tcaGCACTCTCTGCATCCTGCTTGCCACCCTCATCGCCATGCATGAGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCTCAAGGCCGAGGCCTTCGCATACAAc                                                                                          >  5:295664/1‑149 (MQ=21)                                                    aGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGc                                                                                        <  3:213485/149‑1 (MQ=255)                                                      cccTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCt                                                                                      >  7:92025/1‑149 (MQ=255)                                                      cccTCTCTGCATCCTGCTTGCAAGCCTCATCGCCATGCATGCGGTGATCGCCGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAAGCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCt                                                                                      >  6:269983/1‑149 (MQ=255)                                                        ctctctGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCg                                                                                                                   <  5:309670/118‑1 (MQ=255)                                                        ctctctGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCg                                                                                                                   >  6:309670/1‑118 (MQ=255)                                                            ctGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGaccca                                                                                >  3:67860/1‑149 (MQ=255)                                                               cATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGacccaccc                                                                             >  1:316963/1‑149 (MQ=255)                                                                aTCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGtgccgc                                                                                                                                          <  4:283916/87‑1 (MQ=255)                                                                aTCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGtgccgc                                                                                                                                          >  3:283916/1‑87 (MQ=255)                                                                    tGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGa                                                                        >  1:240393/1‑149 (MQ=255)                                                                    tGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAAAACCACGCCTGACCCACCCGAgg                                                                         >  7:10069/1‑148 (MQ=255)                                                                     gCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGAt                                                                       >  1:224089/1‑149 (MQ=255)                                                                     gCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGAt                                                                       >  1:229148/1‑149 (MQ=255)                                                                      cTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGa                                                                        <  2:316963/147‑1 (MQ=255)                                                                      cTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGAt                                                                       <  2:240393/148‑1 (MQ=255)                                                                            aGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATa                                                                                              >  6:39930/1‑119 (MQ=255)                                                                            aGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATa                                                                                              <  5:39930/119‑1 (MQ=255)                                                                              ccTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATg                                                              <  7:23108/149‑1 (MQ=255)                                                                                  aTCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCa                                                            <  2:361753/147‑1 (MQ=255)                                                                                   tCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTc                                                                                                                    <  1:217089/90‑1 (MQ=255)                                                                                   tCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTc                                                                                                                    >  2:217089/1‑90 (MQ=255)                                                                                      ccATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGAt                                                                       <  7:315438/132‑1 (MQ=255)                                                                                      ccATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGAt                                                                       >  8:315438/1‑132 (MQ=255)                                                                                              gcGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCAt                                               >  2:87060/1‑148 (MQ=255)                                                                                              gcGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCAt                                               >  6:298800/1‑148 (MQ=255)                                                                                                ggTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATaaa                                            <  1:153946/149‑1 (MQ=255)                                                                                                ggTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATaaa                                            <  1:14391/149‑1 (MQ=255)                                                                                                 gTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCCGGGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCCCGCCTGACCCACCCGAGGATCCCCACATGCAGACAGCCGTGCATa                                              >  1:74550/1‑147 (MQ=255)                                                                                                     tCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCttt                                       <  4:131212/149‑1 (MQ=255)                                                                                                     tCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACAACGCCTGACCCACCCGATGATCACCACATGCAGACAGCCGTGCATAAAGCttt                                       <  2:40652/149‑1 (MQ=255)                                                                                                        acgacAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAg                                           <  4:32490/142‑1 (MQ=255)                                                                                                        acgacAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAg                                           >  3:32490/1‑142 (MQ=255)                                                                                                          gacAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCg                                     <  5:64330/146‑1 (MQ=255)                                                                                                          gacAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCg                                     >  6:64330/1‑146 (MQ=255)                                                                                                               cGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAg                             >  4:238131/1‑149 (MQ=255)                                                                                                                gACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGc                            <  8:92025/149‑1 (MQ=255)                                                                                                                     cATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATg                       <  2:229148/149‑1 (MQ=255)                                                                                                                     cATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATg                       <  2:224089/149‑1 (MQ=255)                                                                                                                           ccgcGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTc                 <  8:265809/149‑1 (MQ=255)                                                                                                                             gcgATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCg                >  2:98123/1‑148 (MQ=255)                                                                                                                                  tgcggcgcTCC‑GGAGCGCGGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGt          >  2:363489/5‑149 (MQ=255)                                                                                                                                  cgccgcGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGt          >  4:109874/1‑149 (MQ=255)                                                                                                                                     cgcgCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGc       <  3:238131/149‑1 (MQ=255)                                                                                                                                     cgcgCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGc       <  4:67860/149‑1 (MQ=255)                                                                                                                                         cTCC‑GGAGCGCGGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGa   <  1:363489/149‑1 (MQ=255)                                                                                                                                         cTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGa   <  3:109874/149‑1 (MQ=255)                                                                                                                                          tCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGaa  >  5:8709/1‑149 (MQ=255)                                                                                                                                            c‑GGAGCGCGGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATcac                                                                    <  7:151677/81‑1 (MQ=255)                                                                                                                                            c‑GGAGCGCGGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATcac                                                                    >  8:151677/1‑81 (MQ=255)                                                                                                                                                     |                                                                                                                                           GTTGCCCGGCGCGCCGCGATGGTCGCCGAGGCAATGCGTCGCCAGATCATCAGCCCTCTCTGCATCCTGCTTGCCAGCCTCATCGCCATGCATGCGGTGATCGACGACAGCGACCCCATGCGCCGCGATCGCCGCGCTCC‑GGAGCGCCGCACCGCCCAAGTTCGAATCGTTCGCAAGGCCGAGGCCTTCGAATACCACGCCTGACCCACCCGAGGATCACCACATGCAGACAGCCGTGCATAAAGCTTTCGAGGACAAGCGAATGGTGCTCGCAGCGTTGCTTGAA  >  NC_002947/4417456‑4417741

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