Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,030,751	C→T	6.9%	intergenic (+37/‑458)	PP_4433 → / → dadA‑I	amino acid MFS transporter/D‑amino acid dehydrogenase small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,030,751	0	C	T	6.9%	167.8 / 2.2	58	intergenic (+37/‑458)	PP_4433/dadA‑I	amino acid MFS transporter/D‑amino acid dehydrogenase small subunit
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

CCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGCC  >  NC_002947/5030608‑5030898                                                                                                                                                |                                                                                                                                                     ccTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCAc                                                                                                                                                >  6:65705/1‑149 (MQ=255)   tactacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAAGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCAcc                                                                                                                                               >  2:76016/1‑148 (MQ=255)     ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTTATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCAccc                                                                                                                                              <  7:19565/147‑1 (MQ=255)     ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTA‑GCACCCAg                                                                                                                                            <  8:32299/149‑1 (MQ=255)      tacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAgg                                                                                                                                           >  7:231902/1‑149 (MQ=255)      tacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAg                                                                                                                                            >  1:69937/1‑148 (MQ=255)             tGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGc                                                                                                                                    >  8:54765/1‑149 (MQ=255)               gTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCg                                                                                                                                   >  7:108915/1‑148 (MQ=255)                    aTCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATc                                                                                                                                                      <  5:120163/125‑1 (MQ=255)                    aTCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATc                                                                                                                                                      >  6:120163/1‑125 (MQ=255)                    aTCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTc                                                                                                                             <  5:65705/149‑1 (MQ=255)                     tCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTc                                                                                                                             >  1:310010/1‑148 (MQ=255)                               gtattattGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGTTTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCg                                                                                                                   <  6:19/146‑1 (MQ=255)                                gtttattGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTATTAGGCAGATGGTATATAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCgg                                                                                                                  <  1:218983/146‑1 (MQ=255)                                         gttgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGAc                                                                                                        >  5:54993/1‑149 (MQ=255)                                         gttgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGC‑AAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCACGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGAc                                                                                                        <  1:71925/149‑1 (MQ=255)                                             ttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATtcctc                                                                                                                     <  2:316400/132‑1 (MQ=255)                                             ttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATtcctc                                                                                                                     >  1:316400/1‑132 (MQ=255)                                             ttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCgg                                                                                                     <  6:250558/148‑1 (MQ=255)                                               cTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCt                                                                                                   <  8:108915/148‑1 (MQ=255)                                                   gtgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTg                                                                                                                                     >  6:90129/1‑110 (MQ=255)                                                   gtgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTg                                                                                                                                     <  5:90129/110‑1 (MQ=255)                                                   gtgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTc                                                                                              >  8:314901/1‑149 (MQ=255)                                                       tCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGTCCGGCTACGTCGat                                                                                           <  1:130402/148‑1 (MQ=255)                                                          agaCTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATaa                                                                                       <  8:231902/149‑1 (MQ=255)                                                                aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcag                                                                                   >  7:144546/1‑147 (MQ=255)                                                                   gTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAAAAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCAATCCTCCGGGCAAGCCGACCGCCTCCGGCGATATAAGCAGCAgcg                                                                              >  2:242135/1‑149 (MQ=255)                                                                      gCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTc                                                                           >  8:288446/1‑149 (MQ=255)                                                                        tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                 >  3:7624/1‑141 (MQ=255)                                                                        tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                 <  4:7624/141‑1 (MQ=255)                                                                          gATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCAAAGCTCCATTCCTCCGGGCaa                                                                                                              <  1:179832/110‑1 (MQ=255)                                                                          gATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCAAAGCTCCATTCCTCCGGGCaa                                                                                                              >  2:179832/1‑110 (MQ=255)                                                                                  gggATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCttt                                                               <  2:310010/149‑1 (MQ=255)                                                                                         aaTACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCt                                                        >  4:60392/1‑149 (MQ=255)                                                                                          aTACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCtt                                                       <  4:296031/149‑1 (MQ=255)                                                                                              cTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTAtct                                                   <  7:314901/149‑1 (MQ=255)                                                                                               tGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTAtctc                                                  >  5:131532/1‑149 (MQ=255)                                                                                                   gCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCC‑‑TTTATTCTTATCTCACctct                                            >  3:297313/1‑149 (MQ=255)                                                                                                        aTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTgg                                         <  8:144546/149‑1 (MQ=255)                                                                                                                   ggCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGccc                               <  7:288446/148‑1 (MQ=255)                                                                                                                    gCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCt                              >  6:116769/1‑148 (MQ=255)                                                                                                                     cAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            <  3:60392/149‑1 (MQ=255)                                                                                                                      aGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGAt                           <  2:18811/149‑1 (MQ=255)                                                                                                                        aTGGTTTAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATACGCAGCTGCGCTCCAGCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGAt                           <  1:242135/147‑1 (MQ=255)                                                                                                                         tGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCC‑‑TTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAg                      <  4:297313/149‑1 (MQ=255)                                                                                                                           gTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACa                       <  7:54765/148‑1 (MQ=255)                                                                                                                           gTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACa                       <  5:116769/148‑1 (MQ=255)                                                                                                                            tataGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATGTCACCTCTTGGCATCTGGCCCTGATAACAgt                     >  8:215972/1‑149 (MQ=255)                                                                                                                                  gCCTGGGTCCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTAATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCtt                <  7:215972/148‑1 (MQ=255)                                                                                                                                    cTGGGTGCTGATTA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGTGCTCCACCCTTCCTTTTTAtt                                                          <  2:28501/104‑1 (MQ=255)                                                                                                                                    cTGGGTGCTGATTA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGTGCTCCACCCTTCCTTTTTAtt                                                          >  1:28501/1‑104 (MQ=255)                                                                                                                                    cTGGGTGCTGATTA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGc             >  6:198202/1‑149 (MQ=255)                                                                                                                                    cTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGc             >  6:110433/1‑149 (MQ=255)                                                                                                                                            tAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCtt     >  5:107633/1‑149 (MQ=255)                                                                                                                                             aaTCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTg    >  2:47175/1‑149 (MQ=255)                                                                                                                                              aTCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGc             >  6:280791/1‑139 (MQ=255)                                                                                                                                              aTCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGc             <  5:280791/139‑1 (MQ=255)                                                                                                                                              aTCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTg    <  1:76016/148‑1 (MQ=255)                                                                                                                                              aTCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGc   >  3:215362/1‑149 (MQ=255)                                                                                                                                               tCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGcc  >  8:255510/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCA‑GCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGCC  >  NC_002947/5030608‑5030898

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