Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	150,817	A→T	5.5%	C45S (TGC→AGC)	PP_0142 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	150,817	0	A	T	5.5%	222.6 / 2.8	72	C45S (TGC→AGC)	PP_0142	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base T (2/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GTTTTTCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTC  >  NC_002947/150671‑150960                                                                                                                                                   |                                                                                                                                                  gTTTTTCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAAt                                                                                                                                                 >  1:168301/1‑149 (MQ=255) gTTTTTCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAAt                                                                                                                                                 >  7:10422/1‑149 (MQ=255) gTTTTTCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAAt                                                                                                                                                 >  6:68000/1‑149 (MQ=255)      tCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAg                                                                                                                                            >  6:303168/1‑149 (MQ=255)          cagcGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAgc‑g                                                                                                                                         >  1:59551/1‑147 (MQ=255)           agcGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAg                                                                                                                                      <  5:68000/149‑1 (MQ=255)              gcCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGgcg                                                                                                                                   >  4:305322/1‑149 (MQ=255)                   gctggctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGcgctc                                                                                                                                <  5:108139/147‑1 (MQ=255)                     tggctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGcgctcgct                                                                                                                             >  6:309302/1‑148 (MQ=255)                     tggctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCgctcgctc                                                                                                                            >  6:321418/1‑149 (MQ=255)                         ttCAGGCAGCGCGTGGGTGAAGTGCGACAGGCGCTCGCAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCTGCAGTGTCCAGGCGGACCAGCTGGCAGACGTCGGCGACGGTCTCGGCAGGGCAATGCGAGC‑CCAGCCGCTCGATCTcga                                                                                                                        >  8:94663/1‑148 (MQ=255)                            aGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCgcgc                                                                                                                      >  5:151423/1‑148 (MQ=255)                            aGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCgcgc                                                                                                                      <  6:151423/148‑1 (MQ=255)                            aGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCt                                                                                                                     <  1:39240/149‑1 (MQ=255)                             ggCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCtt                                                                                                                    <  2:199255/149‑1 (MQ=255)                               cAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGa                                                                                                                  <  5:2786/149‑1 (MQ=255)                                aGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAg                                                                                                                 <  8:161198/149‑1 (MQ=255)                                  cgcgTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGt                                                                                                               <  3:297250/149‑1 (MQ=255)                                     gTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCt                                                                                                                     >  3:199858/1‑140 (MQ=255)                                     gTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCt                                                                                                                     <  4:199858/140‑1 (MQ=255)                                     gTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTgg                                                                                                            >  2:33881/1‑149 (MQ=255)                                             aGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGCGTCCAGGCAGCCCAGCTGGCTGAGGTCGGCGGCGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAg                                                                                                       <  7:62248/146‑1 (MQ=255)                                             aGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGAGTCCAGGCAGCCCAGCTGGCTGAGGTCGGCGGCGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAg                                                                                                       >  8:62248/1‑146 (MQ=255)                                              gTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGtgtg                                                                                                   >  5:94595/1‑149 (MQ=255)                                                 cGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCTATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGc                                                                                                  >  6:152666/1‑147 (MQ=255)                                                 cGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCTATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGc                                                                                                  >  6:148082/1‑147 (MQ=255)                                                 cGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑a                                                                                               <  6:94595/149‑1 (MQ=255)                                                   aCAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑Agg                                                                                             <  5:309302/149‑1 (MQ=255)                                                           tCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGCCCCAGCTGGCTGAGGCCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCGTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTc                                                                                     >  7:151104/1‑149 (MQ=255)                                                           tCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCCGCAGTGTCCAGGCGGCCCAGCAGGCTGAGGTCGGCGACGCAGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCATGAGGTTGGCGTAGTGTCCC‑AAGGTCCAGt                                                                                      >  7:134672/1‑148 (MQ=255)                                                             ggAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCAc                                                                                   >  5:75471/1‑149 (MQ=255)                                                                  cTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGt                                                                               >  7:188806/1‑148 (MQ=255)                                                                  cTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGt                                                                               >  5:64974/1‑148 (MQ=255)                                                                   tAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTaa                                                                             >  1:151456/1‑149 (MQ=255)                                                                   tAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTaa                                                                             >  5:1999/1‑149 (MQ=255)                                                                     gCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATa                                                                           >  7:42774/1‑149 (MQ=255)                                                                           cGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACgcag                                                                      <  5:321418/148‑1 (MQ=255)                                                                           cGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACgcag                                                                      <  2:238824/148‑1 (MQ=255)                                                                               cagcagCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGCCGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTGGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGcaggcagg                                                                 >  5:57849/1‑149 (MQ=255)                                                                               cagcagCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGcgctcg                                                                                                                               <  3:167222/88‑1 (MQ=255)                                                                               cagcagCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGcgctcg                                                                                                                               >  4:167222/1‑88 (MQ=255)                                                                                 gcagcGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAAGGGTCCAGTCACCGGTAATACGCaggcaggc                                                                <  1:33881/149‑1 (MQ=255)                                                                                  cagcGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCg                                                              <  6:1999/149‑1 (MQ=255)                                                                                     cGAAGCGCCGGCCGGGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGtt                                                           <  6:57849/149‑1 (MQ=255)                                                                                     cGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGtt                                                           <  1:92201/149‑1 (MQ=255)                                                                                          cgcCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGgctg                                                      <  2:59551/149‑1 (MQ=255)                                                                                            ccGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGc                                                                                                           >  2:311092/1‑95 (MQ=255)                                                                                             cGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑TGGGTCCAGTCACCGGTAATACGCTGGCAGGCCGGTTGGctgctg                                                   <  8:10422/149‑1 (MQ=255)                                                                                             cGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGctgctg                                                   <  6:75471/149‑1 (MQ=255)                                                                                             cGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGctgctg                                                   <  6:64974/149‑1 (MQ=255)                                                                                               gCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGgt                                                 >  8:160996/1‑149 (MQ=255)                                                                                               gCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGgt                                                 <  7:94663/149‑1 (MQ=255)                                                                                                  gtgtCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTc                                              <  8:42774/149‑1 (MQ=255)                                                                                                     tCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCa                                            <  2:151456/148‑1 (MQ=255)                                                                                                      ccAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAgg                                          >  3:91450/1‑149 (MQ=255)                                                                                                       cAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAgg                                          <  7:140987/148‑1 (MQ=255)                                                                                                       cAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAgg                                          <  5:303168/148‑1 (MQ=255)                                                                                                               ccAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTg                                 <  3:305322/149‑1 (MQ=255)                                                                                                                cAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTgg                                <  4:160519/149‑1 (MQ=255)                                                                                                                 aGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTggg                               <  7:160996/149‑1 (MQ=255)                                                                                                                      gctgAGGTCGGCGACGGTGTCGGCAGGGCTATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTg                           <  5:148082/148‑1 (MQ=255)                                                                                                                      gctgAGGTCGGCGACGGTGTCGGCAGGGCTATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTg                           <  5:152666/148‑1 (MQ=255)                                                                                                                      gctgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTgg                          <  7:264902/149‑1 (MQ=255)                                                                                                                      gctgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTgg                          <  2:168301/149‑1 (MQ=255)                                                                                                                        tgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGt                         >  8:137882/1‑148 (MQ=255)                                                                                                                           ggTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTc                        <  5:211068/146‑1 (MQ=255)                                                                                                                           ggTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTc                        >  6:211068/1‑146 (MQ=255)                                                                                                                               ggCGACGGTGTCGGCAGGGCAATGCGAGCGGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAgg                   >  1:24230/1‑148 (MQ=255)                                                                                                                                gCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTc                <  7:137882/149‑1 (MQ=255)                                                                                                                                gCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTc                >  2:212356/1‑149 (MQ=255)                                                                                                                                              aGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTc  <  2:24230/149‑1 (MQ=255)                                                                                                                                              aGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTc  >  5:319593/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTTTTTCAGCAGCGCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGC‑GCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCC‑AGGGTCCAGTCACCGGTAATACGCAGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTC  >  NC_002947/150671‑150960

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