Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,142:1	+A	6.5%	coding (805/1902 nt)	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,142	1	.	A	6.5%	246.5 / 11.2	62	coding (805/1902 nt)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base . (26/32);  new base A (2/2);  total (28/34)
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.65e-01

GTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTC  >  NC_002947/1490997‑1491287                                                                                                                                                   |                                                                                                                                                  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                                                                                                                                               <  1:263975/149‑1 (MQ=255)  tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                               >  5:383209/1‑149 (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                                                                                                                                  >  2:309751/1‑146 (MQ=255)                  tGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                                  >  4:188413/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                                                                                              >  6:93891/1‑149 (MQ=255)                                                   cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                              <  3:188413/149‑1 (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                                                  <  2:162571/111‑1 (MQ=255)                                                                                                                                     cATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                                  >  1:162571/1‑111 (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)                                                                                                                                                gCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAtctc  <  4:60963/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTC  >  NC_002947/1490997‑1491287

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