Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,428,461	T→G	7.1%	D102A (GAC→GCC)	sohB ←	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,428,461	0	T	G	7.1%	174.8 / 4.9	56	D102A (GAC→GCC)	sohB	peptidase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCA  >  NC_002947/4428315‑4428602                                                                                                                                                   |                                                                                                                                               gcgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCg                                                                                                                                              >  6:202385/1‑149 (MQ=255) gcgcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCg                                                                                                                                              <  3:10272/149‑1 (MQ=255)   gcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGc                                                                                                                                             >  4:336302/1‑148 (MQ=255)          gCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGACATCaaa                                                                                                                                      <  5:181244/148‑1 (MQ=255)            gACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                   <  7:263954/149‑1 (MQ=255)             aCGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                   >  4:205744/1‑148 (MQ=255)              cGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGACATCAAAgcctc                                                                                                                                 <  4:70140/149‑5 (MQ=255)                 aGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                              <  1:327741/149‑1 (MQ=255)                     aaGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGCGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑aa                                                                                                                         >  6:272206/1‑149 (MQ=255)                        ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATag                                                                                                                      <  4:21836/149‑2 (MQ=255)                        ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATAc                                                                                                                      <  3:26285/149‑1 (MQ=255)                        ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGAACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATAc                                                                                                                      <  5:296092/149‑1 (MQ=255)                               tgtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGctc                                                                                                               >  3:4566/1‑149 (MQ=255)                                gtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGctct                                                                                                              >  4:98151/1‑149 (MQ=255)                                   cACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGt                                                                                                           >  5:161109/1‑149 (MQ=255)                                   cACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGt                                                                                                           <  4:4566/149‑1 (MQ=255)                                    aCCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTc                                                                                                          >  3:241499/1‑149 (MQ=255)                                    aCCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTc                                                                                                          <  3:98151/149‑1 (MQ=255)                                      cAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCt                                                                                                        >  2:127639/1‑149 (MQ=255)                                              cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGTCTCTTGTCCTCGGCCTTg                                                                                                <  6:161109/149‑1 (MQ=255)                                              cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTg                                                                                                >  7:308083/1‑149 (MQ=255)                                               gccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGc                                                                                               <  4:241499/149‑1 (MQ=255)                                                        cAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTtcttct                                                                                      >  4:67246/1‑149 (MQ=255)                                                                    cACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCg                                                                          >  6:50102/1‑149 (MQ=255)                                                                     aCTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGc                                                                         >  3:126599/1‑149 (MQ=255)                                                                     aCTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGc                                                                         >  5:284444/1‑149 (MQ=255)                                                                          gTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCGTGCCCTTCTTCTGCTGTTTTTCCGCCTTgg                                                                    <  5:272206/149‑1 (MQ=255)                                                                           tCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGATTTTCCGCCTTGGc                                                                   <  1:127639/149‑1 (MQ=255)                                                                            cACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGcc                                                                  >  3:316430/1‑149 (MQ=255)                                                                             aCGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGctc                                                                                                               <  6:9745/103‑1 (MQ=255)                                                                             aCGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGctc                                                                                                               >  5:9745/1‑103 (MQ=255)                                                                              cgcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctg                                                                <  6:20169/149‑1 (MQ=255)                                                                                  cgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCtgctgc                                                            <  8:308083/149‑1 (MQ=255)                                                                                        cgagcgTCAGCAGTGCGGTGATTTCGGTGCGCAGGCTTTCGCGGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTgcgc                                                      <  5:294522/149‑1 (MQ=255)                                                                                                agcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                              >  6:263040/1‑149 (MQ=255)                                                                                                agcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                              >  2:184807/1‑149 (MQ=255)                                                                                                       cGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       >  5:154558/1‑149 (MQ=255)                                                                                                       cGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       <  5:50102/149‑1 (MQ=255)                                                                                                           gATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCt                                   >  7:33958/1‑149 (MQ=255)                                                                                                             tttCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCtgc                                 >  2:8816/1‑147 (MQ=255)                                                                                                                cGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTc                               <  8:33958/148‑1 (MQ=255)                                                                                                                 gTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCg                              <  6:154558/148‑1 (MQ=255)                                                                                                                     cgcAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAg                         <  3:205744/149‑1 (MQ=255)                                                                                                                           cTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTCAGCTGGGCCTGCTCGACAAGGcccag                   <  1:8816/149‑4 (MQ=255)                                                                                                                            tttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGAc                  <  3:336302/149‑1 (MQ=255)                                                                                                                            tttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGAc                  >  7:94728/1‑149 (MQ=255)                                                                                                                             ttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTCACCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACt                 >  3:291676/1‑149 (MQ=255)                                                                                                                              tCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTc                >  1:247379/1‑149 (MQ=255)                                                                                                                                 gTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAg             <  5:263040/149‑1 (MQ=255)                                                                                                                                      cGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGgcgc         >  1:220686/1‑149 (MQ=255)                                                                                                                                       gAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGta                               >  1:212150/1‑124 (MQ=255)                                                                                                                                       gAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTc                               <  2:212150/125‑1 (MQ=255)                                                                                                                                           cTTTGATGGCGACATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTCCGCAGGGCCTTGAGCTGTGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    >  5:340161/1‑148 (MQ=255)                                                                                                                                            tttGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    >  3:410/1‑147 (MQ=255)                                                                                                                                            tttGATGGCGACATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       >  2:119434/1‑112 (MQ=255)                                                                                                                                              tGATGTCGACATCAAAGTCCATTACG‑AATACCCGGCTCTTGTCCTCTGCCTTGCCCTTCTTCTGATTTTTTTCCGCCTTTGCCTGCTGCTTGCGCAGTGCCTTGATATTTTCATTTTCGAGCAGGCCGGACTCCATTCGCTCACGCa  <  7:164340/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GCGCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACG‑AATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCA  >  NC_002947/4428315‑4428602

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