Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,078,133	G→C	6.0%	A378G (GCC→GGC)	PP_2688 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,078,133	0	G	C	6.0%	212.9 / 3.7	67	A378G (GCC→GGC)	PP_2688	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

GAACGGGCGCAGGTCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCAGGATTGGTTCTTCGCTCAG  >  NC_002947/3077985‑3078278                                                                                                                                                     |                                                                                                                                                  gAACGGGCGCAGGTCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAgg                                                                                                                                                   >  6:195188/1‑149 (MQ=255)    cGGGCGCAGGTCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCgg                                                                                                                                                <  5:195188/149‑1 (MQ=255)       gcgcAGGTCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGcc                                                                                                                                             <  2:117186/149‑1 (MQ=255)              tCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCaa                                                                                                                                      >  6:195016/1‑149 (MQ=255)                 aTGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACAt                                                                                                                                   <  4:97276/149‑1 (MQ=255)                      aCGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCa                                                                                                                               >  8:60342/1‑148 (MQ=255)                      aCGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCAt                                                                                                                              <  5:208068/149‑1 (MQ=255)                       cgcgGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATa                                                                                                                             <  5:195016/149‑1 (MQ=255)                        gcgGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGGGCCCGCCTCGACCaaag                                                                                                                                    <  5:19148/141‑3 (MQ=255)                        gcgGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGGGCCCGCCTCGACCaaag                                                                                                                                    >  6:19148/1‑139 (MQ=255)                          ggTGCGATGCCGCTGTCGACGAAGCTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCGTCGATCTGGGCTTTGGTCGAGGCGGGCCCGCCCAACATGCCATAgcc                                                                                                                          >  2:54237/1‑149 (MQ=255)                                aTGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAATATGCCATAGCCGCCGGa                                                                                                                    <  3:222909/149‑1 (MQ=255)                                         tCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACAt                                                                                                                                   >  7:271027/1‑125 (MQ=255)                                         tCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACAt                                                                                                                                   <  8:271027/125‑1 (MQ=255)                                         tCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTggg                                                                                                            <  7:60342/148‑1 (MQ=255)                                           gacgaAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGtt                                                                                                          >  8:237276/1‑148 (MQ=255)                                              gaAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATACCATAGCCGCCGGAGCCCTGGGTTTCCt                                                                                                      >  4:5683/1‑149 (MQ=255)                                                 ggTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGa                                                                                                   >  1:211943/1‑149 (MQ=255)                                                  gTCGGGCAGGTCGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTg                                                                                                    <  3:167108/147‑1 (MQ=255)                                                  gTCGGGCAGGTAGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTg                                                                                                    >  4:167108/1‑147 (MQ=255)                                                     gggCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGaccac                                                                                               <  2:211943/149‑1 (MQ=255)                                                       gCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGAccacca                                                                                             <  8:262758/149‑1 (MQ=255)                                                            ttGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGTCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCGATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCg                                                                                        <  1:54237/149‑1 (MQ=255)                                                               aCAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCggg                                                                                      >  3:76780/1‑148 (MQ=255)                                                                gaGGGACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGTGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCtt                                                                                                     <  2:298226/131‑1 (MQ=255)                                                                cAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGCAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCAAAGCCGCCGGAGCCCTGGGTTTCCtt                                                                                                     >  1:298226/1‑132 (MQ=255)                                                                  ggCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAgg                                                                                  >  8:331644/1‑149 (MQ=255)                                                                   gCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACa                                                                                                                                    <  8:75044/98‑1 (MQ=255)                                                                   gCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACa                                                                                                                                    >  7:75044/1‑98 (MQ=255)                                                                   gCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCtt                                                                                                     <  8:128549/129‑1 (MQ=255)                                                                   gCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCtt                                                                                                     >  7:128549/1‑129 (MQ=255)                                                                   gCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGt                                                                                 >  4:335198/1‑149 (MQ=255)                                                                     acaGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCgg                                                                                       >  4:213742/1‑141 (MQ=255)                                                                     acaGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCgg                                                                                       <  3:213742/141‑1 (MQ=255)                                                                       aGCGTCGGTTGAGCGATATAGGCTTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGc                                                                             <  7:237276/149‑1 (MQ=255)                                                                           tCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGaccac                                                                                               <  4:223894/127‑1 (MQ=255)                                                                           tCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGaccac                                                                                               <  4:128126/127‑1 (MQ=255)                                                                           tCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGaccac                                                                                               >  3:223894/1‑127 (MQ=255)                                                                           tCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGaccac                                                                                               >  3:128126/1‑127 (MQ=255)                                                                               ttGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACAt                                                                     >  3:101730/1‑149 (MQ=255)                                                                               ttGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACAt                                                                     >  2:246947/1‑149 (MQ=255)                                                                                  aGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATgg                                                                   >  7:155326/1‑148 (MQ=255)                                                                                      atatAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCgg                                                                                                                                                <  2:197500/67‑1 (MQ=255)                                                                                      atatAGGCGTGCGGGCGGGCTTTGATGCCCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCgg                                                                                                                                                >  1:197500/1‑67 (MQ=255)                                                                                          aGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATcc                                                          >  6:236128/1‑149 (MQ=255)                                                                                              gTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATACCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATccgccg                                                      <  3:5683/149‑1 (MQ=255)                                                                                                   ggcgggcTTTGCTGCGCGCGCGGAAGTCCTCGATCCGGGCTTTGGTCGTGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTg                                                 <  8:155326/149‑1 (MQ=255)                                                                                                    gcgggcTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCa                                                                           >  7:244417/1‑122 (MQ=255)                                                                                                      gggcTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCgg                                              <  7:268279/149‑1 (MQ=255)                                                                                                          tttGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACt                                          >  1:166492/1‑149 (MQ=255)                                                                                                                  gcgcgcgGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTc                                  >  5:3725/1‑149 (MQ=255)                                                                                                                     cgcgGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGc                               <  5:236128/149‑1 (MQ=255)                                                                                                                       cgGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGt                                   >  7:156662/1‑143 (MQ=255)                                                                                                                       cgGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGt                                   <  8:156662/143‑1 (MQ=255)                                                                                                                             tCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACgttgtt                       >  4:305507/1‑149 (MQ=255)                                                                                                                              ccTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATccgccg                                                      <  7:961/117‑1 (MQ=255)                                                                                                                              ccTCGATCTGGGCTTTGGTCCAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTCACCACCAGTTCGGGCAGCTGGGCCAGTACATGGGACAGATccgcag                                                      >  8:961/1‑115 (MQ=255)                                                                                                                               cTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTc                                  <  8:249022/136‑1 (MQ=255)                                                                                                                               cTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTc                                  >  7:249022/1‑136 (MQ=255)                                                                                                                               cTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTc                      >  5:197657/1‑148 (MQ=255)                                                                                                                               cTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTc                      <  6:197657/148‑1 (MQ=255)                                                                                                                               cTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGCACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCa                     <  1:246947/149‑1 (MQ=255)                                                                                                                                  gATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCAgg                   <  7:331644/148‑1 (MQ=255)                                                                                                                                            ttGGTCGGGCCCGCCTCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTc                      >  4:332656/1‑135 (MQ=255)                                                                                                                                            ttGGTCGGGCCCGCCTCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTc                      <  3:332656/135‑1 (MQ=255)                                                                                                                                             tGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCAGGATTGGttcttc        <  4:101730/148‑1 (MQ=255)                                                                                                                                               gTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACAACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCATGCGGCACTGCCAGCTCGGCACGTTGTTCAGGATTGGTTCTTCGCt     >  2:114819/1‑149 (MQ=255)                                                                                                                                                  gAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCAGGATTGGTTCTTCGCTCAg  <  3:305507/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GAACGGGCGCAGGTCGATGTGACGCGGTGCGATGCCGCTGTCGACGAAGGTCGGGCAGGTTGACAGGCACAGCGTCGGTTGAGCGATATAGGCGTGCGGGCGGGCTTTGATGCGCGCGCGGAAGTCCTCGATCTGGGCTTTGGTCGAGGCGGGCCCGACCAACATGCCATAGCCGCCGGAGCCCTGGGTTTCCTTGACCACCAGTTCGGGCAGGTGGGCCAGTACATGGGACAGATCCGCCGGCTTGCGGCACTGCCAGGTCGGCACGTTGTTCAGGATTGGTTCTTCGCTCAG  >  NC_002947/3077985‑3078278

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