Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,443,382	Δ1 bp	5.9%	coding (1446/1989 nt)	PP_4781 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,443,382	0	G	.	5.9%	260.5 / 10.9	85	coding (1446/1989 nt)	PP_4781	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (40/39);  new base . (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

CCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGTTCATCGG  >  NC_002947/5443233‑5443531                                                                                                                                                       |                                                                                                                                                      ccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccg                                                                                                                                                       <  2:266321/150‑1 (MQ=255)   gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                                     >  2:526272/1‑150 (MQ=255)   gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                                     <  1:620370/150‑1 (MQ=255)   gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                                     <  1:1655420/150‑1 (MQ=255)      aCCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                                  >  1:1712134/1‑150 (MQ=255)        cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCa                                                                                                                                                >  2:1241989/1‑150 (MQ=255)                  gcccGTCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAg                                                                                                                                      >  2:1425986/3‑150 (MQ=255)                   tCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAgg                                                                                                                                     <  2:1632607/150‑1 (MQ=255)                    ccGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTTCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAggg                                                                                                                                    <  2:1374541/150‑1 (MQ=255)                      ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTa                                                                                                                                  >  2:1432994/1‑150 (MQ=255)                           gTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGT‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAAt                                                                                                                             <  1:1240639/150‑1 (MQ=255)                            tCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATg                                                                                                                            <  2:1308584/150‑1 (MQ=255)                            tCGGCCGGAATGTCCAGGCTCAGGAACAGGCCTTTGCCCAGCGCCGGCCCGGCACAACCAGCGCGCAAATCCTGCAACCGGCCGCCCAGGCTAAACGGGGCCGCCGGCGCACGCAGG‑CGGCCGGCCTGCAGAGCGAAGAGGGTAaaaagg                                                                                                                            >  2:802769/1‑148 (MQ=255)                              ggCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGTGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGCGTAAGAATggg                                                                                                                          <  2:318052/150‑2 (MQ=255)                                          cAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATc                                                                                                              >  1:1527252/1‑150 (MQ=255)                                          cAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATc                                                                                                              >  1:880323/1‑150 (MQ=255)                                           aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                             >  2:1251688/1‑150 (MQ=255)                                           aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                             >  2:545495/1‑150 (MQ=255)                                           aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCAGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                             >  1:1260388/1‑150 (MQ=255)                                            ggCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCca                                                                                                            <  1:690878/150‑1 (MQ=255)                                            ggCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCca                                                                                                            <  2:447585/150‑1 (MQ=255)                                              cTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCTTGCAACAGGGCGCACAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCcatc                                                                                                          <  2:1038603/150‑1 (MQ=255)                                                     aCAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCtt                                                                                                   >  1:667225/1‑150 (MQ=255)                                                         gCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCACAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATc                                                                                                                       <  2:610454/126‑1 (MQ=255)                                                         gCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCACAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATc                                                                                                                       >  1:610454/1‑126 (MQ=255)                                                          ccTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCGCCCGCTGGGGCGGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGACGGCCTGCAGTTCGGACAGGGTAAGAATGGCCTCGACCATATCCATCATGCCTTGGGCa                                                                                              <  2:880323/150‑1 (MQ=255)                                                           cTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAg                                                                                             >  2:1206144/1‑150 (MQ=255)                                                            tttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                                  <  1:1730242/96‑1 (MQ=255)                                                            tttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                                  >  2:1730242/1‑96 (MQ=255)                                                             ttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCTCGCAGG‑CGGCCGGCCTCAATTGCGGACAGGGGAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGAc                                                                                           <  1:1179547/150‑1 (MQ=255)                                                             ttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTTGGCAGAc                                                                                           <  1:1206144/150‑1 (MQ=255)                                                                 ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                       >  2:956958/1‑150 (MQ=255)                                                                 ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                       <  1:1358512/150‑1 (MQ=255)                                                                 ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                       >  1:1678410/1‑150 (MQ=255)                                                                  cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                                      >  2:292536/1‑150 (MQ=255)                                                                  cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGAACACTGAACGGTGCCGGCTACGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                                      >  1:987510/1‑150 (MQ=255)                                                                        cTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                                <  1:707838/150‑1 (MQ=255)                                                                         tGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgtg                                                                               <  2:2647/150‑1 (MQ=255)                                                                            ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGAcc                                                                                                                   <  1:1287135/111‑1 (MQ=255)                                                                            ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGAcc                                                                                                                   >  2:1287135/1‑111 (MQ=255)                                                                            ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                            <  1:1282316/150‑1 (MQ=255)                                                                             cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                                                >  2:26230/1‑129 (MQ=255)                                                                             cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                                                <  1:26230/129‑1 (MQ=255)                                                                                cATAACCGGCGCGCAAACTCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGt                                                                        <  1:96415/150‑1 (MQ=255)                                                                                 aTAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTa                                                                       <  1:526272/150‑1 (MQ=255)                                                                                     ccGGCGCGCAAACCCTGAAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGa                                                                   <  1:1677561/150‑1 (MQ=255)                                                                                       ggCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                            >  2:1288298/1‑139 (MQ=255)                                                                                       ggCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                            <  1:1288298/139‑1 (MQ=255)                                                                                            actgaaccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                                            >  2:1380154/7‑150 (MQ=255)                                                                                              aaaCCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCa                                                          <  2:1442730/150‑1 (MQ=255)                                                                                                  ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                      >  1:489868/1‑150 (MQ=255)                                                                                                  ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                      >  2:1829483/1‑150 (MQ=255)                                                                                                  ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                      >  1:166054/1‑150 (MQ=255)                                                                                                  ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGAACATAGGACAGGGTAACGACGGCGTCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                      >  1:1176959/1‑150 (MQ=255)                                                                                                    tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                                    <  1:1425986/150‑1 (MQ=255)                                                                                                    tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                                    >  2:1717455/1‑150 (MQ=255)                                                                                                    tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCAAAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                                    >  1:771893/1‑150 (MQ=255)                                                                                                             gcgcgcAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAgg                                           <  2:232073/150‑1 (MQ=255)                                                                                                             gcgcgcAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAgg                                           <  2:676366/150‑1 (MQ=255)                                                                                                                cgcAGGCTGAACGGTGCCGGCTGCGCCTGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtg                                        >  2:437208/1‑150 (MQ=255)                                                                                                                  cAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                                        >  1:184202/1‑132 (MQ=255)                                                                                                                  cAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                                        <  2:184202/132‑1 (MQ=255)                                                                                                                    ggCTGAACGGTGCCGGCTGCGCGCGCCGG‑CGGCCGGCCTGCAGTTCGGAGAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCa                                    <  2:576996/150‑1 (MQ=255)                                                                                                                     gCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                                   >  2:1396951/1‑150 (MQ=255)                                                                                                                     gCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                                   >  2:1396897/1‑150 (MQ=255)                                                                                                                      cTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATc                                  <  2:1305910/150‑1 (MQ=255)                                                                                                                       tGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCa                                 <  2:113745/150‑1 (MQ=255)                                                                                                                       tGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCa                                 <  1:783053/150‑1 (MQ=255)                                                                                                                              tGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgag                          >  1:1524076/1‑150 (MQ=255)                                                                                                                                 cGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                       <  2:489868/150‑1 (MQ=255)                                                                                                                                 cGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                       >  1:31684/1‑150 (MQ=255)                                                                                                                                  ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtgt                                       <  1:1896361/133‑1 (MQ=255)                                                                                                                                  ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtgt                                       >  2:1896361/1‑133 (MQ=255)                                                                                                                                  ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcgag                      >  2:547409/1‑150 (MQ=255)                                                                                                                                    cTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGcc                    <  2:1176959/150‑1 (MQ=255)                                                                                                                                      gcgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                                     <  1:91465/99‑1 (MQ=255)                                                                                                                                      gcgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                                     >  2:91465/1‑99 (MQ=255)                                                                                                                                       gaactgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt                 >  2:1720770/6‑150 (MQ=255)                                                                                                                                       cgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt                 <  1:1174235/150‑1 (MQ=255)                                                                                                                                          ctgCAGGCCGGCC‑GCCTGCGGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                                <  2:1112534/82‑1 (MQ=255)                                                                                                                                          ctgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                       <  2:394950/139‑1 (MQ=255)                                                                                                                                          ctgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                       >  1:394950/3‑141 (MQ=255)                                                                                                                                          ctgCAGG‑CCGGCCGCCTGCGGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                                >  1:1112534/3‑84 (MQ=255)                                                                                                                                           tgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACt             >  1:767703/2‑150 (MQ=255)                                                                                                                                           cGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACt             <  1:1713384/150‑1 (MQ=255)                                                                                                                                            gCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTc            >  2:1510745/1‑150 (MQ=255)                                                                                                                                            gCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGCTGTGCATCAGCTCGACCGAGCCAACCACTc            >  2:330050/1‑150 (MQ=255)                                                                                                                                              aGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCg          >  2:992188/1‑150 (MQ=255)                                                                                                                                                g‑cggcCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGtt        <  2:31684/150‑1 (MQ=255)                                                                                                                                                       ggccTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGTTCATCgg  <  2:771893/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGTTCATCGG  >  NC_002947/5443233‑5443531

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