Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	318,975	G→A	7.8%	T256I (ACC→ATC)	PP_0264 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	318,975	0	G	A	7.8%	185.1 / 4.4	64	T256I (ACC→ATC)	PP_0264	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/3);  total (32/32)
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.30e-01

CGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCA  >  NC_002947/318828‑319119                                                                                                                                                    |                                                                                                                                                 cgTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGt                                                                                                                                                 >  7:24141/1‑149 (MQ=255)  gTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGTCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGt                                                                                                                                                 >  1:118305/1‑148 (MQ=255)     gCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGCCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             >  3:79374/1‑149 (MQ=255)      cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             <  6:298749/148‑1 (MQ=255)      cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             <  7:26553/148‑1 (MQ=255)          aCCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAgac                                                                                                                                        >  7:205179/1‑147 (MQ=255)                ccTGCCATCAACATCACTCCCAGCAGCGTACCGGCACCGATCGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCcgcga                                                                                                                                  <  8:199629/149‑2 (MQ=255)                    ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGTCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTg                                                                                                                              <  2:118305/149‑1 (MQ=255)                    ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGtt                                                                                                                               >  4:32362/1‑148 (MQ=255)                    ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCcgcgacatc                                                                                                                              <  8:205179/149‑6 (MQ=255)                                    cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  4:292135/1‑149 (MQ=255)                                    cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  4:209098/1‑149 (MQ=255)                                    cagcacCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  7:10717/1‑149 (MQ=255)                                        agggTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGc                                                                                                          <  7:256694/146‑1 (MQ=255)                                            tACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTgg                                                                                                      <  1:152785/149‑1 (MQ=255)                                             aCCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTggg                                                                                                     >  1:149889/1‑149 (MQ=255)                                              ccGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTggg                                                                                                     <  8:24141/148‑1 (MQ=255)                                               cGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTc                                                                                                   >  6:155934/1‑149 (MQ=255)                                               cGGCACCGATGGCCATTACCGTCTGCACCAAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGCGCACCAGCGGCTGTGGGGCCTGGGTc                                                                                                   >  4:223710/1‑149 (MQ=255)                                                 gCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCg                                                                                                  >  1:146236/1‑148 (MQ=255)                                                    ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTg                                                                                                       >  4:240085/1‑140 (MQ=255)                                                    ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTg                                                                                                       <  3:240085/140‑1 (MQ=255)                                                    ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGa                                                                                                >  1:142817/1‑147 (MQ=255)                                                       aTGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAggg                                                                                           >  3:61285/1‑149 (MQ=255)                                                        tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                                    >  5:228782/1‑107 (MQ=255)                                                        tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                                    <  6:228782/107‑1 (MQ=255)                                                            cATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGa                                                                                      >  8:207389/1‑149 (MQ=255)                                                             aTTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGAt                                                                                     >  3:157994/1‑149 (MQ=255)                                                                aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                                   <  5:155934/148‑1 (MQ=255)                                                                aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                                   <  4:289102/148‑1 (MQ=255)                                                                   gTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                               <  3:283749/149‑1 (MQ=255)                                                                    tCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                               >  4:165163/1‑148 (MQ=255)                                                                     cTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             >  1:163159/1‑149 (MQ=255)                                                                        cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGNCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  1:30471/146‑1 (MQ=255)                                                                        cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  3:292135/146‑1 (MQ=255)                                                                        cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  4:61285/146‑1 (MQ=255)                                                                           cGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCt                                                                       <  3:165163/149‑1 (MQ=255)                                                                            gAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTg                                                                      <  2:163159/149‑1 (MQ=255)                                                                             aGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTgg                                                                     <  4:79374/149‑1 (MQ=255)                                                                                 aTCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCt                                                                 <  1:123189/149‑1 (MQ=255)                                                                                    aaCCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCATGGCGCTCgg                                                              >  1:271360/1‑149 (MQ=255)                                                                                      ccAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTTCGGGTTGAGAACAAGCGGCTGTCGGGACTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCg                                                            >  7:217566/1‑149 (MQ=255)                                                                                       cAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGt                                                           >  8:241828/1‑149 (MQ=255)                                                                                         gCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCg                                                         <  3:209098/149‑1 (MQ=255)                                                                                            tGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTc                                                      <  2:142817/149‑1 (MQ=255)                                                                                              cgcgGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAg                                                    <  7:241828/149‑1 (MQ=255)                                                                                                  gggCGAGGATGCTGACTAGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTg                                                >  4:60146/1‑149 (MQ=255)                                                                                                   ggCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGTGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGc                                               <  7:207389/149‑1 (MQ=255)                                                                                                    gCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCg                                              <  2:149889/149‑1 (MQ=255)                                                                                                             cTGCCTTGCCAACCAGTTTCCTTTCTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGTTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      <  8:217566/148‑1 (MQ=255)                                                                                                              tGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACg                                     >  8:14855/1‑148 (MQ=255)                                                                                                                aCTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGa                                   <  7:14855/148‑1 (MQ=255)                                                                                                                    gCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgc                                <  3:60146/147‑1 (MQ=255)                                                                                                                     ccAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                              <  2:146236/148‑1 (MQ=255)                                                                                                                      cAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                              >  4:75537/1‑147 (MQ=255)                                                                                                                       aaCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCg                                                               >  8:248844/1‑113 (MQ=255)                                                                                                                       aaCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCg                                                               <  7:248844/113‑1 (MQ=255)                                                                                                                        aCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGt                          >  3:128913/1‑149 (MQ=255)                                                                                                                             tttCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGCGTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGc                     >  1:88167/1‑149 (MQ=255)                                                                                                                                    gATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGaa               <  3:32362/148‑1 (MQ=255)                                                                                                                                    gATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAg              >  2:210225/1‑149 (MQ=255)                                                                                                                                     aTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      <  6:132933/124‑1 (MQ=255)                                                                                                                                     aTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      >  5:132933/1‑124 (MQ=255)                                                                                                                                          gcgGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCAGCGAAGTCAGGa        >  3:187303/1‑149 (MQ=255)                                                                                                                                                   ggTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCa  <  4:157994/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCA  >  NC_002947/318828‑319119

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