Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,362,770	C→A	6.5%	P387Q (CCG→CAG)	phoQ →	sensor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,362,770	0	C	A	6.5%	218.5 / 3.6	70	P387Q (CCG→CAG)	phoQ	sensor protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base A (2/3);  total (34/36)
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.24e-01

GGCGCAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCTTGGC  >  NC_002947/1362624‑1362906                                                                                                                                                   |                                                                                                                                         ggcgCAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGAcgggccgg                                                                                                                                        <  2:273648/149‑1 (MQ=255)  gcgcAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGAcgggccggg                                                                                                                                       >  3:5318/1‑149 (MQ=255)  gcgcAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGAcgggccggg                                                                                                                                       >  6:8174/1‑149 (MQ=255)    gcAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCg                                                                                                                                     >  3:100747/1‑149 (MQ=255)       ggTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGc                                                                                                                                  <  7:119952/149‑1 (MQ=255)       ggTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGc                                                                                                                                  <  1:314132/149‑1 (MQ=255)           ccGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCg                                                                                                                               <  5:8174/148‑1 (MQ=255)           ccGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGt                                                                                                                              >  7:333973/1‑149 (MQ=255)            cGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTc                                                                                                                             <  7:4119/149‑1 (MQ=255)              agggAGCACGGCGCGGAGCTGGAAATGCTCGGCACCCTGCTGGAGAACGCCTACCCGCGGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCACGCCCGGCCCGTCgt                                                                                                                           <  8:324431/147‑2 (MQ=255)              aTGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGa                                                                                                                           <  8:333973/149‑1 (MQ=255)                 gAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCg                                                                                                                            >  6:371192/1‑145 (MQ=255)                 gAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCg                                                                                                                            <  5:371192/145‑1 (MQ=255)                   gCAGGGCGCTTGGTTGGAAATGCTCGGCACCCTGCTGGAGCACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAACGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCACGCCCGGCCCGTCGTcctcg                                                                                                                      <  8:119844/149‑4 (MQ=255)                   gCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGc                                                                                                                      >  6:335657/1‑149 (MQ=255)                   gCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGc                                                                                                                      <  2:134444/149‑1 (MQ=255)                     aGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGt                                                                                                                    <  7:100287/149‑1 (MQ=255)                        gcgcGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGaa                                                                                                                 >  1:281622/1‑149 (MQ=255)                                ttgGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCt                                                                                                         >  7:63354/1‑149 (MQ=255)                                   gAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCACCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTg                                                                                                        >  3:225999/1‑147 (MQ=255)                                        gCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTg                                                                                                 >  6:61150/1‑149 (MQ=255)                                            ggCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTggg                                                                                               >  1:217397/1‑147 (MQ=255)                                                  ctgctgGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTg                                                                                                 <  2:62627/139‑1 (MQ=255)                                                  ctgctgGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTg                                                                                                 >  1:62627/1‑139 (MQ=255)                                                  ctgctgGACAACGCCTACCAGATGAGCCTGGGCAAAGTGCGCGTGAGCCTGATTAAAGCGACTAGGTATCTGACCTTGAGCATCCAGGACCACGGGCCACTCGAGCAAGTCGACCACCGTGAACGCATTCTGGAGCGTGGGGAGCgg                                                                                         >  6:170139/1‑147 (MQ=255)                                                    gctgGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTg                                                                                                 >  1:352954/1‑137 (MQ=255)                                                     ctgGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTg                                                                                      <  4:5318/147‑1 (MQ=255)                                                       ggAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGa                                                                                    >  7:216409/1‑147 (MQ=255)                                                             cGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATGTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGAcagccagc                                                                            >  2:59150/1‑149 (MQ=255)                                                              gCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCa                                                                           >  4:76462/1‑149 (MQ=255)                                                                 tACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTgg                                                                                                       <  8:329185/118‑1 (MQ=255)                                                                 tACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTgg                                                                                                       >  7:329185/1‑118 (MQ=255)                                                                      gCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCACGCCCGGCCCGTCg                                                                                                                            >  2:325927/1‑92 (MQ=255)                                                                      gCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCACGCCCGGCCCGTCg                                                                                                                            <  1:325927/92‑1 (MQ=255)                                                                          aGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGgggcaggg                                                               >  8:66928/1‑149 (MQ=255)                                                                          aGCCGGGGCCAAGTCCGGGTCAGCCTGATGAAAGCGCCGGGGCATCTGCCCTTGTGCATCGAGCACGCCGGGCCGGGCGTGCCCGTCGACCAGCGGGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGgggcag                                                                 <  4:100747/147‑1 (MQ=255)                                                                            ccTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTg                                                                                                        <  3:133882/106‑1 (MQ=255)                                                                            ccTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTg                                                                                                        >  4:133882/1‑106 (MQ=255)                                                                                   cAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCtt                                                      <  2:217397/149‑1 (MQ=255)                                                                                   cAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCtt                                                      >  7:257775/1‑149 (MQ=255)                                                                                             tGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGt                                             >  5:326447/1‑148 (MQ=255)                                                                                                gCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAg                                         <  7:66928/149‑1 (MQ=255)                                                                                                gCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAg                                         <  8:216409/149‑1 (MQ=255)                                                                                                 ccTGATGAAAGCACCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAgg                                        <  4:225999/149‑1 (MQ=255)                                                                                                  cTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTgg                                                    >  5:302964/1‑136 (MQ=255)                                                                                                  cTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTgg                                                    <  6:302964/136‑1 (MQ=255)                                                                                                    gatgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAgg                                        >  2:273330/1‑146 (MQ=255)                                                                                                    gatgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAgg                                        <  1:273330/146‑1 (MQ=255)                                                                                                     atgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCt                                                       >  1:13262/1‑130 (MQ=255)                                                                                                     atgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCt                                                       <  2:13262/130‑1 (MQ=255)                                                                                                     atgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGat                                      <  5:335657/147‑1 (MQ=255)                                                                                                     atgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCAGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGatat                                    >  1:256882/1‑149 (MQ=255)                                                                                                      tgaAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATAtc                                   <  3:260948/149‑1 (MQ=255)                                                                                                        aaaGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATAtc                                   <  2:281622/147‑1 (MQ=255)                                                                                                        aaaGCACCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTgg                                                                                                <  1:312344/86‑1 (MQ=255)                                                                                                        aaaGCACCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTgg                                                                                                >  2:312344/1‑86 (MQ=255)                                                                                                         aaGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATAtcgtc                                >  8:136684/1‑149 (MQ=255)                                                                                                          aGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTgg                                                    <  5:210573/128‑1 (MQ=255)                                                                                                          aGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGCGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTgg                                                    >  6:210573/1‑128 (MQ=255)                                                                                                           gcgcCTTGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGGCCGTCGACCAGCGTGAACGTATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGTCTTGGCAGTGGTCAAGGATATCGTCGa                              <  5:170139/149‑1 (MQ=255)                                                                                                             gcCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGa                              >  7:97611/1‑147 (MQ=255)                                                                                                                  ggTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTAt                       <  3:49307/149‑1 (MQ=255)                                                                                                                         gACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGcc                 >  3:186336/1‑148 (MQ=255)                                                                                                                           ggtggtgGATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGTACGCATTCTGGAGGGTGGTGAGCGGCTGGACATCCAGGACGCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTGGACAGCTATGACGCCGAg              <  7:289178/145‑1 (MQ=255)                                                                                                                               gtgCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCgg                                                                      <  7:73416/89‑1 (MQ=255)                                                                                                                               gtgCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCgg                                                                      >  8:73416/1‑89 (MQ=255)                                                                                                                                ttcATCGAGGACGACGGGCCGGGCGTGGCCGTCGACGAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGa          <  2:256882/146‑1 (MQ=255)                                                                                                                                tgCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGa          <  5:61150/148‑1 (MQ=255)                                                                                                                                  cATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGc        <  8:63354/148‑1 (MQ=255)                                                                                                                                   aTCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCt      >  3:216099/1‑149 (MQ=255)                                                                                                                                    tCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCAGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCtt     >  1:72251/1‑149 (MQ=255)                                                                                                                                      gAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCTTgg   <  7:136684/149‑1 (MQ=255)                                                                                                                                       aGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCTTGGc  >  6:40127/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                         GGCGCAGGTACCGATGGAGCAGGGCGCGTTGTTGGAAATGCTCGGCAACCTGCTGGAGAACGCCTACCGGCTGAGCCTGGGCCAAGTGCGGGTGAGCCTGATGAAAGCGCCGGGGTATCTGACCTTGTGCATCGAGGACGACGGGCCGGGCGTGCCCGTCGACCAGCGTGAACGCATTCTGGAGCGTGGGGAGCGGCTGGACAGCCAGCACCCGGGGCAGGGCATCGGCTTGGCAGTGGTCAAGGATATCGTCGACAGCTATGACGCCGAGCTGAGCCTTGGC  >  NC_002947/1362624‑1362906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑