Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,065,802	T→A	7.6%	Q53L (CAG→CTG)	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,802	0	T	A	7.6%	196.7 / 3.4	65	Q53L (CAG→CTG)	PP_3582	RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCTTCTGATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGG  >  NC_002947/4065654‑4065948                                                                                                                                                      |                                                                                                                                                   ggCTTCTGATTGCACCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACt                                                                                                                                                    <  5:205911/149‑1 (MQ=255)         aTTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCgccgcc                                                                                                                                            >  4:107120/1‑149 (MQ=255)         aTTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCgccgcc                                                                                                                                            >  4:84115/1‑149 (MQ=255)             cgcCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGc                                                                                                                                        >  7:261092/1‑149 (MQ=255)              gcCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTgc                                                                                                                                                  >  6:160214/1‑138 (MQ=255)              gcCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTgc                                                                                                                                                  <  5:160214/138‑1 (MQ=255)                   ggTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCa                                                                                                                                  <  6:161446/149‑1 (MQ=255)                         gTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagca                                                                                                                              <  1:210740/147‑5 (MQ=255)                             ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                        >  7:70516/1‑149 (MQ=255)                             ctgcGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCg                                                                                                                        <  3:107120/149‑1 (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                                                                                                                     <  1:151163/140‑1 (MQ=255)                                         tctTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCggtg                                                                                                                     >  2:151163/1‑140 (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                                                                                                    >  1:253068/1‑149 (MQ=255)                                                 cAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAggg                                                                                                    >  4:265640/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                                                                                  >  4:239986/1‑149 (MQ=255)                                                                   aTCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAgg                                                                                  >  1:8669/1‑149 (MQ=255)                                                                       ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                              >  6:70770/1‑149 (MQ=255)                                                                       ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCAc                                                                              >  5:34316/1‑149 (MQ=255)                                                                       ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                               <  1:72777/100‑4 (MQ=38)                                                                       ccGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaagc                                                                                                                               >  2:72777/1‑97 (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                                                                                                                                 >  2:223162/1‑85 (MQ=37)                                                                                   caccaGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACAGCACGGCGGCGCAGTTCaa                                                                                                                                 <  1:223162/86‑2 (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)                                                                                                                                                      |                                                                                                                                                   GGCTTCTGATTGCGCCACGGTCTGGTTGCTGCGGTTGAGCTCTTCGACCAGCGCATTGAGCCCAGCATCGCCGTAGATTTCCCACCAGGCGCCGCGGGCGA‑TGGCGTCGGACGGCGAGGCCTGGGTCCAGCCTTCGGCGTGCTTGAACTGCGCCGCCGTGCTGTTCAGCTCAGGGCGGTGGTAGTCGGGGCTCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGG  >  NC_002947/4065654‑4065948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑