Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	835,966	G→C	6.8%	T38S (ACC→AGC)	PP_0718 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	835,966	0	G	C	6.8%	183.0 / 4.0	59	T38S (ACC→AGC)	PP_0718	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/30)
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.90e-01

CAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTC  >  NC_002947/835818‑836113                                                                                                                                                      |                                                                                                                                                    cAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCAtggtgg                                                                                                                                                     >  4:350994/1‑149 (MQ=255) cAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCAtggtgg                                                                                                                                                     <  6:60316/149‑1 (MQ=255) cAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCAtggtgg                                                                                                                                                     <  7:74141/149‑1 (MQ=255)           tCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           >  1:143376/1‑149 (MQ=255)           tCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           <  3:350994/149‑1 (MQ=255)             cGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGc                                                                                                                                         >  5:231164/1‑149 (MQ=255)                  ccAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCcggcg                                                                                                                                    >  7:47688/1‑149 (MQ=255)                     gCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 <  1:180724/149‑1 (MQ=255)                       aCCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGAt                                                                                                                               <  2:143376/149‑1 (MQ=255)                       aCCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGAt                                                                                                                               <  3:123284/149‑1 (MQ=255)                                 ggCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTggg                                                                                                                      >  4:218922/1‑148 (MQ=255)                                   cTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCtt                                                                                                                         <  3:90790/143‑1 (MQ=255)                                   cTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCtt                                                                                                                         >  4:90790/1‑143 (MQ=255)                                      ggACCGAACAGGGCATAGGCCAGCAGGGGGATGGTGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctcgt                                                                                                                                     >  4:295056/1‑124 (MQ=255)                                      ggACCGAACAGGGCATAGGCCAGCAGGGGGATGGTGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctcgt                                                                                                                                     <  3:295056/128‑5 (MQ=255)                                       gACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAg                                                                                                               >  8:197997/1‑149 (MQ=255)                                              cAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCt                                                                                                        <  3:112497/149‑1 (MQ=255)                                               aGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAAGCTGCCtt                                                                                                       >  6:317845/1‑149 (MQ=255)                                                 ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 <  4:272445/121‑1 (MQ=255)                                                 ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 >  3:272445/1‑121 (MQ=255)                                                 ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCCGCGATGTCCTTGGGTAGCCAGGCTGCCTTg                                                                                                      >  7:277418/1‑148 (MQ=255)                                                   cATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAg                                                                                                   <  1:21935/149‑1 (MQ=255)                                                      aGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctc                                                                                                                                       <  2:128444/110‑3 (MQ=255)                                                      aGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctc                                                                                                                                       >  1:128444/1‑108 (MQ=255)                                                            gcagGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  6:237875/1‑149 (MQ=255)                                                                  ggATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGAc                                                                                     >  5:195124/1‑148 (MQ=255)                                                                           tAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCggg                                                                            >  3:192480/1‑148 (MQ=255)                                                                            aGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGcc                                                                                                                 <  5:16250/110‑1 (MQ=255)                                                                            aGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGcc                                                                                                                 >  6:16250/1‑110 (MQ=255)                                                                              cGTACAGGCCATAGATACCGGGGACGCCTGAAGC‑CTCGCCGTAGTCGATGGCCACCGGCACGAGCATGTTGGCGAGCACCAGGCCGGCGTCGATGGCCTTGGGTAGCCTGGCTGCCATGTAGTGGAGCAGGGTGACCAGGCCGGGTAAc                                                                        <  4:302156/149‑1 (MQ=255)                                                                              cGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           <  3:147693/82‑1 (MQ=255)                                                                              cGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           >  4:147693/1‑82 (MQ=255)                                                                                  cAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCa                                                                                            >  4:260990/1‑125 (MQ=255)                                                                                  cAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCa                                                                                            <  3:260990/125‑1 (MQ=255)                                                                                     gCCATAGATACCGGGGACGCCGGAAGCGCTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGt                                                                  >  1:16088/1‑149 (MQ=255)                                                                                      ccATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTg                                                                <  2:16088/149‑1 (MQ=255)                                                                                       cATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGc                                                               <  8:47688/149‑1 (MQ=255)                                                                                           gATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGTTAGCCAAGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGt                                                           <  5:317845/149‑1 (MQ=255)                                                                                            aTACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTc                                                          <  7:197997/149‑1 (MQ=255)                                                                                              aCCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          <  5:109291/115‑1 (MQ=255)                                                                                              aCCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  6:109291/1‑115 (MQ=255)                                                                                                   ggACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGa                                                   <  3:191122/149‑1 (MQ=255)                                                                                                               c‑cTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTgggg                                        >  6:146774/1‑148 (MQ=255)                                                                                                                    ggCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTg                                                                                                 >  2:222059/1‑87 (MQ=255)                                                                                                                    ggCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTg                                                                                                 <  1:222059/87‑1 (MQ=255)                                                                                                                       gTAGGCGACGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGgtgt                                >  7:167245/1‑149 (MQ=255)                                                                                                                          ggCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCgg                             >  3:23376/1‑149 (MQ=255)                                                                                                                          ggCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCgg                             >  2:165958/1‑149 (MQ=255)                                                                                                                           gCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGc                            >  3:330057/1‑149 (MQ=255)                                                                                                                           gCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGc                            >  7:74995/1‑149 (MQ=255)                                                                                                                            cGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGct                           <  5:146774/149‑1 (MQ=255)                                                                                                                              acgCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGctct                         <  8:167245/147‑1 (MQ=255)                                                                                                                                gCCCACAGGCACGAGCATGGTGGTGAGCACGAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTggg                                              <  3:248809/126‑1 (MQ=255)                                                                                                                                gCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTggg                                              >  4:248809/1‑126 (MQ=255)                                                                                                                                 cccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTgg                       <  6:195124/148‑1 (MQ=255)                                                                                                                                  ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  2:257706/1‑80 (MQ=255)                                                                                                                                  ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          <  1:257706/80‑1 (MQ=255)                                                                                                                                  ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGc                     >  3:146329/1‑149 (MQ=255)                                                                                                                                  ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGc                     <  4:330057/149‑1 (MQ=255)                                                                                                                                     caGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCgag                  <  8:74995/149‑1 (MQ=255)                                                                                                                                        gCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCTTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAg               <  5:237875/149‑1 (MQ=255)                                                                                                                                           cGAGCATGGTGTTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGgc            <  6:231164/149‑1 (MQ=255)                                                                                                                                                     ggtgAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGACGGGCGCTGAACCTc  <  4:192480/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTC  >  NC_002947/835818‑836113

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