Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,067	C→G	6.2%	L627L (CTG→CTC) ‡	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,067	0	C	G	6.2%	249.9 / 6.6	81	L627L (CTG→CTC) ‡	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/40);  new base G (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTC  >  NC_002947/421920‑422215                                                                                                                                                     |                                                                                                                                                     gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑Tc                                                                                                                                                      <  4:195187/148‑1 (MQ=255) gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCa                                                                                                                                                     <  7:73990/149‑1 (MQ=255)   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                                                                                                                              >  5:77799/1‑148 (MQ=255)                         cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                              >  6:170506/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                                                                                                             >  3:359345/1‑132 (MQ=255)                                                          gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                             <  4:359345/132‑1 (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                                                                                   >  4:329098/1‑149 (MQ=255)                                                                   cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                                   >  7:313245/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                                                            <  7:365945/149‑1 (MQ=255)                                                                                          tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                            <  8:313245/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                                                                                             <  7:204825/106‑1 (MQ=255)                                                                                                    cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                             >  8:204825/1‑106 (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                                         <  5:299022/122‑1 (MQ=255)                                                                                                                                    ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                         >  6:299022/5‑126 (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)                                                                                                                                                     |                                                                                                                                                     GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTC  >  NC_002947/421920‑422215

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