Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,704,637	T→A	11.1%	M1* (TTG→TAG) †‡	PP_1499 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,704,637	0	T	A	11.1%	152.2 / 11.4	63	M1* (TTG→TAG) †‡	PP_1499	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base A (4/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

GGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAA  >  NC_002947/1704493‑1704784                                                                                                                                                  |                                                                                                                                                    ggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGccc                                                                                                                                                 >  1:268436/1‑149 (MQ=255)      acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGCCGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGaa                                                                                                                                            <  7:283/149‑1 (MQ=255)      acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGaa                                                                                                                                            <  3:185171/149‑1 (MQ=255)       cggcgacggcgTGCTGATCGACAGCCCGGCGGAGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGGTCCAGCGCTCGGGCAACTGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCTAag                                                                                                                                           <  5:158089/149‑1 (MQ=255)        ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGaa                                                                                                                                            >  8:289078/1‑147 (MQ=255)        ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAaga                                                                                                                                          <  2:28080/149‑1 (MQ=255)        ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑AAGCCCGCGAaga                                                                                                                                          <  3:18769/149‑1 (MQ=255)        ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑AAGCCCGCGAaga                                                                                                                                          <  5:85572/149‑1 (MQ=255)        ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCCCGGGC‑TTGCCCCCGaaca                                                                                                                                          >  2:225850/1‑147 (MQ=255)                    cTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGAGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAgg                                                                                                                              >  5:24828/1‑149 (MQ=255)                     tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGt                                                                                                                             >  8:240595/1‑149 (MQ=255)                     tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGC‑AAGCCCGCGa                                                                                                                                             <  5:132830/133‑1 (MQ=255)                     tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGC‑AAGCCCGCGa                                                                                                                                             >  6:132830/1‑133 (MQ=255)                       aTCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGt                                                                                                                             <  1:200386/147‑1 (MQ=255)                          gACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCg                                                                                                                         <  7:100572/148‑1 (MQ=255)                            cAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGGTTCgtgt                                                                                                                      >  7:121985/1‑149 (MQ=255)                             aGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGAGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                     <  6:24828/149‑1 (MQ=255)                               cccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTGGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACt                                                                                                                   <  1:225850/149‑1 (MQ=255)                                        gcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAg                                                                                                          >  8:47633/1‑149 (MQ=255)                                        gcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAg                                                                                                          >  8:87483/1‑149 (MQ=255)                                                tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGt                                                                                                  >  1:36207/1‑149 (MQ=255)                                                  gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTTCCGGTgg                                                                                                <  2:190750/149‑1 (MQ=255)                                                  gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTgg                                                                                                >  8:187987/1‑149 (MQ=255)                                                        cTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAAt                                                                                          <  7:247478/149‑1 (MQ=255)                                                           gacgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATAcc                                                                                       >  3:310863/1‑149 (MQ=255)                                                            acgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                      <  3:270626/149‑1 (MQ=255)                                                                cccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGCGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCa                                                                                   >  3:235113/1‑148 (MQ=255)                                                                 ccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGGGCAGAGGCCGGGTCAGGTTTGGTGGACTAGAAAGGAGTTGCCGGTGGGCGAATAACAGCcaa                                                                                 >  4:313768/1‑149 (MQ=255)                                                                   ggCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAc                                                                                >  7:70097/1‑148 (MQ=255)                                                                      cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                              <  7:192693/131‑1 (MQ=255)                                                                      cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                              >  8:192693/1‑131 (MQ=255)                                                                       aaCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATcc                                                                           <  2:36207/149‑1 (MQ=255)                                                                               ccTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTg                                                                   >  1:67274/1‑149 (MQ=255)                                                                                cTGTTCGTCGTTGTGGCGCGGTTCCAGCGCTCGGGAAAATGAGCGTGTTGGCGTCTTGGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTg                                                                   <  8:70097/148‑1 (MQ=255)                                                                                  gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTcc                                                                >  3:216775/1‑149 (MQ=255)                                                                                  gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTCATGGCGTCTTAACGCGCATTACCCACAACAATGCC‑CTACAAGTTTCGAGTACTAGTAAGGAATTAGAGCTGGACTCATACCAGCAAACCAACAgcgcacagtac                                                                >  4:117268/1‑142 (MQ=255)                                                                                       tcgAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCg                                                           >  6:147117/1‑149 (MQ=255)                                                                                          aaGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAt                                                        <  7:187987/149‑1 (MQ=255)                                                                                           aGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtc                                                       <  7:289078/149‑1 (MQ=255)                                                                                                cgcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTcgc                                                  <  4:310863/149‑1 (MQ=255)                                                                                                 gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAACGTcgcg                                                 >  2:281057/1‑149 (MQ=255)                                                                                                  cgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                                <  2:268436/149‑1 (MQ=255)                                                                                                     gtCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑GTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGcac                                             <  7:26753/148‑1 (MQ=255)                                                                                                       ccAGCGCTCCGGCAAATGACCGTGTTGGCGTCTTCCCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTACCCGTGGACGAATACCAGCAAACCATCCGCGCCCAGTCCGACCGCAACGTCGCGGca                                              >  8:82870/1‑146 (MQ=255)                                                                                                        cAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGa                                          >  7:242007/1‑149 (MQ=255)                                                                                                           cgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGc                                       <  2:117739/149‑1 (MQ=255)                                                                                                             cTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑ATGCCCGCCAAGAGGCCCGTACAGGTTTCGAGTACTAGAAAGGAGTTCCCGGTGGACCAATACCAGCAAACCATACCCGCCCAGTCCGACCCCAACCGCGCGGCACACACGCCg                                     >  4:240089/1‑149 (MQ=255)                                                                                                                  gCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATcc                                 <  1:122436/148‑1 (MQ=255)                                                                                                                  gCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGGACAGGTTTAGTGTACTTGAAAGGAGGTCCCGGTGGACGAATACCACCAAACCAAACGCGCCCAGTACGACCGCATCGTAGCGGCACAAACGCCGATcc                                 >  2:122436/1‑148 (MQ=255)                                                                                                                   cAAATGTCCTTGTTGGCTTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTAGAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAGCATCCGCGCCCTGTCGGACCGAATCGTCGCGGCACAGACTCCGATCCgg                               <  3:117268/149‑1 (MQ=255)                                                                                                                   cAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCgg                               <  8:121985/149‑1 (MQ=255)                                                                                                                    aaaTGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATc                                  <  2:241425/145‑1 (MQ=255)                                                                                                                    aaaTGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCCGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATc                                  >  1:241425/1‑145 (MQ=255)                                                                                                                               ttGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgcg                   <  5:147117/149‑1 (MQ=255)                                                                                                                               ttGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgc                    <  7:240595/148‑1 (MQ=255)                                                                                                                               ttGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTACGACCGCAACGTCGCGGCACAGACGCCGATCCGGGTCCATGAcgcg                   >  6:293553/1‑149 (MQ=255)                                                                                                                                 ggCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAACGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCgg                  <  1:281057/148‑1 (MQ=255)                                                                                                                                  gcctCTTCGCGGGC‑AAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt                 >  3:98973/4‑148 (MQ=255)                                                                                                                                      cTTCGCGGGC‑AAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt            >  7:144271/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGC‑AAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg           >  5:13564/1‑149 (MQ=255)                                                                                                                                           cgGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac       <  7:47633/149‑1 (MQ=255)                                                                                                                                            gggC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacg      <  2:67274/149‑1 (MQ=255)                                                                                                                                              gC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacgac    <  8:242007/149‑1 (MQ=255)                                                                                                                                                 ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGacaa  >  8:245725/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGC‑TTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAA  >  NC_002947/1704493‑1704784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑