Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,221,731	C→T	7.6%	A286V (GCC→GTC)	PP_4601 →	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,221,731	0	C	T	7.6%	149.7 / 2.1	54	A286V (GCC→GTC)	PP_4601	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

GCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCAGTGATCGGCAG  >  NC_002947/5221584‑5221878                                                                                                                                                     |                                                                                                                                                    gcgcGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGcc                                                                                                                                                    >  3:148306/1‑149 (MQ=255) gcgcGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGcc                                                                                                                                                    >  6:4682/1‑149 (MQ=255)      aTCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGgc                                                                                                                                               <  4:148306/149‑1 (MQ=255)         aTGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGcgcg                                                                                                                                            >  5:25985/1‑149 (MQ=255)          tGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCgg                                                                                                                                           <  4:356913/149‑1 (MQ=255)             cACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCggggg                                                                                                                                        >  1:250649/1‑149 (MQ=255)                     cgCGGCCTGGCAGTCCTCATCGTGCACCTGCGCATCGTACCGGTGGCATACGGCGCCATCCACCGCCCGGGCGGGCCGGTGCGGCGCC‑GCATCGTTCCGGGCTCGGCGTTTTCCGAGCAATTCTTTTCCGACGCGCGGGGGCTGTTTCa                                                                                                                                <  2:316253/148‑1 (MQ=255)                       cGGCCTGGAAGTCCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCa                                                                                                                                <  3:11894/147‑1 (MQ=255)                       cGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAg                                                                                                                               <  2:250649/148‑1 (MQ=255)                         gCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCTCCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGc                                                                                                                             <  1:113416/149‑1 (MQ=255)                         gCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGcc                                                                                                                            <  6:324396/149‑1 (MQ=255)                         gCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGc                                                                                                                             >  1:285480/1‑148 (MQ=255)                              gAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAg                                                                                                                               >  2:205401/1‑141 (MQ=255)                              gAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAg                                                                                                                               <  1:205401/141‑1 (MQ=255)                               aaGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGa                                                                                                                      >  2:301659/1‑149 (MQ=255)                                   ccTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCGCCGCCGATGCCGAGCAATGCATTGGCGAGGCGCCGGGGCGGGTTGACGCCAGCTGATagg                                                                                                                  >  4:159889/1‑149 (MQ=255)                                       aGCGAGCACCTGCGCATCGTACCGGTGGCAGACACCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATcggcgg                                                                                                               >  7:146525/1‑148 (MQ=255)                                         cGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGTCGAGGCGCGGGGGCTGTTTCAGGCCAGc                                                                                                                         >  2:246175/1‑136 (MQ=255)                                         cGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGTCGAGGCGCGGGGGCTGTTTCAGGCCAGc                                                                                                                         >  2:251561/1‑136 (MQ=255)                                         cGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGTCGAGGCGCGGGGGCTGTTTCAGGCCAGc                                                                                                                         <  1:251561/136‑1 (MQ=255)                                         cGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGTCGAGGCGCGGGGGCTGTTTCAGGCCAGc                                                                                                                         <  1:246175/136‑1 (MQ=255)                                             cACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTg                                                                                                        >  7:79130/1‑149 (MQ=255)                                                cTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCg                                                                                                     <  3:159889/149‑1 (MQ=255)                                                         gTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAACCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCt                                                                                                         >  7:51491/1‑136 (MQ=255)                                                         gTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAACCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCt                                                                                                         <  8:51491/136‑1 (MQ=255)                                                           aCCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAAAGCGTTGCCGAGGCGCGGGGGCTGTTTCACGCCAGCTGATCGGCGGTAGCCGGGCGGGACTGGcgcg                                                                                          >  1:132579/1‑149 (MQ=255)                                                                    aGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTTCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGc                                                                                 >  6:274417/1‑149 (MQ=255)                                                                    aGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATcgg                                                                                                                  <  8:188204/116‑1 (MQ=255)                                                                    aGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATcgg                                                                                                                  >  7:188204/1‑116 (MQ=255)                                                                    aGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGc                                                                                 <  2:285480/149‑1 (MQ=255)                                                                      aCGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCg                                                                               >  8:338686/1‑149 (MQ=255)                                                                        gccgccATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGt                                                                             <  6:25985/149‑1 (MQ=255)                                                                              aTCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgctg                                                                       >  3:220149/1‑149 (MQ=255)                                                                                 cacaGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCtt                                                                    >  2:276653/1‑149 (MQ=255)                                                                                     gCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAAc                                                                <  5:4682/149‑1 (MQ=255)                                                                                      cccGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCgcggc                                                                                        >  4:118667/1‑124 (MQ=255)                                                                                      cccGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCgcggc                                                                                        <  3:118667/124‑1 (MQ=255)                                                                                          gTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgctg                                                                       >  5:388579/1‑137 (MQ=255)                                                                                          gTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgctg                                                                       <  6:388579/137‑1 (MQ=255)                                                                                               gctgctgctGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTg                                                       >  3:170104/1‑148 (MQ=255)                                                                                                               cAGCGAGCCGCGCTCGGCGATTGCCGCGGAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTTGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGGGGTCgg                                       >  7:148286/1‑149 (MQ=255)                                                                                                               cAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAAc                                                                >  5:267533/1‑124 (MQ=255)                                                                                                               cAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAAc                                                                <  6:267533/124‑1 (MQ=255)                                                                                                               cAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCgg                                       <  7:338686/149‑1 (MQ=255)                                                                                                                       cgcgCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCg                               >  3:87459/1‑149 (MQ=255)                                                                                                                         cgcTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCa                             <  1:301659/149‑1 (MQ=255)                                                                                                                           cTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTGCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAgc                           <  8:79130/149‑1 (MQ=255)                                                                                                                            tCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAgcg                          >  8:200441/1‑149 (MQ=255)                                                                                                                                 gATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCg                     >  4:255537/1‑149 (MQ=255)                                                                                                                                      ccGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACt                >  5:267125/1‑149 (MQ=255)                                                                                                                                       cGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACtt               >  7:252040/1‑149 (MQ=255)                                                                                                                                       cGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACtt               >  3:130434/1‑149 (MQ=255)                                                                                                                                         aGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCa             <  7:236985/149‑1 (MQ=255)                                                                                                                                         aGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCa             <  5:274417/149‑1 (MQ=255)                                                                                                                                           cAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCAGt           <  8:146525/149‑1 (MQ=255)                                                                                                                                                 tttGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCAGTGATCgg     >  1:142129/1‑149 (MQ=255)                                                                                                                                                     ccGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCAGTGATCGGCAg  <  4:220149/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCC‑GCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGAACTGGCGGCTGAAGGCGCTGTGGTCGGTGTAGCCGCAGCGCATCGCCACTTCAGTGATCGGCAG  >  NC_002947/5221584‑5221878

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