Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,438,254	A→T	10.5%	K457M (AAG→ATG) ‡	PP_3932 →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,438,254	0	A	T	10.5%	161.3 / 10.2	57	K457M (AAG→ATG) ‡	PP_3932	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base T (2/4);  total (23/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 = 3.95e-01

CCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGC  >  NC_002947/4438110‑4438399                                                                                                                                                 |                                                                                                                                                  ccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGccc                                                                                                                                              <  1:10514/150‑1 (MQ=255)     aCAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggta                                                                                                                                          <  1:384907/150‑3 (MQ=255)       aGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGCGGTGTGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtaca                                                                                                                                        <  1:691993/150‑5 (MQ=255)         gACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGc                                                                                                                                      >  2:557138/1‑150 (MQ=255)            cTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggta                                                                                                                                          <  1:667028/143‑3 (MQ=255)            cTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggta                                                                                                                                          >  2:667028/1‑141 (MQ=255)            cTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAAc                                                                                                                                   >  1:222131/1‑150 (MQ=255)               cACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACcgc                                                                                                                                <  2:150583/150‑1 (MQ=255)                 cGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGt                                                                                                                              >  2:425885/1‑150 (MQ=255)                  gTGACACAGGCCTACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                             >  1:159813/1‑150 (MQ=255)                       acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggta                                                                                                                                          >  2:523325/1‑130 (MQ=255)                       acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggta                                                                                                                                          <  1:523325/132‑3 (MQ=255)                       acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGt                                                                                                                              <  2:361210/144‑1 (MQ=255)                       acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGt                                                                                                                              >  1:361210/1‑144 (MQ=255)                       acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                        <  2:375587/150‑1 (MQ=255)                        caGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAg                                                                                                                       <  2:196479/150‑1 (MQ=255)                           gCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                    >  2:464155/1‑150 (MQ=255)                               aCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGGCGAACCGGCCGAAg                                                                                                                >  2:152585/1‑150 (MQ=255)                                  aGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGcac                                                                                                             >  1:514556/1‑150 (MQ=255)                                        tGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGcc                                                                                                       >  2:835786/1‑150 (MQ=255)                                            gATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGt                                                                                                   >  1:695071/1‑150 (MQ=255)                                               tACCGTGAGTATCGGTTTTGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTACGAGGCGGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGAt                                                                                                <  1:480514/150‑1 (MQ=255)                                                 ccGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCt                                                                                              >  1:778927/1‑150 (MQ=255)                                                  cGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCTTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTg                                                                                             >  1:304747/1‑150 (MQ=255)                                                     gAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAAc                                                                                          <  2:222131/150‑1 (MQ=255)                                                                     ggCTGCGACGCAGACATGCCTGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGt                                                                          <  1:219566/150‑1 (MQ=255)                                                                       cTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTcc                                                                        >  1:717818/1‑150 (MQ=255)                                                                         gCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCa                                                                      <  2:695071/150‑1 (MQ=255)                                                                         gCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCa                                                                      <  2:686475/150‑1 (MQ=255)                                                                          cGACGCAGACATGCCTGCCCCCCGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCCCAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCAc                                                                     <  1:421636/150‑1 (MQ=255)                                                                          cGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCAc                                                                     <  2:159813/150‑1 (MQ=255)                                                                                   cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                            <  2:778927/150‑1 (MQ=255)                                                                                   cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                            <  2:798790/150‑1 (MQ=255)                                                                                    aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGGCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTgg                                                           >  1:692278/1‑150 (MQ=255)                                                                                     tGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGc                                                          <  2:233993/150‑1 (MQ=255)                                                                                         ggCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATc                                                      <  2:63088/150‑1 (MQ=255)                                                                                          gCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCa                                                     <  2:514556/150‑1 (MQ=255)                                                                                           cccccAGCCTCGACGAACTGCTGGCCAGTGCCGCCCAGGCGCTGTACCGGGCCAATGCCCGGGGCCGCAACCGCGTCGAACAGGACGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAg                                                    <  1:313072/150‑1 (MQ=255)                                                                                             cccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGca                                                  <  1:455798/150‑1 (MQ=255)                                                                                                         gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCa                                      >  2:735762/1‑150 (MQ=255)                                                                                                            ctgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCa                                   <  2:717818/150‑1 (MQ=255)                                                                                                             tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                            <  2:167636/124‑1 (MQ=255)                                                                                                             tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                            >  1:167636/1‑124 (MQ=255)                                                                                                              gctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAgg                                 <  1:425885/150‑1 (MQ=255)                                                                                                                tgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGaa                               <  1:661365/150‑1 (MQ=255)                                                                                                                  gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                   >  2:177539/1‑112 (MQ=255)                                                                                                                  gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                   <  1:177539/112‑1 (MQ=255)                                                                                                                             ccAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACa                  >  1:111161/1‑150 (MQ=255)                                                                                                                              cAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                       >  2:153015/1‑144 (MQ=255)                                                                                                                              cAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                       <  1:153015/144‑1 (MQ=255)                                                                                                                                   gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCCCCTGGTGCTCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACCACAGCACCAGGCACAGCCCGCGGCCAAACCCg            >  2:389365/1‑150 (MQ=255)                                                                                                                                     tGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGc          <  2:111161/150‑1 (MQ=255)                                                                                                                                         ccGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCtt                                                             <  2:147041/95‑1 (MQ=255)                                                                                                                                         ccGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCtt                                                             >  1:147041/1‑95 (MQ=255)                                                                                                                                         ccGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAccc      <  2:468524/150‑1 (MQ=255)                                                                                                                                           gggCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCtt    <  1:835786/150‑1 (MQ=255)                                                                                                                                             gCCAAGGCCCGGGGCCACAACCGCGTCACCCCAATCCAAGCCAACCACGTGGTTACCTAAACTATCCTGCCAGCAGCTCCCACGACAGCTTGTCGATCAGCACACACAAAAGCACCAGCAACAGCCCGCGGACAAACCCGGCACCCTTGc  <  2:270280/150‑1 (MQ=255)                                                                                                                                                aaGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGc                                                          <  2:484046/91‑1 (MQ=255)                                                                                                                                                aaGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGc                                                          >  1:484046/1‑91 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGC  >  NC_002947/4438110‑4438399

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