Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,144	G→C	7.0%	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,144	0	G	C	7.0%	174.5 / 4.8	58	S267T (AGT→ACT)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (21/32);  new base C (2/2);  total (24/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 = 6.83e-01

TGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATC  >  NC_002947/1491002‑1491286                                                                                                                                                  |                                                                                                                                                 tGGTACGCCGGCCCGACGATGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                          >  2:457377/1‑150 (MQ=255)   gTACGCCGGCACGAAACTGATATAATTCCTTGGAACCGTCTAGGTTTCGGCTGACCGCCAC‑TCATTTATCCCGCTTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCgatggt                                                                                                                                        <  1:299797/150‑6 (MQ=255)       gCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCCAAAGTTATATCGGGTACTTCCGCGGCTGCGCGGGCACCATCGCCAGTGGCG‑TGGTGcac                                                                                                                                    >  2:29585/1‑150 (MQ=255)           gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCACCAGTGGCG‑TGGTGCACAAAg                                                                                                                                <  2:1077297/150‑1 (MQ=255)             ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCCGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGt                                                                                                                              <  2:299388/150‑1 (MQ=255)               gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGa                                                                                                                            <  2:185562/150‑1 (MQ=255)                    tGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACgaga                                                                                                                       <  1:457377/150‑1 (MQ=255)                      aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGc                                                                                                                                              <  1:86615/126‑1 (MQ=255)                      aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGc                                                                                                                                              >  2:86615/1‑126 (MQ=255)                       tGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGGCGAGATTg                                                                                                                    >  1:86363/1‑150 (MQ=255)                             tCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtgc                                                                                                              >  2:726800/1‑150 (MQ=255)                                tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGca                                                                                                                                     >  2:312299/1‑124 (MQ=255)                                tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGca                                                                                                                                     <  1:312299/124‑1 (MQ=255)                                     acagTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTAC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCg                                                                                                      <  2:371774/147‑1 (MQ=255)                                     aCCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCGAtggtg                                                                                                                                       <  2:767390/117‑1 (MQ=255)                                     aCCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCGAtggtg                                                                                                                                       >  1:767390/1‑117 (MQ=255)                                      ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGACCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                     <  2:726696/150‑1 (MQ=255)                                          cGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCaa                                                                                                 <  2:139350/150‑1 (MQ=255)                                              gTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGagc                                                                                             >  1:339611/1‑150 (MQ=255)                                              gTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGagc                                                                                             <  1:835147/150‑1 (MQ=255)                                               tttCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTCCTCCCGTTCA‑T‑CGTTATGCTTCTCAACCTGAACTTCCGCTGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGagca                                                                                            >  1:861869/1‑150 (MQ=255)                                                   ggCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCATCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCg                                                                                        >  2:571473/1‑150 (MQ=255)                                                   ggCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                  >  2:17362/1‑124 (MQ=255)                                                   ggCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTc                                                                                                                  <  1:17362/124‑1 (MQ=255)                                                         ccGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCaa                                                                                 >  1:658505/1‑150 (MQ=255)                                                          cGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACGATGAACATCGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGACGTCGGGAAAGAGCAGCCG‑GGTCAAg                                                                                >  1:619288/1‑150 (MQ=255)                                                           gCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGt                                                                               <  2:861869/150‑1 (MQ=255)                                                             aaCTTCACCGA‑CCTGCGTTTC‑‑CAGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTcc                                                                             >  2:719253/1‑150 (MQ=255)                                                                ttCCCCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCCCAGCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATc                                                                          <  1:11286/150‑1 (MQ=255)                                                                ttCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATc                                                                          <  1:669403/150‑1 (MQ=255)                                                                   aCCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCAAAACAAATCTCTCAGTTCTTCTCTGGCACCATCGCAAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTc                                                                       <  1:806317/150‑1 (MQ=255)                                                                    ccGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCa                                                                      <  2:3664/150‑1 (MQ=255)                                                                             cGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAG‑AGCCGTGGTCAAGTCCATCGTCACCTACGaa                                                              <  1:29585/150‑1 (MQ=255)                                                                                     ccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGAtt                                                                                                                     >  1:19536/1‑90 (MQ=255)                                                                                     ccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGAtt                                                                                                                     <  2:19536/90‑1 (MQ=255)                                                                                     ccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGaa                                                        >  1:643620/1‑150 (MQ=255)                                                                                     ccGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGaa                                                        <  2:643620/150‑1 (MQ=255)                                                                                      cGGTGCAGTACGTCAACCGTCCCAACCGGAACGTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                       <  2:342595/150‑1 (MQ=255)                                                                                      cGGTGCAGTACGTCAAACGTGCAAACCTGATCTTCTGCGCCTCCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                       <  2:10153/150‑1 (MQ=255)                                                                                       ggTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGGACt                                                      <  1:1002824/150‑1 (MQ=255)                                                                                       ggTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACt                                                      <  2:339611/150‑1 (MQ=255)                                                                                                  tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                           >  2:108282/1‑150 (MQ=255)                                                                                                      ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                        >  2:788897/1‑54 (MQ=255)                                                                                                      ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                        <  1:788897/54‑1 (MQ=255)                                                                                                            cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggcca                                 >  2:862608/1‑150 (MQ=255)                                                                                                                cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCaa                                                                                 <  2:821219/98‑1 (MQ=255)                                                                                                                cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCaa                                                                                 >  1:821219/1‑98 (MQ=255)                                                                                                                      ggccGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGTGCAAGAGCAGCGG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGAcc                       <  2:152476/148‑1 (MQ=255)                                                                                                                          gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGa                   <  2:999624/150‑1 (MQ=255)                                                                                                                              cTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCAt               >  2:723562/1‑150 (MQ=255)                                                                                                                               ttCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATg              <  1:571473/150‑1 (MQ=255)                                                                                                                                 tttgtgCACCA‑CGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATAGAGGGATACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGa            >  1:738991/6‑150 (MQ=255)                                                                                                                                  gcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa           <  1:726800/150‑1 (MQ=255)                                                                                                                                  gcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa           <  1:553570/150‑1 (MQ=255)                                                                                                                                      gCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTc                                                                       >  2:566904/1‑86 (MQ=255)                                                                                                                                      gCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTc                                                                       <  1:566904/86‑1 (MQ=255)                                                                                                                                       gaCCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTAAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAGga      <  2:948588/149‑1 (MQ=255)                                                                                                                                           aTCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGGCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATc  >  1:1105364/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGA‑CCTGCGTTTC‑‑CCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCG‑GGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATC  >  NC_002947/1491002‑1491286

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