Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,137	G→C	5.9%	S267T (AGT→ACT)	cysNC →	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,491,137	0	G	C	5.9%	212.6 / 3.9	68	S267T (AGT→ACT)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base C (2/2);  total (30/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 = 9.36e-01

GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATC  >  NC_002947/1490991‑1491285                                                                                                                                                   |                                                                                                                                                      gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                                                                                                                                     >  1:200621/1‑149 (MQ=255) gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                                                                                                                                     <  6:218242/149‑1 (MQ=255)  ccGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCATCCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                                    >  4:98342/1‑149 (MQ=255)  ccGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                                    >  8:233423/1‑149 (MQ=255)       gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tgg                                                                                                                                              <  5:168593/149‑1 (MQ=255)       gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                               <  3:356913/148‑1 (MQ=255)        tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCTGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                             <  7:129039/149‑1 (MQ=255)        tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                             >  5:383209/1‑149 (MQ=255)        tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                             <  1:263975/149‑1 (MQ=255)          cGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGc                                                                                                                                           <  6:243301/149‑1 (MQ=255)                cccGACGCCGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAg                                                                                                                                     <  1:352466/149‑1 (MQ=255)                     cGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCCTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                                <  8:176090/149‑1 (MQ=255)                        tGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                                >  4:188413/1‑146 (MQ=255)                        tGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                                >  2:309751/1‑146 (MQ=255)                              aaaTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                       >  5:217528/1‑149 (MQ=255)                               aaTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCg                                                                                                                      >  1:55628/1‑149 (MQ=255)                                         aCCGTCGAGGTTTCGGCTGACCGCAACTTCACCGCCCTGCGTTTCCCGGGGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGGCCGTGGCG‑GGGTGCGCAAAGGTGACGAGATTGTCGTGCTGCCggc                                                                                                            >  8:58590/1‑147 (MQ=255)                                           cGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          <  3:310008/149‑1 (MQ=255)                                           cGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          <  2:312969/149‑1 (MQ=255)                                                 ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          >  4:288914/1‑143 (MQ=255)                                                 ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          <  3:288914/143‑1 (MQ=255)                                                 ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAaga                                                                                                    <  2:11137/149‑1 (MQ=255)                                                  gTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAagag                                                                                                   >  5:95284/1‑149 (MQ=255)                                                  gTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGGGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAagag                                                                                                   <  5:199148/149‑1 (MQ=255)                                                     tCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCATCCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGagca                                                                                                 <  3:98342/148‑1 (MQ=255)                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAg                                                                                                                                     <  3:144332/109‑1 (MQ=255)                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAg                                                                                                                                     >  4:144332/1‑109 (MQ=255)                                                         cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                            <  3:188413/149‑1 (MQ=255)                                                         cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                            >  6:93891/1‑149 (MQ=255)                                                         cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAACGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                            >  5:320980/1‑149 (MQ=255)                                                         cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTAAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                            <  5:15795/149‑1 (MQ=255)                                                           gACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGt                                                                                          >  3:60963/1‑149 (MQ=255)                                                            aCCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTc                                                                                         <  6:320980/149‑1 (MQ=255)                                                               gCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAg                                                                                      <  6:383209/149‑1 (MQ=255)                                                                cAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                       >  2:30855/1‑115 (MQ=255)                                                                cAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAg                                                                                      >  2:361247/1‑148 (MQ=255)                                                                cAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGAGGCG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                       <  1:30855/115‑1 (MQ=255)                                                                 aaCTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTc                                                                                    >  5:195731/1‑149 (MQ=255)                                                                                 gTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGaa                                                                    <  2:55628/149‑1 (MQ=255)                                                                                  tttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                                   <  8:102285/149‑1 (MQ=255)                                                                                  tttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                                   <  7:233423/149‑1 (MQ=255)                                                                                  tttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                                   <  1:210247/149‑1 (MQ=255)                                                                                     cccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGt                                                                 <  6:195731/148‑1 (MQ=255)                                                                                      ccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGa                                                               <  4:124388/149‑1 (MQ=255)                                                                                        ggTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                             <  6:217528/149‑1 (MQ=255)                                                                                                      aCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGACAACGCTGGc                                               >  4:95201/1‑149 (MQ=255)                                                                                                      aCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGc                                               >  5:119203/1‑149 (MQ=255)                                                                                                       ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGcc                                              >  1:263813/1‑149 (MQ=255)                                                                                                             cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggcc                                        >  7:290986/1‑149 (MQ=255)                                                                                                             cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                         >  2:213142/1‑148 (MQ=255)                                                                                                                 cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCcaggccaggc                                    >  7:170919/1‑149 (MQ=255)                                                                                                                     aCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                 <  6:119203/148‑1 (MQ=255)                                                                                                                          cgcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGAccc                            <  7:13462/148‑1 (MQ=255)                                                                                                                          cgcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCt                           >  7:239045/1‑149 (MQ=255)                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTg                          <  1:361247/149‑1 (MQ=255)                                                                                                                                 ctttgtgCACCA‑CGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATg                    >  8:239998/7‑149 (MQ=255)                                                                                                                                  tttgtgCACCA‑CGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATg                    >  8:229977/6‑148 (MQ=255)                                                                                                                                  cgcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATg                    <  2:263813/148‑1 (MQ=255)                                                                                                                                   gcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTTGTGCAGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATgg                   <  1:309751/148‑1 (MQ=255)                                                                                                                                   gcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATgg                   <  6:95284/148‑1 (MQ=255)                                                                                                                                     gtgCACCA‑CGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa                 <  7:239998/146‑1 (MQ=255)                                                                                                                                          ccATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc              >  1:42846/1‑146 (MQ=255)                                                                                                                                           cATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                                >  1:162571/1‑111 (MQ=255)                                                                                                                                           cATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                                <  2:162571/111‑1 (MQ=255)                                                                                                                                            cacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgaga          <  7:229977/147‑1 (MQ=255)                                                                                                                                               gCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGa      <  2:200621/149‑1 (MQ=255)                                                                                                                                                 cAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa                 >  1:136603/1‑136 (MQ=255)                                                                                                                                                 cAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa                 <  2:136603/136‑1 (MQ=255)                                                                                                                                                   gTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAtc  >  4:195877/1‑149 (MQ=255)                                                                                                                                                   gTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCAATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAtc  >  5:118403/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATC  >  NC_002947/1490991‑1491285

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