Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	318,970	T→C	10.3%	T258A (ACC→GCC)	PP_0264 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	318,970	0	T	C	10.3%	153.8 / 7.3	59	T258A (ACC→GCC)	PP_0264	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base C (2/4);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGA  >  NC_002947/318823‑319116                                                                                                                                                      |                                                                                                                                                   gccgccGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTc                                                                                                                                                   >  5:219236/1‑149 (MQ=255)     ccgTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCAC‑‑ACCAGTTTCCTTGATGTCGCGGGTCTggg                                                                                                                                               >  3:182932/1‑149 (MQ=255)        ttGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGa                                                                                                                                            <  5:157351/149‑1 (MQ=255)        ttGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTggatca                                                                                                                                            >  7:144267/1‑145 (MQ=255)          gCATCACCAGGCCTGCCAGTAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                          >  4:86016/1‑149 (MQ=255)           cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                          >  7:242247/1‑148 (MQ=255)           cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                          >  4:242293/1‑148 (MQ=255)           cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGAATTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                          >  2:91004/1‑148 (MQ=255)            aTCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTggatcacc                                                                                                                                          >  2:175769/1‑141 (MQ=255)              cACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAgg                                                                                                                                      >  7:251060/1‑149 (MQ=255)              cACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAgg                                                                                                                                      >  3:56169/1‑149 (MQ=255)                ccAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGcc                                                                                                                                    <  8:251060/149‑1 (MQ=255)                  aGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGAcccg                                                                                                                                  <  1:175769/149‑3 (MQ=255)                  aGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGAcccg                                                                                                                                  <  8:144267/149‑3 (MQ=255)                   ggCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGc                                                                                                                                 <  4:56169/149‑1 (MQ=255)                   ggCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                                 <  7:256840/149‑1 (MQ=255)                    gCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCg                                                                                                                                >  8:213613/1‑149 (MQ=255)                      cTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTgg                                                                                                                                                >  3:37448/1‑131 (MQ=255)                      cTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTgg                                                                                                                                                <  4:37448/131‑1 (MQ=255)                         ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTg                                                                                                                           <  1:93101/149‑1 (MQ=255)                         ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCcgcgacat                                                                                                                            <  7:30462/148‑5 (MQ=255)                          cAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGTACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTga                                                                                                                          <  2:311479/149‑1 (MQ=255)                          cAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTga                                                                                                                          <  8:242247/149‑1 (MQ=255)                                    aGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCg                                                                                                                 <  3:242293/148‑1 (MQ=255)                                     gTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGATTTGCCACCACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCg                                                                                                                 <  1:91004/149‑1 (MQ=255)                                       aCCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtg                                                                                                             >  4:88804/1‑149 (MQ=255)                                        ccAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgt                                                                                                            <  8:244316/149‑1 (MQ=255)                                          agcagcGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCgg                                                                                                                               <  2:27079/128‑1 (MQ=255)                                          agcagcGTACCGGCACCGATGGCCATTACCCTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCgg                                                                                                                               >  1:27079/1‑128 (MQ=255)                                           gcagcGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTggg                                                                                                         >  6:63024/1‑149 (MQ=255)                                           gcagcGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTggg                                                                                                         <  3:305826/149‑1 (MQ=255)                                                      gCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCg                                                                                               >  3:309502/1‑148 (MQ=255)                                                        aCCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGAt                                                                                            <  3:86016/149‑1 (MQ=255)                                                          cGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTc                                                                                                >  2:89215/1‑143 (MQ=255)                                                          cGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTc                                                                                                <  1:89215/143‑1 (MQ=255)                                                               gCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTg                                                                                     <  8:146743/149‑1 (MQ=255)                                                                           tGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCggg                                                                                                                              >  7:58240/1‑96 (MQ=255)                                                                           tGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCggg                                                                                                                              <  8:58240/96‑1 (MQ=255)                                                                             cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATgg                                                                       >  8:136918/1‑149 (MQ=255)                                                                                 gAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTg                                                                   <  5:63024/149‑1 (MQ=255)                                                                                  aGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTgg                                                                  >  2:161553/1‑149 (MQ=255)                                                                                            cAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGt                                                        >  4:246265/1‑149 (MQ=255)                                                                                            cAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGt                                                        <  7:136918/149‑1 (MQ=255)                                                                                                                 gCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                   >  2:209570/1‑149 (MQ=255)                                                                                                                   tGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACg                                  >  3:126553/1‑148 (MQ=255)                                                                                                                       ttGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTGTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgc                             >  1:25328/1‑149 (MQ=255)                                                                                                                          ccA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                           <  1:209570/148‑1 (MQ=255)                                                                                                                          ccA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                           <  2:25328/148‑1 (MQ=255)                                                                                                                                ccAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGaa                               <  6:312936/140‑1 (MQ=255)                                                                                                                                ccAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGaa                               >  5:312936/1‑140 (MQ=255)                                                                                                                                       ccTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                            >  4:256341/1‑88 (MQ=255)                                                                                                                                       ccTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                            <  3:256341/88‑1 (MQ=255)                                                                                                                                        cTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCcgc               >  6:325735/1‑148 (MQ=255)                                                                                                                                         ttGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCcgcg              >  6:105319/1‑148 (MQ=255)                                                                                                                                           taTGTCGCGGGGCGGGGTGATCCAGGCCGGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGaa            <  4:182932/147‑1 (MQ=255)                                                                                                                                            aTGGCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGt          <  7:213613/149‑1 (MQ=255)                                                                                                                                                 gcgGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAgg      <  3:88804/148‑1 (MQ=255)                                                                                                                                                  cgGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGAt    <  4:126553/149‑1 (MQ=255)                                                                                                                                                    ggTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGa  <  6:219236/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCA‑‑ACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGA  >  NC_002947/318823‑319116

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