Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,734,772	G→C	6.6%	R298G (CGT→GGT)	hutH ←	histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,734,772	0	G	C	6.6%	196.3 / 4.4	61	R298G (CGT→GGT)	hutH	histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (2/2);  total (32/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 = 9.80e-01

GTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCT  >  NC_002947/5734641‑5734915                                                                                                                                    |                                                                                                                                                  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                                                                                                                                  >  6:100799/1‑130 (MQ=255)                    aCCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                  <  5:100799/130‑1 (MQ=255)                       ggTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgca                                                                                                            <  4:205309/149‑1 (MQ=255)                        gTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcag                                                                                                           <  4:222900/149‑1 (MQ=255)                        gTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcag                                                                                                           >  2:186/1‑149 (MQ=255)                         ttCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                               <  2:95896/128‑1 (MQ=255)                         ttCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                               >  1:95896/1‑128 (MQ=255)                             gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg‑cag                                                                                                     >  7:192644/1‑149 (MQ=255)                             gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg‑cag                                                                                                     >  5:158585/1‑149 (MQ=255)                                        gccgccGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAg                                                                                                                                                >  2:20035/1‑96 (MQ=255)                                        gccgccGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAg                                                                                                                                                <  1:20035/96‑1 (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                                                                <  1:186/149‑1 (MQ=255)                                                                 gCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTtgtg                                                                <  5:98743/149‑1 (MQ=255)                                                                  cAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgc                                                               <  5:117640/149‑1 (MQ=255)                                                                  cAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgc                                                               <  6:118951/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:124932/1‑149 (MQ=255)                                                                            cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTgcgacagcgaca                                                     >  3:342392/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                                                                                      >  4:37851/1‑95 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTg                                                                                      <  3:37851/95‑1 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                               <  5:303661/149‑1 (MQ=255)                                                                                                  ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                               <  4:124932/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                                        >  3:337098/5‑117 (MQ=255)                                                                                                                         gcctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTc                                        <  4:337098/113‑1 (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)                                                                                                                                    |                                                                                                                                                  GTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCG‑CAGCGAGTACGGGTCCTGGACCT‑TGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCT  >  NC_002947/5734641‑5734915

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