Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,065,807	C→G	8.2%	A51A (GCG→GCC)	PP_3582 ←	RND transporter outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,065,807	0	C	G	8.2%	194.0 / 6.7	61	A51A (GCG→GCC)	PP_3582	RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base G (2/3);  total (30/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 = 3.87e-01

ATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAA  >  NC_002947/4065662‑4065954                                                                                                                                                   |                                                                                                                                                    aTTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCgccgcc                                                                                                                                                  >  4:107120/1‑149 (MQ=255) aTTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCgccgcc                                                                                                                                                  >  4:84115/1‑149 (MQ=255)     cgcCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGc                                                                                                                                              >  7:261092/1‑149 (MQ=255)           ggTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCa                                                                                                                                        <  6:161446/149‑1 (MQ=255)                 gTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagca                                                                                                                                    <  1:210740/147‑5 (MQ=255)                     ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                              <  3:107120/149‑1 (MQ=255)                     ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                              >  7:70516/1‑149 (MQ=255)                       gcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggt                                                                                                                            <  5:3558/149‑1 (MQ=255)                        cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                           <  3:202666/149‑1 (MQ=255)                        cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                           <  2:260550/149‑1 (MQ=255)                        cGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggt                                                                                                                            <  2:138185/148‑1 (MQ=255)                                ctctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTg                                                                                                                                               >  8:10678/1‑121 (MQ=255)                                ctctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTg                                                                                                                                               <  7:10678/121‑1 (MQ=255)                                 tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                           >  2:151163/1‑140 (MQ=255)                                 tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                           <  1:151163/140‑1 (MQ=255)                                 tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCggg                                                                                                                  >  6:118303/1‑149 (MQ=255)                                   ttCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGc                                                                                                                >  3:35498/1‑149 (MQ=255)                                     cGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGc                                                                                                                <  2:242200/147‑1 (MQ=255)                                      gACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCa                                                                                                             <  8:70516/149‑1 (MQ=255)                                        ccAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAgg                                                                                                           >  1:172998/1‑149 (MQ=255)                                        ccAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAgg                                                                                                           <  2:33365/149‑1 (MQ=255)                                         cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                          <  3:84115/149‑1 (MQ=255)                                         cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                          >  4:265640/1‑149 (MQ=255)                                         cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                          >  1:253068/1‑149 (MQ=255)                                          aGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGt                                                                                                         <  3:265640/149‑1 (MQ=255)                                             gcATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACa                                                                                                      <  8:308197/149‑1 (MQ=255)                                                          cATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAg                                                                                         <  4:35498/149‑1 (MQ=255)                                                           aTCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAgg                                                                                        >  1:8669/1‑149 (MQ=255)                                                           aTCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAgg                                                                                        >  4:239986/1‑149 (MQ=255)                                                               ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                    >  6:70770/1‑149 (MQ=255)                                                               ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                    >  5:34316/1‑149 (MQ=255)                                                               ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                                     >  2:72777/1‑97 (MQ=38)                                                               ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                                     <  1:72777/100‑4 (MQ=38)                                                                cGTAGATTTCCCACCAGGCGCCGCGGGCGAGTGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                                    >  1:202916/1‑149 (MQ=255)                                                                cGTAGATTTACCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACg                                                                                   >  5:291212/1‑149 (MQ=255)                                                                     aTTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACaa                                                                              >  3:313645/1‑149 (MQ=255)                                                                           caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaa                                                                                                                                       <  1:223162/86‑2 (MQ=37)                                                                           caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaa                                                                                                                                       >  2:223162/1‑85 (MQ=37)                                                                           caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCT‑GGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCa                                                                       >  6:78919/1‑149 (MQ=255)                                                                                ggCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTTCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGa                                                                   <  2:172998/149‑1 (MQ=255)                                                                                     cgcgGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAgg                                                              >  7:243689/1‑149 (MQ=255)                                                                                          gCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCTGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGa                                                         <  5:193199/149‑1 (MQ=255)                                                                                            gA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgc                                                       <  2:253068/149‑1 (MQ=255)                                                                                                 gCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacg                                                   >  4:296604/1‑149 (MQ=255)                                                                                                            cGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCa                                        <  3:239986/149‑1 (MQ=255)                                                                                                            cGAGGCCTGGATCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGt                                          >  1:34058/1‑147 (MQ=255)                                                                                                             gAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAg                                       >  8:153327/1‑149 (MQ=255)                                                                                                              aGGCCT‑GGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTg                                     <  5:78919/149‑1 (MQ=255)                                                                                                                gCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGc                                    >  8:144551/1‑149 (MQ=255)                                                                                                                gCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGc                                    >  7:62655/1‑149 (MQ=255)                                                                                                                  cTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacggc                                                 >  3:277837/1‑134 (MQ=255)                                                                                                                  cTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAgccacggc                                                 <  4:277837/134‑1 (MQ=255)                                                                                                                    gggTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGc                                <  8:243689/149‑1 (MQ=255)                                                                                                                     ggTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCa                               <  6:291212/149‑1 (MQ=255)                                                                                                                         cAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAg                              >  6:337720/1‑146 (MQ=255)                                                                                                                         cAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAg                              <  5:337720/146‑1 (MQ=255)                                                                                                                             cTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACgg                       >  6:184843/1‑149 (MQ=255)                                                                                                                              ttCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGt                      >  4:182979/1‑149 (MQ=255)                                                                                                                                      gCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTggg              >  6:20865/1‑149 (MQ=255)                                                                                                                                          gAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGt            <  5:118303/147‑1 (MQ=255)                                                                                                                                           aaCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTg         <  2:8669/149‑1 (MQ=255)                                                                                                                                            aCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTgg        <  8:62655/149‑1 (MQ=255)                                                                                                                                            aCTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGATAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCtt                                  <  6:122926/123‑1 (MQ=255)                                                                                                                                                  ccgccgTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCaaaa  <  5:184843/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAA  >  NC_002947/4065662‑4065954

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