Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,787,400	Δ1 bp	50.8%	coding (1235/1392 nt)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,787,400	0	T	.	50.8%	95.5 / 146.3	59	coding (1235/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base . (13/17);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01

AACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCAC  >  NC_002947/3787253‑3787542                                                                                                                                                    |                                                                                                                                               aaCGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATg                                                                                                                                               <  6:254524/149‑1 (MQ=255)    gTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGa                                                                                                                                           >  7:188838/1‑149 (MQ=255)       gcgcGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGAc                                                                                                                                          <  7:267281/148‑1 (MQ=255)         gcgTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGt                                                                                                                                      <  8:188838/149‑1 (MQ=255)          cgTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTg                                                                                                                                     >  1:205323/1‑149 (MQ=255)                 agaCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCcgcg                                                                                                                                <  3:163857/147‑1 (MQ=255)                   aCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgt                                                                                                                             <  2:119937/149‑1 (MQ=255)                   aCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgt                                                                                                                             <  5:113260/149‑1 (MQ=255)                     aaaCGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgtgt                                                                                                                           >  6:154848/1‑149 (MQ=255)                                 ttCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCt                                                                                                              >  8:42088/1‑149 (MQ=255)                                  tCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCt                                                                                                              >  2:298903/1‑149 (MQ=255)                                   cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGTCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTc                                                                                                             <  5:256464/148‑1 (MQ=255)                                   cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGc                                                                                                               <  5:154848/147‑1 (MQ=255)                                   cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGc                                                                                                               <  5:101944/146‑1 (MQ=255)                                   cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCt                                                                                                              <  3:247238/147‑1 (MQ=255)                                    gccgccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTc                                                                                                             <  2:253903/148‑1 (MQ=255)                                    gccgccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCAt                                                                                                           <  2:205323/149‑1 (MQ=255)                                      cgccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATg                                                                                                          >  7:328908/1‑149 (MQ=255)                                        ccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCg                                                                                                       >  3:283282/1‑149 (MQ=255)                                           gCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGc                                                                                                    <  4:132934/149‑1 (MQ=255)                                              tcctTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATg                                                                                                 >  2:127530/1‑149 (MQ=255)                                                 ttGATGGCAGTGACCACGTCGTAGTACTTCAGCTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCCACATCGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                                                                                >  4:20332/1‑148 (MQ=255)                                                        ccGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGt                                                                                       <  1:127530/149‑1 (MQ=255)                                                            gCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCg                                                                                                       <  3:264994/130‑1 (MQ=255)                                                            gCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCg                                                                                                       >  4:264994/1‑130 (MQ=255)                                                             cccACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCTTCTTcc                                                                                  <  4:283282/149‑1 (MQ=255)                                                              ccACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTcc                                                                                  <  8:328908/149‑1 (MQ=255)                                                               gaCGTCGTAGTACTTCAGGTAAAGGTCCAGGGTCATGTTGGGAGTGGCAAGCCAGCGCGTGATGGCTCGGTAGATGGTACCGACA‑GTGGAAGCGTGCAGCGGGTGTTGGCCATCATGCTCATTCCGTGCATGTACACGCCGTCTTCCTc                                                                                <  6:357913/148‑1 (MQ=255)                                                                aCGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCa                                                                               >  5:217512/1‑149 (MQ=255)                                                                 cgtcgtAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCAc                                                                              <  3:173460/149‑1 (MQ=255)                                                                            ttCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCca                                                                    >  4:24282/1‑149 (MQ=255)                                                                            ttCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCAGTCTTCCTCACGCTGTTTCca                                                                    >  5:356742/1‑149 (MQ=255)                                                                                   cAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTc                                                              >  8:32890/1‑148 (MQ=255)                                                                                      ggTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATg                                                                                                 <  2:50175/109‑1 (MQ=255)                                                                                      ggTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCAGGGCATg                                                                                                 >  1:50175/1‑109 (MQ=255)                                                                                        tCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATTGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGAt                                                        <  3:20332/149‑1 (MQ=255)                                                                                                    ttGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTTCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGaa                                                           <  3:282660/134‑1 (MQ=255)                                                                                                    ttGGGAGTGGCAAGCCAGCGCGTGCTGCCGCGGTAGATGTTACCGACATGGAGACCAGTCACGCGCGAGTTGCCCAGCATGCTCAAGCCGCTCATGTACACTCCGACTTCCTCACCCACTTTCCACACTTCGaa                                                           >  4:282660/1‑134 (MQ=255)                                                                                                     tGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATgc                                          <  4:298530/149‑1 (MQ=255)                                                                                                           tGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGc                                     <  3:24282/149‑1 (MQ=255)                                                                                                           tGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGc                                     >  6:179841/1‑149 (MQ=255)                                                                                                           tGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCg                                    >  1:53218/1‑149 (MQ=255)                                                                                                             gCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGt                                   >  5:119014/1‑149 (MQ=255)                                                                                                               aaGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTgc                                 <  1:298903/149‑1 (MQ=255)                                                                                                               aaGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTgcg                                >  6:343271/1‑149 (MQ=255)                                                                                                                   cAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAg                            >  4:163890/1‑149 (MQ=255)                                                                                                                       gcgTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGc                         >  8:267400/1‑148 (MQ=255)                                                                                                                          tGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCgaga                     <  3:337620/149‑1 (MQ=255)                                                                                                                          tGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCgag                      >  4:337620/1‑148 (MQ=255)                                                                                                                             tGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGc                  >  4:143819/1‑149 (MQ=255)                                                                                                                                  cgGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAg             <  1:275763/149‑1 (MQ=255)                                                                                                                                       gATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTc                                                                                      >  3:310469/1‑72 (MQ=255)                                                                                                                                       gATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTc                                                                                      <  4:310469/72‑1 (MQ=255)                                                                                                                                       gATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCg         <  6:119014/149‑1 (MQ=255)                                                                                                                                       gATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCg         >  2:285159/1‑149 (MQ=255)                                                                                                                                       gATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCgg        <  3:163890/149‑1 (MQ=255)                                                                                                                                         tGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGc       <  5:179841/149‑1 (MQ=255)                                                                                                                                            tACCGACATGGCGACCCGAGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGccc     >  2:127376/1‑148 (MQ=255)                                                                                                                                              ccGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCccac  <  3:143819/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               AACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCAC  >  NC_002947/3787253‑3787542

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