Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,625,964	C→T	7.2%	R513H (CGC→CAC)	PP_4944 ←	carbamoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,625,964	0	C	T	7.2%	161.9 / 2.1	55	R513H (CGC→CAC)	PP_4944	carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/2);  total (30/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 = 1.00e+00

CGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTCGAG  >  NC_002947/5625819‑5626109                                                                                                                                                  |                                                                                                                                                   cGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACAACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTACCGCGGGTTGTACTCGCGCt                                                                                                                                                 >  2:74005/1‑149 (MQ=255)          gaaCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCt                                                                                                                                        >  4:367562/1‑149 (MQ=255)          gaaCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCt                                                                                                                                        <  3:56352/149‑1 (MQ=255)          gaaCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCt                                                                                                                                        <  3:320709/149‑1 (MQ=255)            aCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTgg                                                                                                                                      <  3:317327/149‑1 (MQ=255)            aCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTgg                                                                                                                                      <  7:32490/149‑1 (MQ=255)                     gggCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGa‑a                                                                                                                            <  7:276951/149‑1 (MQ=255)                        cATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTg                                                                                                                         <  1:74005/149‑1 (MQ=255)                         aTCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGc                                                                                                                        >  5:121246/1‑149 (MQ=255)                              tcggCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTcc                                                                                                                   >  4:170695/1‑149 (MQ=255)                                ggCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTc                                                                                                                 <  8:319879/149‑1 (MQ=255)                                   gAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGCTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTg                                                                                                              >  5:203025/1‑149 (MQ=255)                                     gCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGca                                                                                                            <  2:176130/149‑1 (MQ=255)                                       aGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGCTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGcacc                                                                                                          <  6:203025/149‑1 (MQ=255)                                           catcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACtt                                                                                                      >  4:33255/1‑149 (MQ=255)                                            atcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTc                                                                                                     >  6:202535/1‑149 (MQ=255)                                            atcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTc                                                                                                     >  6:134880/1‑149 (MQ=255)                                                gTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGc                                                                                                 <  6:121246/149‑1 (MQ=255)                                                     ccacGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCGGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCAcgcg                                                                                            >  3:275495/1‑149 (MQ=255)                                                      cacGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCgg                                                                                           <  3:367562/149‑1 (MQ=255)                                                         gacgGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGtt                                                                                        >  7:252190/1‑149 (MQ=255)                                                            ggTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTc                                                                                     <  8:252190/149‑1 (MQ=255)                                                                  gCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGTGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCActct                                                                               >  7:184724/1‑149 (MQ=255)                                                                  gCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGTGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCActct                                                                               >  7:179459/1‑149 (MQ=255)                                                                     aGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACGTCCGGCCCGCGGGGTTTCCCctcgt                                                                              >  8:235333/1‑145 (MQ=255)                                                                          ttCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGAc                                                                       >  6:131995/1‑149 (MQ=255)                                                                           tCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACt                                                                      >  2:309057/1‑149 (MQ=255)                                                                             aGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCgcgc                                                                                                                                                  <  2:350336/72‑1 (MQ=255)                                                                                 acacacCGTTGCCGGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaag                                                                <  4:275495/149‑1 (MQ=255)                                                                                     acCGTTGCCAGTCAGGTCTTCCAGTCCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTCAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtg                                                             >  2:64435/1‑148 (MQ=255)                                                                                            ccAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCAt                                                     >  4:83413/1‑149 (MQ=255)                                                                                             cAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCAt                                                     >  3:108934/1‑148 (MQ=255)                                                                                              aGTCAGGTCTTCCAGTGCCATCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTACCACTCTTCAGCGACTTCGAAGGTGAAGTTCATg                                                    >  5:123455/1‑148 (MQ=255)                                                                                                 cAGGGCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAAc                                                 <  3:33255/148‑1 (MQ=255)                                                                                                  aGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACgg                                               >  2:3058/1‑149 (MQ=255)                                                                                                     tCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGTGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACGACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGgcgc                                            <  8:184724/149‑1 (MQ=255)                                                                                                     tCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGTGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGgcgc                                            <  8:179459/149‑1 (MQ=255)                                                                                                       ttCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACAAGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCATCGAGCACAACTTCCGGCACGCCGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGAGCgg                                          >  5:224652/1‑149 (MQ=255)                                                                                                            gTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATgg                                     >  2:60226/1‑149 (MQ=255)                                                                                                                 ttCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGAACCACTTACGGCACGCGGGTTTTCAACTCTTCAGCGACTTCGAAGGTGAAGGTAATGAACGGCGCGGGATGGTCGAc                                >  7:176227/1‑149 (MQ=255)                                                                                                                       atGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGat                          >  6:115418/1‑149 (MQ=255)                                                                                                                       atGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGa                           <  5:134880/148‑1 (MQ=255)                                                                                                                       atGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGa                           <  3:170695/148‑1 (MQ=255)                                                                                                                       atGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGa                           <  3:83413/148‑1 (MQ=255)                                                                                                                            gtagtaGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGATAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGTTCGACCTTGatcata                      >  1:261533/1‑148 (MQ=255)                                                                                                                               gtagCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCg                  >  2:75644/1‑149 (MQ=255)                                                                                                                               gtagCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCg                  >  2:34289/1‑149 (MQ=255)                                                                                                                                  gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCgggg               <  4:108934/149‑1 (MQ=255)                                                                                                                                  gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCgggg               <  1:349831/149‑1 (MQ=255)                                                                                                                                   cgcgGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGc              <  4:331279/149‑1 (MQ=255)                                                                                                                                           gTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGTACGGCGCGGGATGGTCGACCTTGa                           >  1:132381/1‑128 (MQ=255)                                                                                                                                           gTACTCGCGCTTCAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGTACGGCGCGGGATGGTCGACCTTGa                           <  2:132381/128‑1 (MQ=255)                                                                                                                                              cTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTCGa   <  1:309057/149‑1 (MQ=255)                                                                                                                                               tCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCtt                             <  3:63370/122‑1 (MQ=255)                                                                                                                                               tCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTCGAg  >  8:39388/1‑149 (MQ=255)                                                                                                                                               tCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCAGGAACGGCGCGGGATGGTCGACCtt                             >  4:63370/1‑122 (MQ=255)                                                                                                                                                 gcgcTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTc     <  6:120952/144‑1 (MQ=255)                                                                                                                                                 gcgcTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTc     >  5:120952/1‑144 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGA‑AGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACAGTGTCGAG  >  NC_002947/5625819‑5626109

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