Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,753,541	T→G	16.4%	intergenic (+46/‑320)	glnA → / → glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,753,541	0	T	G	16.4%	115.9 / 21.5	60	intergenic (+46/‑320)	glnA/glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base G (5/5);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CGACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCT  >  NC_002947/5753393‑5753688                                                                                                                                                     |                                                                                                                                                    cgacgaCTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGcctccctc                                                                                                                                                    >  1:1213945/1‑150 (MQ=255)   acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                  <  1:1131559/150‑1 (MQ=255)   acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGATGATGATAACATAGCGGAAGCGCGTCGAGTACGGGATCAAAACGGCCTCCCTAgg                                                                                                                                                  >  2:12040/1‑150 (MQ=255)        ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGc                                                                                                                                             >  1:657494/1‑150 (MQ=255)        ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAAGACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGCGGGAGGc                                                                                                                                             >  1:731533/1‑150 (MQ=255)        ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACATTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                             >  2:248539/1‑150 (MQ=255)             cGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                        <  2:657494/150‑1 (MQ=255)               aTGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttttt                                                                                                                                      <  2:1122208/150‑1 (MQ=255)                tGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTg                                                                                                                                     <  2:513873/150‑1 (MQ=255)                 gCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGt                                                                                                                                    >  1:1039030/1‑150 (MQ=255)                 gCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttttgat                                                                                                                                    <  2:327812/150‑4 (MQ=255)                  ccTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGtt                                                                                                                                   <  1:811690/150‑1 (MQ=255)                   cTTCATCTAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGttt                                                                                                                                  >  2:855447/1‑150 (MQ=255)                   cTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTAGACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGttt                                                                                                                                  >  2:384098/1‑150 (MQ=255)                    ttCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtttt                                                                                                                                       <  2:267386/144‑1 (MQ=255)                    ttCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtttt                                                                                                                                       >  1:267386/1‑144 (MQ=255)                              tGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCtt                                                                                                                       <  1:855447/150‑1 (MQ=255)                                aaGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGc                                                                                                                     <  2:1213945/150‑1 (MQ=255)                                    gCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGaa                                                                                                                 <  2:1039030/150‑1 (MQ=255)                                     cGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCAATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttttt                                                                                                                                      <  2:439323/128‑1 (MQ=255)                                     cGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCAATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttttt                                                                                                                                      >  1:439323/1‑128 (MQ=255)                                         gaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgcc                                                                                                            >  2:774050/1‑150 (MQ=255)                                           agaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccgc                                                                                                          <  2:865453/150‑1 (MQ=255)                                                           aCCTTCGTACACCCGCTGGAATACGAGNTGTACTACAGCTNCTGATCNGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACg                                                                                          >  1:1040045/1‑150 (MQ=255)                                                             cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttttgat                                                                                                                                    <  1:1028772/106‑4 (MQ=39)                                                             cTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttttgat                                                                                                                                    >  2:1028772/1‑103 (MQ=39)                                                                 gTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGc                                                                                    >  2:1081351/1‑150 (MQ=255)                                                                   acacCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                  >  2:320154/1‑150 (MQ=255)                                                                         gCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                             <  2:941920/85‑1 (MQ=255)                                                                         gCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                             >  1:941920/1‑85 (MQ=255)                                                                          cTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCATCCCTCGGGATGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                           <  1:559761/150‑1 (MQ=255)                                                                             gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTg                                                                        <  1:774050/150‑1 (MQ=255)                                                                               aTACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCa                                                                      >  2:321732/1‑150 (MQ=255)                                                                                  gcaGCTGTAGTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATg                                                                                     >  1:1162696/3‑132 (MQ=255)                                                                                  gcaGCTGTAGTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATg                                                                                     <  2:1162696/130‑1 (MQ=255)                                                                                    cgCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctct                                                                                 >  2:878936/2‑134 (MQ=255)                                                                                    aGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctct                                                                                 <  1:878936/134‑1 (MQ=255)                                                                                       tGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTca                                                              >  2:667379/1‑150 (MQ=255)                                                                                         tactacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTcaca                                                            <  1:320154/150‑1 (MQ=255)                                                                                          actacAGCTGCTGATCTGATCGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCacac                                                           <  1:321732/150‑1 (MQ=255)                                                                                               aGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                           >  2:529206/1‑129 (MQ=255)                                                                                               aGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                           <  1:529206/129‑1 (MQ=255)                                                                                                   cctgatctgatTGGCGCCGCTCGCCGATGACATGACCACCACTGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCACAAGCGCCGCTGATTCGCATGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTc                                                  <  1:12040/149‑1 (MQ=255)                                                                                                                gcgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCaa                                     >  1:272439/1‑150 (MQ=255)                                                                                                                 cgcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAAt                                    <  1:667379/150‑1 (MQ=255)                                                                                                                   ccTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCt                                  >  2:312976/1‑150 (MQ=255)                                                                                                                     tCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTgg                                <  2:1040045/150‑1 (MQ=255)                                                                                                                     tCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTgg                                >  1:872288/1‑150 (MQ=255)                                                                                                                         gcCGATTACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCa                            <  1:384098/150‑1 (MQ=255)                                                                                                                         gcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCt                                                                                                <  1:313353/82‑1 (MQ=255)                                                                                                                         gcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCt                                                                                                >  2:313353/1‑82 (MQ=255)                                                                                                                         gcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGt                                                                                              >  2:717711/1‑84 (MQ=255)                                                                                                                         gcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGt                                                                                              <  1:717711/84‑1 (MQ=255)                                                                                                                          ccGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAg                           <  1:312976/150‑1 (MQ=255)                                                                                                                                  aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTaa                                                                  <  2:103839/103‑1 (MQ=255)                                                                                                                                  aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTaa                                                                  >  1:103839/1‑103 (MQ=255)                                                                                                                                   tGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGt                  <  2:731533/150‑1 (MQ=255)                                                                                                                                   tGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGt                  <  2:872288/150‑1 (MQ=255)                                                                                                                                        aaaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAc             >  2:1079103/1‑150 (MQ=255)                                                                                                                                          aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAt                                                                 >  2:457647/1‑96 (MQ=255)                                                                                                                                          aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAt                                                                 <  1:457647/96‑1 (MQ=255)                                                                                                                                           ccgatctccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCCATCCGGTTCAGGTTAAACGTTTCACCaa          <  1:814753/145‑1 (MQ=255)                                                                                                                                                ctccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCa              >  2:571323/1‑141 (MQ=255)                                                                                                                                                ctccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCa              <  1:571323/141‑1 (MQ=255)                                                                                                                                                 tccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACttt    <  2:1216754/150‑1 (MQ=255)                                                                                                                                                  ccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc   <  1:1081351/150‑1 (MQ=255)                                                                                                                                                   cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCCTGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGagatcgg  >  2:814753/1‑144 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCT  >  NC_002947/5753393‑5753688

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