Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,472,713	C→G	6.0%	S38S (TCG→TCC)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,472,713	0	C	G	6.0%	216.5 / 3.9	67	S38S (TCG→TCC)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base G (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCGAGG  >  NC_002947/1472567‑1472860                                                                                                                                                   |                                                                                                                                                    gtgtCATTACCCATTTTGTGCACCTTGCGGCGCAGGTTCGGGTAGACTCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGGCCATGGAGTAGCGCCAGATGCCGACTGCGCCGAGCAGGAAGATGAAGTGCTTCGACTGCGa                                                                                                                                                   <  8:228488/149‑1 (MQ=255)   ggCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCCTGACGCGGAGGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGt                                                                                                                                                 >  5:282609/1‑149 (MQ=255)          cccATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGa                                                                                                                                           <  2:2783/148‑1 (MQ=255)          cccATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGa                                                                                                                                           <  2:402759/148‑1 (MQ=255)          cccATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGa                                                                                                                                           <  2:8817/148‑1 (MQ=255)            cATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGatat                                                                                                                                        >  8:386161/1‑149 (MQ=255)                  cTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTg                                                                                                                                  >  5:196220/1‑149 (MQ=255)                   tGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTgg                                                                                                                                 >  2:174747/1‑149 (MQ=255)                     cACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTgg                                                                                                                                 >  6:402454/1‑147 (MQ=255)                     cACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGc                                                                                                                                <  1:143244/148‑1 (MQ=255)                      aCCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGaa                                                                                                                                             >  8:357342/1‑134 (MQ=255)                      aCCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGaa                                                                                                                                             <  7:357342/134‑1 (MQ=255)                           gcggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggc                                                                                                                           <  5:319696/147‑1 (MQ=255)                            gggcgCAGGAACGGTCACACG‑CGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggccagg                                                                                                                       <  6:312498/148‑1 (MQ=255)                             ggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCa                                                                                                                                      <  5:12112/134‑1 (MQ=255)                             ggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCa                                                                                                                                      >  6:12112/1‑134 (MQ=255)                             ggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggccagg                                                                                                                       >  6:263775/1‑149 (MQ=255)                              gcgcAGGTACGGGTACACGCCGTAAAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggccaggg                                                                                                                      <  6:177588/149‑1 (MQ=255)                                gcAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCg                                                                                                                    <  8:72080/149‑1 (MQ=255)                                  aGGTACGGGTACACGGCGTAGAGGAAGATCATGCCTCGGATGAAATTGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAATTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGAt                                                                                                                  <  4:178048/149‑1 (MQ=255)                                       cGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGc                                                                                                              <  5:263775/148‑1 (MQ=255)                                        gggTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCAt                                                                                                            <  6:196220/149‑1 (MQ=255)                                         ggTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCAt                                                                                                            >  6:183589/1‑148 (MQ=255)                                          gTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCAt                                                                                                            >  1:214682/1‑147 (MQ=255)                                          gTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGa                                                                                                          <  5:402454/149‑1 (MQ=255)                                             cacGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCACATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGc                                                                                                        <  7:317882/148‑1 (MQ=255)                                             cacGCCGTAGAGGAAGAGCATGCCGCGGATGAAATGGGGGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCa                                                                                                       <  5:364357/149‑1 (MQ=255)                                               cGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAgg                                                                                                     <  1:174747/149‑1 (MQ=255)                                                gCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGc                                                                                                    >  1:204303/1‑149 (MQ=255)                                                    tAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCAt                                                                                                <  2:204303/149‑1 (MQ=255)                                                                aTGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCATATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGTCGAAg                                                                                                                                            <  5:49466/93‑1 (MQ=255)                                                                aTGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGTCGAAg                                                                                                                                            >  6:49466/1‑93 (MQ=255)                                                                aTGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGTCGAAg                                                                                                                                            >  6:43455/1‑93 (MQ=255)                                                                aTGACGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGTCGAAg                                                                                                                                            <  5:43455/93‑1 (MQ=255)                                                                   ccGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCAcc                                                                                  <  4:124174/148‑1 (MQ=255)                                                                    cgcgGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCg                                                                                <  4:72482/149‑1 (MQ=255)                                                                    cgcgGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCg                                                                                >  4:369459/1‑149 (MQ=255)                                                                     gcgGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCgg                                                                               <  6:282609/149‑1 (MQ=255)                                                                       ggATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGgcg                                                                             <  3:369459/149‑1 (MQ=255)                                                                             aaTGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTg                                                                        >  2:244613/1‑148 (MQ=255)                                                                                 ggTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACa                                                                   <  7:386161/149‑1 (MQ=255)                                                                                     gcgcCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCa                                                                <  5:384388/148‑1 (MQ=255)                                                                                           aTGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAACCACCCGGCGCACTGCAACAGCACTGTCTg                                                         >  6:154674/1‑149 (MQ=255)                                                                                              gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggccagg                                                                                                                       <  6:294054/84‑1 (MQ=255)                                                                                              gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGcagggccagg                                                                                                                       >  5:294054/1‑84 (MQ=255)                                                                                                gTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGcc                                                    <  5:183589/149‑1 (MQ=255)                                                                                                     gcCAGATGCCGACTGCGACGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGATCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTg                                                >  5:275167/1‑148 (MQ=255)                                                                                                       cAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGCTCTCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCAt                                             >  5:282914/1‑149 (MQ=255)                                                                                                          aTGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTc                                                           >  3:109920/1‑132 (MQ=255)                                                                                                          aTGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTc                                                           <  4:109920/132‑1 (MQ=255)                                                                                                          aTGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAACCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGtt                                          <  5:154674/149‑1 (MQ=255)                                                                                                             ccGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCa                                        <  5:207863/148‑1 (MQ=255)                                                                                                              cGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAg                                       >  2:252452/1‑148 (MQ=255)                                                                                                                 cTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCAt                                             <  8:178507/139‑1 (MQ=255)                                                                                                                 cTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCAt                                             >  7:178507/1‑139 (MQ=255)                                                                                                                    cgcCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGtt                                          <  3:372771/139‑1 (MQ=255)                                                                                                                    cgcCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGtt                                          >  4:372771/1‑139 (MQ=255)                                                                                                                       cgcccAGGAAGATGAAGTGCTGCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCATCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCga                             <  6:282914/146‑1 (MQ=255)                                                                                                                           cAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAAt                         >  5:41298/1‑149 (MQ=255)                                                                                                                           cAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAAt                         >  4:349356/1‑149 (MQ=255)                                                                                                                              gAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCg                      >  2:134396/1‑149 (MQ=255)                                                                                                                                aGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCggg                    <  2:214682/149‑1 (MQ=255)                                                                                                                                    gAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAg                >  7:161310/1‑149 (MQ=255)                                                                                                                                    gAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAg                >  5:303494/1‑149 (MQ=255)                                                                                                                                                cTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCga    <  3:349356/149‑1 (MQ=255)                                                                                                                                                cTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCga    <  6:41298/149‑1 (MQ=255)                                                                                                                                                 tGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCgag   >  6:397694/1‑149 (MQ=255)                                                                                                                                                   cGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATcc                                  >  1:51049/1‑116 (MQ=255)                                                                                                                                                   cGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATcc                                  <  2:51049/116‑1 (MQ=255)                                                                                                                                                   cGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCgagg  <  1:244613/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTGGGCAGCGAGG  >  NC_002947/1472567‑1472860

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