Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	318,967	G→A	7.8%	R259C (CGC→TGC)	PP_0264 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	318,967	0	G	A	7.8%	186.4 / 4.7	64	R259C (CGC→TGC)	PP_0264	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/3);  total (31/33)
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.97e-01

GCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGA  >  NC_002947/318826‑319113                                                                                                                                              |                                                                                                                                                   gccgTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTgg                                                                                                                                             <  7:25947/149‑1 (MQ=255)   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:209098/1‑149 (MQ=255)                                      cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                        >  4:292135/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                                                                                                 <  3:240085/140‑1 (MQ=255)                                                      ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTg                                                                                                 >  4:240085/1‑140 (MQ=255)                                                      ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGa                                                                                          >  1:142817/1‑147 (MQ=255)                                                         aTGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAggg                                                                                     >  3:61285/1‑149 (MQ=255)                                                          tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                              <  6:228782/107‑1 (MQ=255)                                                          tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                              >  5:228782/1‑107 (MQ=255)                                                              cATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGa                                                                                >  8:207389/1‑149 (MQ=255)                                                               aTTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGAt                                                                               >  3:157994/1‑149 (MQ=255)                                                                  aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                             <  4:289102/148‑1 (MQ=255)                                                                  aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                             <  5:155934/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                                >  5:132933/1‑124 (MQ=255)                                                                                                                                       aTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                <  6:132933/124‑1 (MQ=255)                                                                                                                                            gcgGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCAGCGAAGTCAGGa  >  3:187303/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                   GCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGA  >  NC_002947/318826‑319113

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