Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	665,421	G→C	6.6%	G763R (GGC→CGC)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	665,421	0	G	C	6.6%	197.4 / 4.2	61	G763R (GGC→CGC)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (2/2);  total (29/32)
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.94e-01

TCCAGTGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTGGTTAT  >  NC_002947/665274‑665563                                                                                                                                                     |                                                                                                                                               tCCAGTGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCgg                                                                                                                                                <  1:194399/148‑1 (MQ=255)   cAGTGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCggg                                                                                                                                               >  7:38657/1‑147 (MQ=255)    aGTGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAg                                                                                                                                            >  4:13907/1‑149 (MQ=255)      tGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCg                                                                                                                                          >  3:72359/1‑149 (MQ=255)          cgcTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCCCGCAGCgggcg                                                                                                                                      <  4:11431/149‑4 (MQ=255)          cgcTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCCCGCAGCgggcg                                                                                                                                      <  7:312784/149‑4 (MQ=255)                  ggCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGc                                                                                                                                           >  4:133771/1‑136 (MQ=255)                  ggCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGc                                                                                                                                           <  3:133771/136‑1 (MQ=255)                  ggCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCtgat                                                                                                                                <  2:192510/147‑1 (MQ=255)                  ggCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCtgat                                                                                                                                <  4:72359/147‑1 (MQ=255)                    ccGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATGTGCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTgatgat                                                                                                                             <  1:97483/149‑1 (MQ=255)                    ccGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAg                                                                                                                                            >  2:243031/1‑133 (MQ=255)                    ccGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAg                                                                                                                                            <  1:243031/133‑1 (MQ=255)                    ccGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATc                                                                                                                            <  5:266227/149‑1 (MQ=255)                    ccGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATc                                                                                                                            <  6:39849/149‑1 (MQ=255)                      ggCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGa                                                                                                                          <  5:136618/149‑1 (MQ=255)                       gCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAg                                                                                                                         <  4:36825/149‑1 (MQ=255)                       gCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAg                                                                                                                         >  4:173726/1‑149 (MQ=255)                       gCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAg                                                                                                                         >  4:169185/1‑149 (MQ=255)                           ccGTGCCCACGACAACGCTGAGCGCGGCGCGGTGCGGCCGGACCACGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGAAG‑GCGCTACCGCCAAGGACACCGTCACTACGCCCGCTGCGGGCAGCGCCTACCTGATGATCGAGTcaa                                                                                                                     >  4:14616/1‑147 (MQ=255)                              tGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGc                                                                                                                   >  5:284280/1‑148 (MQ=255)                                 ccAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATg                                                                                                               >  3:327378/1‑149 (MQ=255)                                   aGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGc                                                                                                             >  2:50127/1‑149 (MQ=255)                                             tGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTgtggt                                                                                                   >  4:122440/1‑149 (MQ=255)                                               aGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTgtggt                                                                                                   <  1:50127/147‑1 (MQ=255)                                                    ggcgCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGga                                                                                            <  8:38657/149‑1 (MQ=255)                                                         ggTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCg                                                                                       <  3:173726/149‑1 (MQ=255)                                                         ggTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCg                                                                                       <  3:169185/149‑1 (MQ=255)                                                          gTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGa                                                                                      >  3:189504/1‑149 (MQ=255)                                                           tGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAt                                                                                     <  3:13907/149‑1 (MQ=255)                                                           tGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAt                                                                                     <  5:81211/149‑1 (MQ=255)                                                           tGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAt                                                                                     <  8:102526/149‑1 (MQ=255)                                                           tGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAt                                                                                     <  4:189504/149‑1 (MQ=255)                                                             cGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAtgt                                                                                   <  7:241274/149‑1 (MQ=255)                                                             cGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAtgt                                                                                   <  7:236814/149‑1 (MQ=255)                                                              ggCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGAtgtg                                                                                  <  2:120029/149‑1 (MQ=255)                                                                ccGGACCAGGTCAAGCTGGCGCTGGACAACCCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGAGCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGAGCGATGTGcc                                                                                >  3:98170/1‑149 (MQ=255)                                                                    aCCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCccgcc                                                                            >  4:30596/1‑149 (MQ=255)                                                                       aGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAg                                                                         <  4:327378/149‑1 (MQ=255)                                                                        ggTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAgg                                                                                             <  2:204394/128‑1 (MQ=255)                                                                           cAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGc                                                                      >  7:182822/1‑148 (MQ=255)                                                                             aGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGc                                                                      >  2:167387/1‑146 (MQ=255)                                                                                        gACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGt                                                        >  1:310261/1‑149 (MQ=255)                                                                                             aCCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCt                                                                                                        <  3:301028/96‑1 (MQ=255)                                                                                             aCCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCt                                                                                                        >  4:301028/1‑96 (MQ=255)                                                                                               cctcctACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTg                                                 <  3:122440/149‑1 (MQ=255)                                                                                                     acacCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGt                                                        <  1:265639/136‑1 (MQ=255)                                                                                                        ccGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGcc                                        <  3:30596/149‑1 (MQ=255)                                                                                                         cGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTgtggtg                                                                                                  <  4:57183/90‑1 (MQ=255)                                                                                                           aTG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTgg                                     <  4:164097/149‑1 (MQ=255)                                                                                                                cgcTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGTCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGTCAGGAGATCGATGTGCCCGCCGAGGGCAAGACGTTTGACGTGCAACTGGACAAGGCCTGGGccc                                 <  3:14616/149‑1 (MQ=255)                                                                                                                       gCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGa                          >  2:89405/1‑149 (MQ=255)                                                                                                                        ccAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGAc                         >  4:56509/1‑149 (MQ=255)                                                                                                                          aaGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCt                       >  5:182312/1‑149 (MQ=255)                                                                                                                                        tCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGt         >  6:196292/1‑149 (MQ=255)                                                                                                                                        tCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGt         >  5:120802/1‑149 (MQ=255)                                                                                                                                        tCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGt         >  1:134382/1‑149 (MQ=255)                                                                                                                                         gtagccgcTGCCCGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGGTGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGtt        >  4:299449/5‑149 (MQ=255)                                                                                                                                         ccgcccgcTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGtt        >  6:138083/1‑149 (MQ=255)                                                                                                                                          tagccgcTGCCCGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTg       >  5:324042/4‑149 (MQ=255)                                                                                                                                          cgcccgcTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTg       <  6:182312/149‑1 (MQ=255)                                                                                                                                             ccgcTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTGGt     >  5:40354/1‑148 (MQ=255)                                                                                                                                               gcTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTGGTTAt  >  5:254841/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               TCCAGTGAGCGCTTCTGGGCCGGCTACCGTGCCCAGGACAACGCTGAGGGCGGCGCGGTGCGGCCGGACCAGGTCAAGCTGGCGCTGGACAAACCCTCCTACACCGATG‑GCGCTACCGCCAAGGTCACCGTCACTCCGCCCGCTGCGGGCAGCGGCTACCTGATGATCGAGTCCAGCGATGGCCCGCTGTGGTGGCAGGAGATCGATGTGCCCGCCGAGGGCAAGACCTTTGACGTGCAACTGGACAAGGCCTGGGCCCGTCACGACCTGTACATCAGCGCGTTGGTTAT  >  NC_002947/665274‑665563

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