Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,069	G→T	6.1%	L627M (CTG→ATG) ‡	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,069	0	G	T	6.1%	259.5 / 6.0	82	L627M (CTG→ATG) ‡	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base T (3/2);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

TTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTCT  >  NC_002947/421922‑422216                                                                                                                                                     |                                                                                                                                                    ttGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAgc                                                                                                                                                    <  3:389078/149‑1 (MQ=255)     cATGCGGCCAATCGACAGCACCAGCAAAGCGCGGATGGGCAGAAGCAGCCTGATGAAATCCAGCAAACGCCAGCTCAGAGCCATGACCGACTCCGCAGTCACTTCCATGGACGTAATGCTGTCGGCGTCGGACTTGGCCA‑TCAGAGccc                                                                                                                                                >  6:226080/1‑149 (MQ=255)       tGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCgg                                                                                                                                              <  5:99910/149‑1 (MQ=255)         cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGc                                                                                                                                            >  8:371890/1‑149 (MQ=255)         cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGC‑TGATCGCCCagg                                                                                                                                             <  5:363951/148‑1 (MQ=255)          ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGcc                                                                                                                                           >  6:338896/1‑149 (MQ=255)            ccAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTc                                                                                                                                         >  1:261171/1‑149 (MQ=255)             cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCg                                                                                                                                        >  2:330802/1‑149 (MQ=255)              actCGACAGCACCCGCAACGCGCGGCTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                       <  8:300397/147‑1 (MQ=255)              aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                       >  2:76212/1‑149 (MQ=255)              aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGGCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                       <  5:48294/149‑1 (MQ=255)                     aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                                <  1:213276/149‑1 (MQ=255)                     aGCACCAGCAACGCGCGGGTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGGTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                                <  6:209178/149‑1 (MQ=255)                       cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                               >  6:170506/1‑148 (MQ=255)                       cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                               >  5:77799/1‑148 (MQ=255)                            gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCg                                                                                                                         <  2:261171/149‑1 (MQ=255)                                  cgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaa                                                                                                                   <  7:371890/149‑1 (MQ=255)                                   gcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                                    >  6:357459/1‑147 (MQ=255)                                     ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                                >  2:279675/1‑149 (MQ=255)                                     ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                                <  5:338896/149‑1 (MQ=255)                                     ggTTGGGCAGCAGCAGGATGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCTCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCTA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                                <  5:259167/149‑1 (MQ=255)                                           gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                          >  8:256229/1‑149 (MQ=255)                                             agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                                        <  6:39341/149‑1 (MQ=255)                                                      cTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCggg                                                                                               <  5:170506/149‑1 (MQ=255)                                                        gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                              <  4:359345/132‑1 (MQ=255)                                                        gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                              >  3:359345/1‑132 (MQ=255)                                                         atgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGt                                                                                            <  5:357459/149‑1 (MQ=255)                                                          tgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                           >  7:29160/1‑149 (MQ=255)                                                            aaaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTgg                                                                                          <  7:528/148‑1 (MQ=255)                                                             aaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGt                                                                                         >  8:365945/1‑148 (MQ=255)                                                              aTCCAGCAAACGCCTGCTCAGCGCCATGACCGACTCCGCATTCACTTCCGTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCTAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                                        <  5:226080/148‑1 (MQ=255)                                                              aTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                                        >  5:378088/1‑148 (MQ=255)                                                                ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                                  <  7:30933/120‑1 (MQ=255)                                                                ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                                  >  8:30933/1‑120 (MQ=255)                                                                 cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                                    >  7:313245/1‑149 (MQ=255)                                                                 cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                                    >  4:329098/1‑149 (MQ=255)                                                                 cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                                     <  1:76212/148‑1 (MQ=255)                                                                  aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                                     >  6:390894/1‑147 (MQ=255)                                                                  aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                                     <  5:390894/147‑1 (MQ=255)                                                                   gCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGc                                                                                  <  5:219926/149‑1 (MQ=255)                                                                        cGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGc                                                                              >  1:274771/1‑148 (MQ=255)                                                                          ccAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCa                                                                           <  2:274771/149‑1 (MQ=255)                                                                               tCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                                      >  7:88021/1‑149 (MQ=255)                                                                                 aGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                                      <  1:330802/147‑1 (MQ=255)                                                                                  tcgcCATGACCGACTCCGCATTCACTTCGCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                                   <  8:346828/148‑1 (MQ=255)                                                                                     ccATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTg                                                                <  7:256229/149‑1 (MQ=255)                                                                                      cATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACTTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGt                                                               >  3:127739/1‑149 (MQ=255)                                                                                       aTGACCGACTCCGCATTCACTTCCCTGGTCGACATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTc                                                              >  7:298382/1‑149 (MQ=255)                                                                                        tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                             <  8:313245/149‑1 (MQ=255)                                                                                        tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                             <  7:365945/149‑1 (MQ=255)                                                                                           ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                              <  6:147911/81‑1 (MQ=255)                                                                                           ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                              >  5:147911/1‑81 (MQ=255)                                                                                            cGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGc                                                          >  6:377328/1‑148 (MQ=255)                                                                                                tCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCg                                                     >  1:318349/1‑149 (MQ=255)                                                                                                  cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                              >  8:204825/1‑106 (MQ=255)                                                                                                  cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                              <  7:204825/106‑1 (MQ=255)                                                                                                      ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCAc                                               <  3:124779/149‑1 (MQ=255)                                                                                                         aCTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGCCCTCCGCCATCTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGCCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCTGCAACCGCTGCACCGt                                            >  5:78585/1‑149 (MQ=255)                                                                                                          cTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTGAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                           <  1:279675/149‑1 (MQ=255)                                                                                                          cTTCCCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                           <  8:29160/149‑1 (MQ=255)                                                                                                           ttCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                          >  3:46803/1‑149 (MQ=255)                                                                                                             cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGGCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTGGGGCGGGGGCAAATGCGTGCTCAGGCGCCCCTTTTCCAGCAACCGCTGCACCGTGGCg                                        >  2:188409/1‑149 (MQ=255)                                                                                                             cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                         <  6:164764/148‑1 (MQ=255)                                                                                                             cccTGGGCGTCATGCTGACGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                         <  8:273413/148‑1 (MQ=255)                                                                                                                 ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGAGTCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                                       >  1:47358/1‑147 (MQ=255)                                                                                                                  gtcgtcATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGGCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAg                                    <  6:77799/148‑1 (MQ=255)                                                                                                                         tGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAAt                            >  5:8515/1‑149 (MQ=255)                                                                                                                           cTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATgg                          <  6:8515/149‑1 (MQ=255)                                                                                                                            tGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGc                         <  3:235763/149‑1 (MQ=255)                                                                                                                               cggcgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACGTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAgg                      >  6:167553/1‑149 (MQ=255)                                                                                                                                  cgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcgc                   >  2:329589/1‑149 (MQ=255)                                                                                                                                  ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                          >  6:299022/5‑126 (MQ=255)                                                                                                                                  ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                          <  5:299022/122‑1 (MQ=255)                                                                                                                                   gtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCt                  <  2:318349/149‑1 (MQ=255)                                                                                                                                   cgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGAGCTCAGGCGCCCCTTGTCCCGCAACAGCTGCACCGTGGCGACAGGCAGCATTGGCAGGCGCt                  >  4:238413/4‑149 (MQ=255)                                                                                                                                      ggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTgg                <  4:46803/148‑1 (MQ=255)                                                                                                                                       gcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc              >  7:219290/1‑149 (MQ=255)                                                                                                                                         cTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta            >  6:205116/1‑149 (MQ=255)                                                                                                                                          tGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtac           >  3:377633/1‑149 (MQ=255)                                                                                                                                            ggCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACTTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtacgt         <  4:127739/149‑1 (MQ=255)                                                                                                                                                   cAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGctc   >  3:80682/1‑149 (MQ=255)                                                                                                                                                   cAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGct    <  5:377328/148‑1 (MQ=255)                                                                                                                                                    aGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGCCCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGctct  >  6:82581/1‑149 (MQ=255)                                                                                                                                                     gcgcCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGctct  >  7:29572/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTCT  >  NC_002947/421922‑422216

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