Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,586,285	C→G	6.2%	*61S (TGA→TCA)	PP_1390 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,586,285	0	C	G	6.2%	206.9 / 3.8	65	*61S (TGA→TCA)	PP_1390	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CCTGCGGCTGCAGCGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGA  >  NC_002947/1586138‑1586426                                                                                                                                                    |                                                                                                                                               cctgcGGCTGCAGCGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGttcattc                                                                                                                                                >  5:189381/1‑148 (MQ=255) cctgcGGCTGCAGCGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGttcattc                                                                                                                                                >  6:296738/1‑148 (MQ=255) cctgcGGCTGCAGCGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGttcattc                                                                                                                                                >  2:3953/1‑148 (MQ=255)              cGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTCCGTTGAATGAACGgcgccc                                                                                                                                  <  2:237413/149‑3 (MQ=255)               gTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCa                                                                                                                                 <  5:296738/149‑1 (MQ=255)                     ggCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCTCTGGGCCAAGGAATACATGAGAGTCGAACAGTGATCGGCAGGGCTCCGTTAATACAACGGACGGCTCTAAAG‑Gtag                                                                                                                          >  4:214822/1‑147 (MQ=255)                     ggCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTg                                                                                                                          >  7:162899/1‑149 (MQ=255)                      tcataaCAAGTCGTAAAACCTCAGTGCCACCGAGCTCTAGCACAACTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTTATCGGCAGGGCTCCGTTCATTCAACGGACGTCGCTCAAG‑GTTGc                                                                                                                         <  4:259259/145‑1 (MQ=255)                      gCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTg                                                                                                                          >  3:259259/1‑148 (MQ=255)                       cAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCa                                                                                                                        >  7:276235/1‑149 (MQ=255)                       cAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCa                                                                                                                        <  4:293200/149‑1 (MQ=255)                         gAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACt                                                                                                                      <  4:328464/149‑1 (MQ=255)                            cAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACtgtg                                                                                                                   <  5:101748/149‑1 (MQ=255)                                  ggAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                             >  5:64298/1‑149 (MQ=255)                                       aCCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTTGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAt                                                                                                        <  7:185980/149‑1 (MQ=255)                                         cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                             <  5:342407/142‑1 (MQ=255)                                         cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                             >  6:342407/1‑142 (MQ=255)                                         cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgt                                                                                                      >  2:356296/1‑149 (MQ=255)                                            aGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAGTGTTGCACTGTGAGGCACGCAATGTGc                                                                                                    <  1:3953/149‑1 (MQ=255)                                                        ctctCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTa                                                                                       <  5:108223/149‑1 (MQ=255)                                                         tctcGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTAt                                                                                      >  3:160968/1‑149 (MQ=255)                                                             gCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGAGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTgg                                                                                   <  3:214822/148‑1 (MQ=255)                                                                     acacCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACgctgg                                                                          >  4:377654/1‑149 (MQ=255)                                                                       acCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACgctgggc                                                                        >  8:147824/1‑149 (MQ=255)                                                                           agccagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGTCGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGggctggc                                                                    <  8:238305/149‑1 (MQ=255)                                                                           agccagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGggctggc                                                                    >  6:182746/1‑149 (MQ=255)                                                                              cagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GGTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGc                                                                 <  7:216265/149‑1 (MQ=255)                                                                               agaagaCTACATTGGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCg                                                                <  8:276235/149‑1 (MQ=255)                                                                                gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgc                                                                                                      <  6:140248/110‑2 (MQ=255)                                                                                gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgc                                                                                                      >  5:140248/1‑109 (MQ=255)                                                                                gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                               <  8:162899/149‑1 (MQ=255)                                                                                  agaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAg                                                             <  3:270743/149‑1 (MQ=255)                                                                                    aCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                           <  7:147824/149‑1 (MQ=255)                                                                                     cTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                           <  6:64298/148‑1 (MQ=255)                                                                                      tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                         >  8:337331/1‑149 (MQ=255)                                                                                      tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                         >  7:161985/1‑149 (MQ=255)                                                                                      tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                         >  8:232750/1‑149 (MQ=255)                                                                                       aCATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGtt                                                        >  2:220715/1‑149 (MQ=255)                                                                                        cATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAacgga                                                                                                                                          >  4:191659/1‑67 (MQ=255)                                                                                        cATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAacgga                                                                                                                                          <  3:191659/67‑1 (MQ=255)                                                                                         aTTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCg                                                      <  1:356296/149‑1 (MQ=255)                                                                                            cGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGgc                                                    <  5:182746/148‑1 (MQ=255)                                                                                             gCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGgcgc                                                  >  1:79295/1‑149 (MQ=255)                                                                                                   ccAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCcag                                             >  3:302441/1‑148 (MQ=255)                                                                                                     aaGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCcagcag                                          >  2:201588/1‑149 (MQ=255)                                                                                                       ggAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGc                                        >  6:339743/1‑149 (MQ=255)                                                                                                        gAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCt                                       >  6:265698/1‑149 (MQ=255)                                                                                                              aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     <  1:16398/145‑1 (MQ=255)                                                                                                              aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     >  2:16398/1‑145 (MQ=255)                                                                                                              aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     <  1:22586/145‑1 (MQ=255)                                                                                                              aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     >  2:22586/1‑145 (MQ=255)                                                                                                                       gAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCt                        >  3:328554/1‑149 (MQ=255)                                                                                                                       gAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGACTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCt                        >  1:48452/1‑149 (MQ=255)                                                                                                                         agagTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGc                      <  4:160968/149‑1 (MQ=255)                                                                                                                               aTCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     <  6:113000/128‑1 (MQ=255)                                                                                                                               aTCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                     >  5:113000/1‑128 (MQ=255)                                                                                                                                  ggCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCgg              >  4:302739/1‑148 (MQ=255)                                                                                                                                   gCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGt             >  7:112520/1‑148 (MQ=255)                                                                                                                                    ttgagcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGt             >  8:122241/5‑147 (MQ=255)                                                                                                                                     aGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                           <  5:248905/100‑1 (MQ=255)                                                                                                                                     aGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                           >  6:248905/1‑100 (MQ=255)                                                                                                                                      gggCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAgg         >  5:1726/1‑149 (MQ=255)                                                                                                                                       ggCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAggg        <  6:189381/149‑1 (MQ=255)                                                                                                                                       agcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                               >  2:147759/2‑94 (MQ=38)                                                                                                                                       agcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                               <  1:147759/93‑1 (MQ=38)                                                                                                                                           ccGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCAt    <  7:337331/149‑1 (MQ=255)                                                                                                                                           ccGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCAt    >  8:291369/1‑149 (MQ=255)                                                                                                                                             gTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGa  <  2:79295/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               CCTGCGGCTGCAGCGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGA  >  NC_002947/1586138‑1586426

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