Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,443,931	A→C	6.6%	P299P (CCT→CCG)	PP_4781 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,443,931	0	A	C	6.6%	160.1 / 2.1	62	P299P (CCT→CCG)	PP_4781	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

ACCAGGCCACCAAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGAACAGCCCGGCTGCACCAATGAACA  >  NC_002947/5443782‑5444076                                                                                                                                                      |                                                                                                                                                   aCCAGGCCACCAAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTaa                                                                                                                                                    <  1:956857/150‑1 (MQ=255) aCCAGGCCACCAAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTaa                                                                                                                                                    <  1:9421/150‑1 (MQ=255)        caccaAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCAACCAGGCATACAGGCAGGAGCAGAAGACGCTTGTGGTGAACAGCAGCGCAAGAAGGGTGGCAATAGGCAGGGAGAAGGCGTAAGGCATGc                                                                                                                                             >  1:510195/1‑150 (MQ=255)        cacaaAGCAGGAACACCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGc                                                                                                                                             >  2:714515/1‑150 (MQ=255)         accaAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAg                                                                                                                                                   >  2:1111129/1‑143 (MQ=255)         accaAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAg                                                                                                                                                   <  1:1111129/143‑1 (MQ=255)         accaAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATgtt                                                                                                                                            >  1:122368/1‑148 (MQ=255)           caAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCa                                                                                                                                          >  1:901028/1‑150 (MQ=255)            aaGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCTTACGGCAttgcgct                                                                                                                                         <  2:37522/150‑7 (MQ=255)            aaGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCAc                                                                                                                                         >  1:1046202/1‑150 (MQ=255)                aGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGc                                                                                                                                      <  2:856994/149‑1 (MQ=255)                aGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGc                                                                                                                                      >  1:856994/1‑149 (MQ=255)                 ggAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGAAATGCCGTAAGGCATGCTCACCGCCa                                                                                                                                    <  1:214037/150‑1 (MQ=255)                       ccGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGAAGCGCAACGAGCATGCGCCTGTGGCGAGCAATGGTGTATTGCGTGCTCACCGCCAGTaacc                                                                                                                              >  2:918789/1‑147 (MQ=255)                         gTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGGCAGTACCATg                                                                                                                            >  1:248618/1‑150 (MQ=255)                           ccAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGAc                                                                                                                          <  1:288387/150‑1 (MQ=255)                           ccAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGAc                                                                                                                          <  1:43748/150‑1 (MQ=255)                                 gatgatGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAg‑cgc                                                                                                                   >  1:225587/1‑150 (MQ=255)                                    gatgaACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGc                                                                                                                                                 <  2:1027270/118‑1 (MQ=255)                                    gatgaACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGc                                                                                                                                                 >  1:1027270/1‑118 (MQ=255)                                     atgaACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGccccc                                                                                                               <  2:398387/150‑1 (MQ=255)                                        aaCCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGcccccccc                                                                                                            <  1:416003/150‑1 (MQ=255)                                              cgcgcCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATc                                                                                                      <  2:901028/150‑1 (MQ=255)                                                   cACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCCGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGAGCAG‑CGCCCCCCCCAGCATcagcag                                                                                                 >  1:162435/1‑150 (MQ=255)                                                   cACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGAc                                                                                                                          >  2:659081/1‑126 (MQ=255)                                                   cACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGAc                                                                                                                          <  1:659081/126‑1 (MQ=255)                                                    aCTCGCAAGCCTAGCCACAAAGCATACAGGCCGGCGCTGAAGACGCTTATGGTAAACAGCAGCGCCAGTAGCGGGGCCATGCGCAGCGCAATGGAGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCACGCCCCAGCATCAGAAGt                                                                                                >  2:386870/1‑150 (MQ=255)                                                     cTCGCAAGCCGCGCTACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGACGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGtt                                                                                               <  2:310150/150‑1 (MQ=255)                                                        gCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCTTACGGCAt                                                                                                                                               >  1:115020/1‑100 (MQ=255)                                                        gCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCTTACGGCAt                                                                                                                                               <  2:115020/100‑1 (MQ=255)                                                        gCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCa                                                                                            >  2:1082156/1‑150 (MQ=255)                                                           aGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACa                                                                                         <  2:539325/150‑1 (MQ=255)                                                                 gcCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAAGGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGAt                                                                                   <  1:442370/150‑1 (MQ=255)                                                                  ccaccaGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATc                                                                                  >  2:328863/1‑150 (MQ=255)                                                                     ccaGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCAcc                                                                                                                                        <  2:142285/94‑1 (MQ=255)                                                                     ccaGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCAcc                                                                                                                                        >  1:142285/1‑94 (MQ=255)                                                                      caGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAg‑c                                                                                                                     >  2:1106825/1‑111 (MQ=255)                                                                      caGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAg‑c                                                                                                                     <  1:1106825/111‑1 (MQ=255)                                                                         gCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCATCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCa                                                                           >  1:295631/1‑150 (MQ=255)                                                                                 gCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAg                                                                   <  1:556896/150‑1 (MQ=255)                                                                                    ggCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCt                                                                >  2:496554/1‑150 (MQ=255)                                                                                      cgcTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTGACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGc                                                              <  2:150275/150‑1 (MQ=255)                                                                                       ccTGAACACCCTTATGCTCAACAGCACCGCCACCAGCGCGGCCATGCGCAGCGCAATGCCGTAAGCCACTCTCACCGCCACTACCATGCCCAGCCCCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGc                                                              <  2:510195/149‑1 (MQ=255)                                                                                            aGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAAtt                                                        <  2:1128883/150‑1 (MQ=255)                                                                                             gACGCTTATGGTGAACAGCAGCGCCAGGCGGGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTg                                                       >  1:634916/1‑150 (MQ=255)                                                                                                gCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAg                                                    <  1:787233/150‑1 (MQ=255)                                                                                                gCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAg                                                    <  1:1082156/150‑1 (MQ=255)                                                                                                       gTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATca                                                                                                     <  2:548539/94‑1 (MQ=255)                                                                                                       gTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATca                                                                                                     >  1:548539/1‑94 (MQ=255)                                                                                                       gTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATACCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCa                                                                                                           >  1:599223/1‑88 (MQ=255)                                                                                                       gTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATACCGTAAGGCATGCTCACCGCCAGTACCATGAACAG‑CGCCCCCCCCa                                                                                                           <  2:599223/88‑1 (MQ=255)                                                                                                        tGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGcccccccc                                                                                                            <  2:812493/86‑1 (MQ=255)                                                                                                        tGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGcccccccc                                                                                                            >  1:812493/1‑86 (MQ=255)                                                                                                             agcagcGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCg                                       >  1:792904/1‑150 (MQ=255)                                                                                                                  cgcCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCACCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTg                                  <  1:714515/150‑1 (MQ=255)                                                                                                                   gcCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCACGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGCGGATGCTCAGTTGCAGGAAGTGCGGGCCGAGCTgg                                 >  1:1156327/1‑150 (MQ=255)                                                                                                                           cGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCCCACCGCCAGAACGTAGACCAG‑GTCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGaa                         <  2:634916/150‑1 (MQ=255)                                                                                                                                ccATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCgg                                          >  1:613453/1‑128 (MQ=255)                                                                                                                                ccATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCgg                                          <  2:613453/128‑1 (MQ=255)                                                                                                                                 cATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGaa                         >  2:424853/1‑144 (MQ=255)                                                                                                                                 cATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGaa                         <  1:424853/144‑1 (MQ=255)                                                                                                                                       ggtgagcaTGCCTTACGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGAACAGCCCGGCTGc             >  2:16799/8‑150 (MQ=255)                                                                                                                                               tGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGAACAGCCCGGCTGCACCAATGa     >  1:269868/1‑150 (MQ=255)                                                                                                                                                  cGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCGGGGCGAACTGGCAACCGAACAGCCCGGGCGCACCAAGGAACa  >  2:320382/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ACCAGGCCACCAAGCAGGAACGCCGTCCAGGCGATGATGAACCAGCGCGCCACTCGCAAGCCGCGCCACCAGGCATACAGGCCGGCGCTGAAGACGCTTATGGTGAACAGCAGCGCCAGCAGCGTGGCCATGCGCAGCGCAATGCCGTAAGGCATGCTCACCGCCAGTACCATGACCAG‑CGCCCCCCCCAGCATCAGCAGTTGCAACAGGCGATCGAAGCCACGGTTGAGGCTGCCCAATTGCAGGAAGTGCCGGGCGAACTGGCAACCGAACAGCCCGGCTGCACCAATGAACA  >  NC_002947/5443782‑5444076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑