Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,065,810	C→G	8.2%	T50T (ACG→ACC) ‡	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,810	0	C	G	8.2%	191.3 / 6.4	61	T50T (ACG→ACC) ‡	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 = 6.19e-01

ATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGT  >  NC_002947/4065662‑4065956                                                                                                                                                      |                                                                                                                                                   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)                                                                                                                                                    gccgTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGt  <  8:290875/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGT  >  NC_002947/4065662‑4065956

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