Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,059,329	G→A	6.3%	R203H (CGC→CAC)	PP_2671 →	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,059,329	0	G	A	6.3%	187.7 / 2.1	63	R203H (CGC→CAC)	PP_2671	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/2);  total (32/31)
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.16e-01

CATCCGCTGGGCGGTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTGGCAGGCATGCGTGCC  >  NC_002947/3059186‑3059471                                                                                                                                                |                                                                                                                                                cATCCGCTGGGCGGTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATgg                                                                                                                                            >  5:23618/1‑149 (MQ=255)            cGGTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTg                                                                                                                                 <  7:279353/149‑1 (MQ=255)             ggTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTgg                                                                                                                                <  2:277723/149‑1 (MQ=255)              gTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTTCACCAGGTCTCAGGCGGGCTATTGACGGGACGCCTGCGCGAGGGCGGCCTGCCGGAAGCTCGCCATTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGa                                                                                                                               <  5:324906/149‑1 (MQ=255)              gTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGa                                                                                                                               <  5:330132/149‑1 (MQ=255)               tACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAg                                                                                                                              >  2:106941/1‑149 (MQ=255)                  ggcggGCCATGGGCCTGCTTGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAg                                                                                                                           <  1:112650/149‑1 (MQ=255)                     gggCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCt                                                                                                                        >  3:100989/1‑149 (MQ=255)                              gCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGa                                                                                                                >  7:249495/1‑148 (MQ=255)                                 tgctgGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAAt                                                                                                            >  6:50729/1‑149 (MQ=255)                                   ctgGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGc                                                                                                          >  8:165133/1‑149 (MQ=255)                                           gctgctGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGa                                                                                                 <  6:23618/149‑1 (MQ=255)                                                  ggggCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAgg                                                                                           <  5:50729/148‑1 (MQ=255)                                                     gCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCct                                                                                       >  1:236759/1‑149 (MQ=255)                                                       cTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCctgc                                                                                     <  8:249495/149‑1 (MQ=255)                                                          gCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTgg                                                                                  >  6:228293/1‑149 (MQ=255)                                                            aCCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCa                                                                                <  7:165133/149‑1 (MQ=255)                                                              cAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGt                                                                              <  6:5202/149‑1 (MQ=255)                                                                ggTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGc                                                                                 >  4:10158/1‑144 (MQ=255)                                                                ggTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGc                                                                                 <  3:10158/144‑1 (MQ=255)                                                                          gggCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTgcaggagcag                                                                    >  2:20969/1‑147 (MQ=255)                                                                            gCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCaggagcaggagc                                                                <  4:100989/149‑1 (MQ=255)                                                                             cAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCACATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCg                                                               >  7:131411/1‑149 (MQ=255)                                                                              aaTTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCg                                                               >  5:273396/1‑148 (MQ=255)                                                                                  gACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTaca                                                           >  1:261253/1‑148 (MQ=255)                                                                                  gACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTaca                                                           >  4:92436/1‑148 (MQ=255)                                                                                     ggCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTAcaca                                                         >  6:258661/1‑147 (MQ=255)                                                                                     ggCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTAcaca                                                         <  5:258661/147‑1 (MQ=255)                                                                                     ggCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGc                                                       >  6:55899/1‑149 (MQ=255)                                                                                       cACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTg                                                     <  6:273396/149‑1 (MQ=255)                                                                                        aCGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCctgctg                                                                                   <  2:75438/118‑1 (MQ=255)                                                                                        aCGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCctgctg                                                                                   >  1:75438/1‑118 (MQ=255)                                                                                            cTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGt                                                              <  6:226239/135‑1 (MQ=255)                                                                                            cTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGt                                                              >  5:226239/1‑135 (MQ=255)                                                                                                  gAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCt                                          >  1:41568/1‑149 (MQ=255)                                                                                                      ccggccTGCCGGTAGCGCGGCATTTGGCCAAACACTTCAACCGCAAGGGTGCAACCCTTGATGAGGCTCATGCCGACAATGCCCTAGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCa                                       <  1:106941/149‑1 (MQ=255)                                                                                                        ggccTGCCGGAAGCGCGCCAGTTGGCCAAACCCTTCACCCGCAGCGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACga                                    <  4:150142/149‑1 (MQ=255)                                                                                                          ccTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACgatg                                  <  5:228293/149‑1 (MQ=255)                                                                                                             gCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCACATGGCTGCAACCCTGGAGCAGGCTCAGGcc                                                                                                                  >  2:253209/1‑67 (MQ=255)                                                                                                             gCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCACATGGCTGCAACCCTGGAGCAGGCTCAGGcc                                                                                                                  <  1:253209/67‑1 (MQ=255)                                                                                                                 gAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGc                                        <  7:93549/136‑1 (MQ=255)                                                                                                                 gAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGc                                        >  8:93549/1‑136 (MQ=255)                                                                                                                  aaGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTggg                          <  2:236759/149‑1 (MQ=255)                                                                                                                  aaGCGCGCCAGTTGGCCAAACACTTCAACCACATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTggg                          <  8:131411/149‑1 (MQ=255)                                                                                                                      gcgcCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTg                            <  4:328369/143‑1 (MQ=255)                                                                                                                      gcgcCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTg                            >  3:328369/1‑143 (MQ=255)                                                                                                                       cgcCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATa                     >  2:175903/1‑149 (MQ=255)                                                                                                                             ttGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTggc               <  1:35753/149‑1 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGcc                                                                                         >  5:350839/1‑70 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGcc                                                                                         <  6:350839/70‑1 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTgg                           <  4:134036/132‑1 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTgg                           >  3:134036/1‑132 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTggcaggca          >  6:222264/1‑149 (MQ=255)                                                                                                                                  cAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCAGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTggcaggca          >  4:178243/1‑149 (MQ=255)                                                                                                                                   aaaCACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTGGCAGGCAt         >  3:329359/1‑149 (MQ=255)                                                                                                                                         ttCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTggcag             <  6:166311/139‑1 (MQ=255)                                                                                                                                         ttCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTggcag             >  5:166311/1‑139 (MQ=255)                                                                                                                                         ttCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTGGCAGGCATGCGTGc   >  1:206511/1‑149 (MQ=255)                                                                                                                                           cAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCAGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTGGCAGGCATGCGTGcc  <  3:178243/148‑1 (MQ=255)                                                                                                                                             aCCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTgg                           >  5:164047/1‑121 (MQ=255)                                                                                                                                             aCCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTgg                           <  6:164047/121‑1 (MQ=255)                                                                                                                                               cGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCt                                          <  5:182156/104‑1 (MQ=255)                                                                                                                                               cGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCt                                          >  6:182156/1‑104 (MQ=255)                                                                                                                                                |                                                                                                                                                CATCCGCTGGGCGGTACGGCGGGCCATGGGCCTGCTGGATGAGCTGCTGGGGGCCTTGCACCAGGTCTCAGGCGGGCAATTGACGGCACGCCTGCGCGAGGCCGGCCTGCCGGAAGCGCGCCAGTTGGCCAAACACTTCAACCGCATGGCTGCAACCCTGGAGCAGGCTCAGGCCGACAATGCCC‑AGCTGACCCAGGCCCTGCTGGCAGTGCAGGAGCAGGAGCGTACACAGCTGGCGCAGACCCTGCACGATGATCTTGGGCAATACCTGGCAGGCATGCGTGCC  >  NC_002947/3059186‑3059471

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