Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,734,779	G→C	6.8%	T295T (ACC→ACG)	hutH ←	histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,734,779	0	G	C	6.8%	183.5 / 4.4	59	T295T (ACC→ACG)	hutH	histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base C (2/2);  total (31/28)
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.50e-01

GTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTG  >  NC_002947/5734641‑5734925                                                                                                                                           |                                                                                                                                                     gTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                            >  3:222900/1‑149 (MQ=255)  ttGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                            >  6:138625/1‑148 (MQ=255)     tCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAcc                                                                                                                                        >  5:62584/1‑149 (MQ=255)      cggcggCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCTATGCCCAGTACCTCGGCGGCCTGACGCAGCTGTGTCAGGCAGGCAccc                                                                                                                                       <  6:226379/149‑1 (MQ=255)       ggcggcCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAccc                                                                                                                                       <  6:132195/148‑1 (MQ=255)       ggcggcCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCa                                                                                                                                      >  7:106153/1‑149 (MQ=255)       ggcggcCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCa                                                                                                                                      >  6:117640/1‑149 (MQ=255)                   cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCt                                                                                                                          <  4:14094/149‑1 (MQ=255)                   cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCt                                                                                                                          >  5:118951/1‑149 (MQ=255)                   cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCt                                                                                                                          >  3:205309/1‑149 (MQ=255)                    aCCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                            <  5:100799/130‑1 (MQ=255)                    aCCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                            >  6:100799/1‑130 (MQ=255)                       ggTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgca                                                                                                                      <  4:205309/149‑1 (MQ=255)                        gTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcag                                                                                                                     >  2:186/1‑149 (MQ=255)                        gTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcag                                                                                                                     <  4:222900/149‑1 (MQ=255)                         ttCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                                         <  2:95896/128‑1 (MQ=255)                         ttCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                                         >  1:95896/1‑128 (MQ=255)                             gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg‑cag                                                                                                               >  5:158585/1‑149 (MQ=255)                             gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg‑cag                                                                                                               >  7:192644/1‑149 (MQ=255)                                          cgccgGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTcc                                                                                                  >  2:3668/1‑149 (MQ=255)                                                   tcacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCt‑t                                                                                         <  4:68618/148‑1 (MQ=255)                                                     acatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTc                                                                                      >  6:93452/1‑149 (MQ=255)                                                           ccTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGtt                                                                                <  8:106153/149‑1 (MQ=255)                                                                 gCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTtgtg                                                                          <  5:98743/149‑1 (MQ=255)                                                                 gCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTtgtg                                                                          <  1:186/149‑1 (MQ=255)                                                                  cAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgc                                                                         <  6:118951/149‑1 (MQ=255)                                                                  cAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgc                                                                         <  5:117640/149‑1 (MQ=255)                                                                   aGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgc                                                                         >  8:95860/1‑148 (MQ=255)                                                                     cGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgcga                                                                       >  2:98568/1‑148 (MQ=255)                                                                        acacCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgcgacagc                                                                   <  7:329483/149‑1 (MQ=255)                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCTCTGTCGCAGTTCTTGTgcgacagcgac                                                                <  1:3668/149‑1 (MQ=255)                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgcgacagcgaca                                                               >  3:342392/1‑149 (MQ=255)                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgcgacagcgaca                                                               >  3:124932/1‑149 (MQ=255)                                                                              gCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACt                                                             >  2:277518/1‑149 (MQ=255)                                                                                 ggTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCg                                                          <  6:62584/149‑1 (MQ=255)                                                                                     gTCCGACACAGCGTTGGCTTCGATGCCCATTACCTCGGCGGCCTGACGCATCTGGGGCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCCAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                        <  1:98568/148‑1 (MQ=255)                                                                                     gTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                       <  4:342392/148‑1 (MQ=255)                                                                                          acacAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGCCGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCg                                                 >  3:10123/1‑149 (MQ=255)                                                                                          acacAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCg                                                 >  6:303661/1‑149 (MQ=255)                                                                                              aGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGcc                                               >  8:188053/1‑147 (MQ=255)                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACTGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                        <  1:261547/135‑1 (MQ=255)                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                        >  2:261547/1‑135 (MQ=255)                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagca                                          <  1:277518/149‑1 (MQ=255)                                                                                                 gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcag                                            >  5:49683/1‑147 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTg                                                                                                <  3:37851/95‑1 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTg                                                                                                >  4:37851/1‑95 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                                         <  4:124932/149‑1 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                                         <  5:303661/149‑1 (MQ=255)                                                                                                      cTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGAGAGTACTGGTACTGGACCT‑TGTCGCAGTTATTCTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcg                                     >  5:64088/1‑149 (MQ=255)                                                                                                         cGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgcg                                   <  7:95860/148‑1 (MQ=255)                                                                                                             gCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAg                               >  3:137479/1‑148 (MQ=255)                                                                                                                cAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                           >  6:96342/1‑149 (MQ=255)                                                                                                                     ccTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCcgcg                      <  6:158585/149‑1 (MQ=255)                                                                                                                         ggcggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgtgt                  <  5:93452/149‑1 (MQ=255)                                                                                                                         gcctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTc                                                  <  4:337098/113‑1 (MQ=255)                                                                                                                         gcctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTc                                                  >  3:337098/5‑117 (MQ=255)                                                                                                                           ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTc                 >  3:268483/3‑148 (MQ=255)                                                                                                                               cTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCt            >  1:132257/1‑149 (MQ=255)                                                                                                                               cTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCt            <  4:268483/149‑1 (MQ=255)                                                                                                                                           ggTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTg  <  7:188053/147‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     GTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTG  >  NC_002947/5734641‑5734925

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