Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	413,107	C→G	8.3%	A686A (GCC→GCG)	glnE →	glutamate‑ammonia‑ligase adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	413,107	0	C	G	8.3%	180.1 / 8.1	60	A686A (GCC→GCG)	glnE	glutamate‑ammonia‑ligase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base G (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAA  >  NC_002947/412959‑413255                                                                                                                                                     |                                                                                                                                                      gcCATTTGAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  1:422014/150‑1 (MQ=255) gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  2:772808/150‑1 (MQ=255) gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      >  2:970986/1‑150 (MQ=255) gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  1:285112/150‑1 (MQ=255)          aGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtc                                                                                                                                             <  1:797188/150‑3 (MQ=255)            cTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtctg                                                                                                                                           <  2:907977/150‑5 (MQ=255)                  cacaGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAg                                                                                                                                     >  1:544378/1‑150 (MQ=255)                               tGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCa                                                                                                                       <  2:544378/150‑1 (MQ=255)                                       tgggAAATGAGCGGCAACCTGCCGTGGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGc                                                                                                               >  2:923492/3‑150 (MQ=255)                                       tCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCGGGTGCTGGCCCTGGCCTGGCGCCAGCCCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑c                                                                                                                           <  2:478042/138‑1 (MQ=255)                                       tCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑c                                                                                                                           >  1:478042/1‑138 (MQ=255)                                                    gCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtca                                                                                                  >  2:41078/1‑150 (MQ=255)                                                     cAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtcat                                                                                                 <  1:970986/150‑1 (MQ=255)                                                       aCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtcatca                                                                                               >  1:256964/1‑150 (MQ=255)                                                        ccTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTcatcat                                                                                              <  1:477136/150‑1 (MQ=255)                                                             cGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                         <  1:41078/150‑1 (MQ=255)                                                              gTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCtgtg                                                                                                                <  1:1047746/126‑1 (MQ=255)                                                              gTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCtgtg                                                                                                                >  2:1047746/1‑126 (MQ=255)                                                                   tgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCGGTGCGACCCGGGTTTCATCATTATCGGTTACg                                                                                   >  2:353473/1‑150 (MQ=255)                                                                       aGTCAGCGCCTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCaa                                                                               <  1:288851/150‑1 (MQ=255)                                                                       aGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCaa                                                                               >  2:425251/1‑150 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCa                                                                                                                       <  1:1064432/108‑1 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCa                                                                                                                       >  2:1064432/1‑108 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCtgt                                                                                                                 >  2:244062/1‑114 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCtgt                                                                                                                 <  1:244062/114‑1 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtca                                                                                                  >  2:691075/1‑129 (MQ=255)                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtca                                                                                                  <  1:691075/129‑1 (MQ=255)                                                                             cGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATggg                                                                         <  1:425251/150‑1 (MQ=255)                                                                              gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCa                                                                                                                       <  1:749590/103‑1 (MQ=255)                                                                              gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCa                                                                                                                       >  2:749590/1‑103 (MQ=255)                                                                              gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATggga                                                                        >  1:806637/1‑149 (MQ=255)                                                                              gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggc                                                                        >  1:69244/1‑150 (MQ=255)                                                                                 tACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggcggc                                                                     <  2:69244/150‑1 (MQ=255)                                                                                           ggctggctGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                           <  2:642212/150‑1 (MQ=255)                                                                                            gctggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGACAGCCTGTGCGACCCg                                                                                                         >  1:900913/1‑103 (MQ=255)                                                                                            gctggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGACAGCCTGTGCGACCCg                                                                                                         <  2:900913/103‑1 (MQ=255)                                                                                              tggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGcc                                                        >  1:724632/1‑150 (MQ=255)                                                                                                  tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACgg                                                    >  1:712674/1‑150 (MQ=255)                                                                                                  tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACgg                                                    >  1:636418/1‑150 (MQ=255)                                                                                                      gCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCg                                                >  1:560175/1‑150 (MQ=255)                                                                                                         aTCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac                                             >  2:1074778/1‑150 (MQ=255)                                                                                                              cGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgga                                        >  2:780581/1‑150 (MQ=255)                                                                                                                    ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtc                                                                                                   <  2:436830/85‑1 (MQ=255)                                                                                                                    ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTtc                                                                                                   >  1:436830/1‑85 (MQ=255)                                                                                                                        ctggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTc                              >  2:982319/1‑150 (MQ=255)                                                                                                                         tggccctgggtGGGCGCTAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGC‑AGGCCGTGGGTGTGGCCCGGGTTCACTCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCa                             <  1:990862/150‑1 (MQ=255)                                                                                                                          ggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCAt                            <  2:712674/150‑1 (MQ=255)                                                                                                                          ggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACATGGACCTGGTCTTCAt                            <  1:923492/150‑1 (MQ=255)                                                                                                                           gccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATc                           <  1:101745/150‑1 (MQ=255)                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacct                                           <  2:1070920/133‑1 (MQ=255)                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacct                                           >  1:1070920/1‑133 (MQ=255)                                                                                                                             cctggcctggCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGacctggacctgg                                   <  2:642900/140‑1 (MQ=255)                                                                                                                             cctggcctggCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGacctggacctgg                                   >  1:642900/1‑140 (MQ=255)                                                                                                                               tggcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATg                       >  1:923198/1‑150 (MQ=255)                                                                                                                                       cgatctaCCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGAt               <  2:510857/144‑1 (MQ=255)                                                                                                                                        gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATc              <  1:780581/150‑1 (MQ=255)                                                                                                                                         ccAGACCGTGGCCCGCCACGGCCGACCCAAGCGCAGCGA‑CGGCGGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATcc             >  2:253181/1‑150 (MQ=255)                                                                                                                                          ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGCCCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCg            >  1:207204/5‑150 (MQ=255)                                                                                                                                          ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                           <  2:1173060/99‑1 (MQ=255)                                                                                                                                          ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                           >  1:1173060/5‑103 (MQ=255)                                                                                                                                            gACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg                                    <  1:975896/124‑1 (MQ=255)                                                                                                                                            gACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg                                    >  2:975896/1‑124 (MQ=255)                                                                                                                                              ccGTGGCCCGCCACGGCCAACCCAAGCGAAGCGATGGGCCGCC‑GTAAAGGGCATGATGCTTCATTATCTATAGCTGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATcagatcg        >  1:510857/1‑144 (MQ=255)                                                                                                                                                   gCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCggca   >  2:58728/1‑148 (MQ=255)                                                                                                                                                    cccGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGaaa  >  1:501421/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGA‑CGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAA  >  NC_002947/412959‑413255

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