Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,710,174	T→C	6.6%	A102A (GCA→GCG)	PP_4168 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,710,174	0	T	C	6.6%	174.3 / 2.4	60	A102A (GCA→GCG)	PP_4168	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base C (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCA  >  NC_002947/4710032‑4710316                                                                                                                                               |                                                                                                                                                  gCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGCCGCGCATCCGGTGcgtga                                                                                                                                             <  3:109190/149‑5 (MQ=255)     tgtgCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGt                                                                                                                                         >  6:60681/1‑149 (MQ=255)       tgCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑gc                                                                                                                                    >  3:114532/1‑149 (MQ=255)          agcaCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCgct                                                                                                                                 <  7:296927/149‑1 (MQ=255)           gcaCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCgctg                                                                                                                                >  1:190792/1‑149 (MQ=255)               cTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCtt                                                                                                                             >  8:241721/1‑148 (MQ=255)               cTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCtt                                                                                                                             >  5:274060/1‑148 (MQ=255)               cTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTc                                                                                                                            <  5:300009/149‑1 (MQ=255)                tGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCa                                                                                                                           >  6:191321/1‑149 (MQ=255)                 gCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACGGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCAc                                                                                                                          >  5:200041/1‑149 (MQ=255)                   tCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCt                                                                                                                        >  4:279301/1‑149 (MQ=255)                              ggCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCAc                                                                                                              <  2:254812/148‑1 (MQ=255)                               gCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCAcgc                                                                                                            >  2:108751/1‑149 (MQ=255)                                cAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCAcgcg                                                                                                           <  6:7349/149‑1 (MQ=255)                                     aaCCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGt                                                                                                                                                   <  8:133644/107‑1 (MQ=255)                                        cGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACGGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGggccggc                                                                                                   <  6:200041/149‑1 (MQ=255)                                            gcTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGt                                                                                               <  5:106830/149‑1 (MQ=255)                                                ggCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGc                                                                                           >  3:103773/1‑149 (MQ=255)                                                    tCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAAc                                                                                       <  4:123434/149‑1 (MQ=255)                                                             gCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACgg                                                                               <  1:108751/148‑1 (MQ=255)                                                             gCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGCGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGc                                                                              >  8:291191/1‑149 (MQ=255)                                                              cccGCGGCTCAGCTTCACCGTGGCGCAGCGCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCt                                                                             >  8:230543/1‑149 (MQ=255)                                                               ccGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCtt                                                                            >  2:245055/1‑149 (MQ=255)                                                                   ggCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTgg                                                                        <  7:241721/149‑1 (MQ=255)                                                                   ggCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTgg                                                                        <  6:274060/149‑1 (MQ=255)                                                                      tCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCg                                                                     >  4:103151/1‑149 (MQ=255)                                                                       cAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGc                                                                    >  2:217877/1‑149 (MQ=255)                                                                        aGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCt                                                                   <  4:114532/149‑1 (MQ=255)                                                                        aGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCt                                                                   <  4:103773/149‑1 (MQ=255)                                                                           ttCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCtgtg                                                                >  3:16576/1‑149 (MQ=255)                                                                                cGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAAc                                                           <  3:103151/149‑1 (MQ=255)                                                                                      gcagcaCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCg                                                      >  5:304684/1‑148 (MQ=255)                                                                                        agcaCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGAGGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGcc                                                       >  1:263719/1‑148 (MQ=255)                                                                                          caCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCg                                                 >  3:275734/1‑149 (MQ=255)                                                                                             ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCt                                                                                             >  8:43465/1‑102 (MQ=255)                                                                                             ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCt                                                                                             <  7:43465/102‑1 (MQ=255)                                                                                             ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCtgt                                                                 >  6:322963/1‑130 (MQ=255)                                                                                             ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGc                                               >  5:209308/1‑148 (MQ=255)                                                                                             ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGc                                               <  6:209308/148‑1 (MQ=255)                                                                                               acagcttcacagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGa                                            >  3:227188/1‑149 (MQ=255)                                                                                                       acagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGa                                            >  8:262247/1‑141 (MQ=255)                                                                                                       acagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGa                                            <  7:262247/141‑1 (MQ=255)                                                                                                       acagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGc                                     >  1:273117/1‑148 (MQ=255)                                                                                                       acagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGCGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGAt                                           >  4:302925/1‑142 (MQ=255)                                                                                                       acagcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGCGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGAt                                           <  3:302925/142‑1 (MQ=255)                                                                                                         agcttGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAg                                  <  7:230543/149‑1 (MQ=255)                                                                                                              gggTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTAc                              <  1:245055/148‑1 (MQ=255)                                                                                                                gTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCt                                         <  3:316724/135‑1 (MQ=255)                                                                                                                gTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCt                                         >  4:316724/1‑135 (MQ=255)                                                                                                                gTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTgg                           <  4:275734/149‑1 (MQ=255)                                                                                                                 ttCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGc                          >  4:28174/1‑149 (MQ=255)                                                                                                                   ccTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                         >  3:232860/1‑148 (MQ=255)                                                                                                                    cTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCa                       <  3:279301/149‑1 (MQ=255)                                                                                                                      cGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                     <  5:60681/149‑1 (MQ=255)                                                                                                                       gTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                     <  3:19104/148‑1 (MQ=255)                                                                                                                             cACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGc              <  4:16576/149‑1 (MQ=255)                                                                                                                                 ggATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGaa          <  2:273117/149‑1 (MQ=255)                                                                                                                                 ggATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGaa          <  3:28174/149‑1 (MQ=255)                                                                                                                                 ggATGCGCGGCCGGCGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGaa          <  7:291191/149‑1 (MQ=255)                                                                                                                                      cgcgGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtc      <  2:35973/148‑1 (MQ=255)                                                                                                                                         ggCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtcatca  <  1:217877/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  GCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGCATCCGTCTTGGCCCGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGT‑‑‑GCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCA  >  NC_002947/4710032‑4710316

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