Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	664,887	A→G	6.9%	S585G (AGC→GGC)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	664,887	0	A	G	6.9%	164.8 / 2.4	58	S585G (AGC→GGC)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/2);  total (29/29)
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.59e-01

CCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCG  >  NC_002947/664742‑665034                                                                                                                                                   |                                                                                                                                                    ccAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGcc                                                                                                                                                  >  1:149530/1‑149 (MQ=255)    aTTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTg                                                                                                                                                >  3:34768/1‑148 (MQ=255)         gCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAg                                                                                                                                          >  3:317488/1‑149 (MQ=255)          cTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTTGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCTGGCCTGCacggtc                                                                                                                                         <  2:345154/149‑6 (MQ=255)                 aCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAg                                                                                                                                             <  3:142647/138‑1 (MQ=255)                 aCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAg                                                                                                                                             <  3:137339/138‑1 (MQ=255)                 aCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAg                                                                                                                                             >  4:137339/1‑138 (MQ=255)                 aCCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAg                                                                                                                                             >  4:142647/1‑138 (MQ=255)                  ccGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCggg                                                                                                                                     >  2:304563/1‑145 (MQ=255)                     aaaCCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccgg                                                                                                                              <  7:44175/149‑1 (MQ=255)                           aTCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGAt                                                                                                                         >  7:231063/1‑148 (MQ=255)                               gAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCa                                                                                                                                                     >  5:15827/1‑116 (MQ=255)                               gAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCa                                                                                                                                                     <  6:15827/116‑1 (MQ=255)                               gAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAccc                                                                                                                    >  3:20290/1‑149 (MQ=255)                                caCCACGACCTGGAACTGGACGAAGTCACCCTTCACCAGGCCGGCAAGGCGGTGGTCGATATCGAAACCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTCACCCTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCg                                                                                                                   <  2:136624/148‑1 (MQ=255)                                aaCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAcc                                                                                                                     <  3:193089/147‑1 (MQ=255)                                aaCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCg                                                                                                                   <  2:131231/149‑1 (MQ=255)                                aaCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCg                                                                                                                   >  8:289104/1‑149 (MQ=255)                                 aCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCTGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGACTCGGGTGGCCGGCCGATCACCCgg                                                                                                                  >  8:246699/1‑149 (MQ=255)                                  ccAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGATCGGCCGGCCACCCgac                                                                                                                 <  1:71201/149‑3 (MQ=255)                                   ccggcccTGGAACTGGACGACGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCTATATCGAACGCCGT‑TGGGCCGAATCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcg                                                                                                                <  4:34768/144‑1 (MQ=255)                                     ggACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccggc                                                                                                                             <  6:153673/134‑1 (MQ=255)                                     ggACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccggc                                                                                                                             >  5:153673/1‑134 (MQ=255)                                        ccTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTg                                                                                                            >  5:165100/1‑148 (MQ=255)                                         cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTg                                                                                                            >  6:74254/1‑147 (MQ=255)                                          tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                         >  2:134623/1‑149 (MQ=255)                                              aCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                     <  4:317488/149‑1 (MQ=255)                                              aCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                     <  1:269173/149‑1 (MQ=255)                                                             ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTGTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                >  1:69831/1‑140 (MQ=255)                                                             ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCt                                                                                                   >  1:208439/1‑136 (MQ=255)                                                             ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCt                                                                                                   <  2:208439/136‑1 (MQ=255)                                                             ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                <  2:69831/139‑1 (MQ=255)                                                             ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCt                                                                                       >  3:202198/1‑148 (MQ=255)                                                              cTTCACCAGGCCGGCAAGCCGGTGCTCGATATCCAAAGCCGT‑CGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGa                                                                                     <  4:29519/149‑1 (MQ=255)                                                                  aCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAgcgc                                                                                 >  7:32258/1‑149 (MQ=255)                                                                  aCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAgcgc                                                                                 <  5:74254/149‑1 (MQ=255)                                                                         cGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                          <  4:20290/149‑1 (MQ=255)                                                                          ggCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGc                                                                             <  2:9451/145‑1 (MQ=255)                                                                          ggCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGc                                                                             >  1:9451/1‑145 (MQ=255)                                                                             aaGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCt                                                                      <  2:149530/149‑1 (MQ=255)                                                                                   gtggtCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggcc                                                                                                                                >  3:202247/1‑85 (MQ=255)                                                                                   gtggtCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggcc                                                                                                                                <  4:202247/85‑1 (MQ=255)                                                                                       tCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTg                                                            <  7:289104/149‑1 (MQ=255)                                                                                         gATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTgat                                                          >  4:222397/1‑147 (MQ=255)                                                                                                             ccGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGAt                                       <  8:231063/149‑1 (MQ=255)                                                                                                              cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCTGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                      <  2:305576/149‑1 (MQ=255)                                                                                                                  gCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGa                                  >  7:62304/1‑149 (MQ=255)                                                                                                                   cccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGAc                                 <  4:202198/149‑1 (MQ=255)                                                                                                                      gCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACggg                              >  2:193960/1‑149 (MQ=255)                                                                                                                                    cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCGGGCGCCTGGAGCAGCCTATGTGGCCGGCTGAGAGCCTGGCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                <  7:246699/149‑1 (MQ=255)                                                                                                                                    cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTc                  >  2:64697/1‑147 (MQ=255)                                                                                                                                       aCCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGtt             >  6:271815/1‑149 (MQ=255)                                                                                                                                        tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgaggag                                             <  2:97773/113‑1 (MQ=255)                                                                                                                                        tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgaggag                                             >  1:97773/4‑116 (MQ=255)                                                                                                                                        tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTc            >  2:255397/4‑149 (MQ=255)                                                                                                                                         cGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCg                 >  2:275517/1‑143 (MQ=255)                                                                                                                                         cGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCg                 <  1:275517/143‑1 (MQ=255)                                                                                                                                             cAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggt       <  1:64697/149‑1 (MQ=255)                                                                                                                                               ggCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtgg     <  3:222397/149‑1 (MQ=255)                                                                                                                                                  cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGTCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCg  >  4:30358/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGT‑TGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCG  >  NC_002947/664742‑665034

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