Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,562,060	A→G	7.1%	L158L (TTG→CTG)	PP_3145 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,562,060	0	A	G	7.1%	159.9 / 2.3	56	L158L (TTG→CTG)	PP_3145	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/30);  new base G (2/2);  total (24/32)
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.57e-01

GTCACGGCCACCCATGGCGGTATTCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATA  >  NC_002947/3561914‑3562208                                                                                                                                                    |                                                                                                                                                     gTCACGGCCACCCATGGCGGTATTCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAAc                                                                                                                                                     <  5:320696/148‑1 (MQ=255)                        tCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCgcgg                                                                                                                             >  7:5055/1‑149 (MQ=255)                          acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                  <  6:2842/142‑4 (MQ=255)                          acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                  >  5:2842/1‑139 (MQ=255)                                 ccGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGcc                                                                                                                    >  3:93573/1‑149 (MQ=255)                                    cgATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCg                                                                                                                 >  2:159385/1‑149 (MQ=255)                                       tAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGtt                                                                                                               >  7:131758/1‑148 (MQ=255)                                         aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  6:162985/149‑1 (MQ=255)                                         aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  7:20793/149‑1 (MQ=255)                                         aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGa                                                                                                             <  8:132144/148‑1 (MQ=255)                                          gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                       <  8:40139/137‑1 (MQ=255)                                          gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                       >  7:40139/1‑137 (MQ=255)                                              gATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                       <  1:159385‑M1/149‑1 (MQ=255)                                              gATAGGGATCACCTCGTAAGGCCGGAGGCCCTCGTTACTGGGGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                       <  7:143815‑M1/149‑1 (MQ=255)                                                   ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                 <  8:170960/118‑5 (MQ=255)                                                   ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                 >  7:170960/1‑114 (MQ=255)                                                    gATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttat                                                                                                 >  5:1028‑M1/1‑145 (MQ=255)                                                     aTCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGAt                                                                                                <  4:93573/149‑1 (MQ=255)                                                        aCCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑CCCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTc                                                                                             >  2:318708/1‑149 (MQ=255)                                                         ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                   >  3:271667/1‑142 (MQ=255)                                                         ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                   <  4:271667/142‑1 (MQ=255)                                                          cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                           <  7:11254/149‑1 (MQ=255)                                                          cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                           >  1:216897/1‑149 (MQ=255)                                                              tAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCGACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTg                                                                                          >  1:241985/1‑147 (MQ=255)                                                               aaGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                      >  6:285762/1‑149 (MQ=255)                                                                aGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCGCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                      <  5:285762/148‑1 (MQ=255)                                                                  gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            <  7:224920/124‑1 (MQ=255)                                                                  gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                            >  8:224920/1‑124 (MQ=255)                                                                         gCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCtgctgc                                                                             >  7:28615/1‑148 (MQ=255)                                                                            cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                      >  4:252782‑M1/1‑120 (MQ=255)                                                                            cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                      <  3:252782‑M1/120‑1 (MQ=255)                                                                               gTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATa                                                                       <  8:5055/148‑1 (MQ=255)                                                                                  aCTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAAccg                                                                   <  5:247132‑M2/149‑4 (MQ=255)                                                                                    tggtggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGAGAAccgcg                                                                 <  8:80284‑M2/149‑6 (MQ=255)                                                                                      gtggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGt                                                               <  3:287060/149‑1 (MQ=255)                                                                                       tggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGtt                                                              <  2:241985/149‑1 (MQ=255)                                                                                         gACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTgg                                                            <  3:153293/149‑1 (MQ=255)                                                                                          acacTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcacccgctacttca                                                           <  6:1028‑M1/149‑43 (MQ=255)                                                                                              tGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                      <  2:92087‑M1/134‑32 (MQ=255)                                                                                              tGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                      >  1:92087‑M1/1‑103 (MQ=255)                                                                                                ccGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCgg                                                      <  1:309204/148‑1 (MQ=255)                                                                                                   tGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  >  2:228375/1‑149 (MQ=255)                                                                                                        gCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTa                                             <  8:28615/149‑1 (MQ=255)                                                                                                         ccGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                               >  7:183866/1‑146 (MQ=255)                                                                                                         ccGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                               >  7:43944/1‑146 (MQ=255)                                                                                                           gggCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGaa                                          >  6:253517/1‑149 (MQ=255)                                                                                                            ggCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAAc                                         <  2:216897/149‑1 (MQ=255)                                                                                                               gc‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  >  8:101645/1‑137 (MQ=255)                                                                                                               gc‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                  <  7:101645/137‑1 (MQ=255)                                                                                                                               ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGaa                        >  5:6853/1‑148 (MQ=255)                                                                                                                               ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAac                       <  8:131758/149‑1 (MQ=255)                                                                                                                                 gggCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                     <  8:183866/149‑1 (MQ=255)                                                                                                                                     cAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTg                 <  1:228375/149‑1 (MQ=255)                                                                                                                                          cccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTa            >  8:193783/1‑149 (MQ=255)                                                                                                                                           ccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTAc           <  6:247132/149‑1 (MQ=255)                                                                                                                                                 ccAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTg     <  8:43944/149‑1 (MQ=255)                                                                                                                                                    aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata  <  1:318708/149‑1 (MQ=255)                                                                                                                                                    aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata  <  6:6853/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTCACGGCCACCCATGGCGGTATTCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATA  >  NC_002947/3561914‑3562208

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