Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,776,111	C→G	7.2%	A153G (GCG→GGG)	sotB →	sugar efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,776,111	0	C	G	7.2%	160.7 / 3.6	56	A153G (GCG→GGG)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (22/29);  new base G (2/2);  total (24/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 = 8.33e-01

ATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2775963‑2776256                                                                                                                                                     |                                                                                                                                                   aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCATCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgc                                                                                                                                                    >  3:178561/1‑149 (MQ=255)  ttGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                   <  1:218456/149‑1 (MQ=255)  ttGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                   >  2:218456/1‑149 (MQ=255)  ttGCCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTCGTCAGCCCCATCGGCATTGCCCTGGCGCATGCGGAGTTATGGTCCATCACCGCCTCGCTGCCGGAGCGTGAGGCGCCGCCGGGCCAGCAGGACAACgcc                                                                                                                                                   >  4:16085/1‑148 (MQ=255)       aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                               <  1:294009/148‑1 (MQ=255)       aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                               >  2:294009/1‑148 (MQ=255)            cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                         <  4:268506/149‑1 (MQ=255)             ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                        >  5:196983/1‑149 (MQ=255)               ccTGTCCTGGCTTTCGCAAAGCGTCTCCCTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCAGGCTGTGTTCTGGTCCATCACCGCCGCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                        <  6:257036/146‑1 (MQ=255)                cTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgc                                                                                                                                         >  7:60989/1‑145 (MQ=255)                cTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgc                                                                                                                                         <  8:60989/145‑1 (MQ=255)                 tGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGc                                                                                                                                     >  2:300132/1‑148 (MQ=255)                   tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAc                                                                                                                                  >  1:125704/1‑149 (MQ=255)                      tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                               >  7:89470/1‑149 (MQ=255)                        gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                               <  8:73197/147‑1 (MQ=255)                            tCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtact                                                                                                                           <  3:41874/147‑1 (MQ=255)                             cGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtactac                                                                                                                         <  3:198998/148‑1 (MQ=255)                                aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCGGCCAAGGCGCAGGGCCTGCAGGCCACCGGTACTACCC‑Tg                                                                                                                    >  2:135700/1‑149 (MQ=255)                                    cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCT‑GCCACCGGTACTACCC‑TGGCCAt                                                                                                               >  2:83665/1‑149 (MQ=255)                                     ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCAGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCAt                                                                                                               >  6:28484/1‑149 (MQ=255)                                        gCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATgg                                                                                                             <  5:141405/148‑1 (MQ=255)                                         cgaTGCGGCTGACCAGCAGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGTAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTc                                                                                                           >  8:263758/3‑149 (MQ=255)                                          cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATgg                                                                                                             <  4:56073/146‑1 (MQ=255)                                           aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCt                                                                                                         >  1:214244/1‑149 (MQ=255)                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTggg                                                                                                      <  4:307143/149‑1 (MQ=255)                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTgg                                                                                                       <  5:194518/148‑1 (MQ=255)                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTgg                                                                                                       <  3:202261/148‑1 (MQ=255)                                              ctgctgGTCAGCCGCATCGACACTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGTTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTgg                                                                                                       <  6:75188/148‑1 (MQ=255)                                                      cAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGGGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGc                                                                                              >  1:312718/1‑149 (MQ=255)                                                       aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGGCAAGGCGCTGGGCCAGCTGGCCCCCGGTACAACCC‑TGGCCCATGTCCTGGGCCACCCGcg                                                                                             >  7:153247/1‑148 (MQ=255)                                                            cATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGcc                                                                                        >  7:199462/1‑149 (MQ=255)                                                            cATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGcc                                                                                        <  2:125704/149‑1 (MQ=255)                                                                ggCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCgtgt                                                                                    >  6:109413/1‑149 (MQ=255)                                                                ggCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCgtgt                                                                                    <  5:101380/149‑1 (MQ=255)                                                                   aTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTgtggt                                                                                 <  5:139362/149‑1 (MQ=255)                                                                   aTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAACC‑AGGCCATGGTCCAGGGCCACCCACGGGGCCGggtggt                                                                                 >  4:24985/1‑149 (MQ=255)                                                                    ttGCCGTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTgtggtg                                                                                <  6:196983/149‑1 (MQ=255)                                                                       cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCGGCCACAGCGCTGGGCCAGCGGGCCACCGGTACTACCCATGGCCA‑CGTCCGGGGCATCCCGCTGGGCCGGGGGGTGGGc                                                                             >  5:206282/1‑149 (MQ=255)                                                                                     gtgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTg                                                               <  2:214244/149‑1 (MQ=255)                                                                                       gtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  5:109413/149‑1 (MQ=255)                                                                                        ttCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  3:314632/148‑1 (MQ=255)                                                                                        ttCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  8:89470/148‑1 (MQ=255)                                                                                        ggcggggccATCCCCGCCTCGTTGGCTGTGCTGGTGGCGCCGCCGGGCCAGCAGTCCAAGGAGCTGGGCCTGCTGGCCACCGGTACTCCCC‑TGGCCATGTGCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  6:206282/141‑1 (MQ=255)                                                                                                 aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGt                                                    >  2:3413/1‑148 (MQ=255)                                                                                                    aCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTcc                                                 <  1:300132/148‑1 (MQ=255)                                                                                                      cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  6:254515/149‑1 (MQ=255)                                                                                                      cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  2:312718/149‑1 (MQ=255)                                                                                                        ccTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGc                                            <  3:24985/149‑1 (MQ=255)                                                                                                                gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                    >  2:32029/1‑149 (MQ=255)                                                                                                                    tGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCCTGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgggg                                  >  3:58914/1‑147 (MQ=255)                                                                                                                          tgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                          <  8:199462/149‑1 (MQ=255)                                                                                                                          tgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                          <  4:178561/149‑1 (MQ=255)                                                                                                                                 gccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca                   <  1:135700/149‑1 (MQ=255)                                                                                                                                  ccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACAC‑TGGCCATGGTCCTGGGCATACCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCAt                  >  2:78856/1‑149 (MQ=255)                                                                                                                                     ggCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCt               >  8:219271/1‑149 (MQ=255)                                                                                                                                        cagcagGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                   >  4:289393/1‑126 (MQ=255)                                                                                                                                        cagcagGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                   <  3:289393/126‑1 (MQ=255)                                                                                                                                         agcagGCCATGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTTGTGGGCGAGGCACTGGGCTTGCGTATCACTTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           <  5:276943/149‑1 (MQ=255)                                                                                                                                          gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          <  5:28484/149‑1 (MQ=255)                                                                                                                                                  aGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCaaa  >  5:38450/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCC‑TGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2775963‑2776256

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