Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,453,178	A→G	6.6%	H447R (CAT→CGT)	PP_1271 →	putative Multidrug efflux MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,453,178	0	A	G	6.6%	175.0 / 2.4	60	H447R (CAT→CGT)	PP_1271	putative Multidrug efflux MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (35/21);  new base G (2/2);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CTGCCGTTGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGCCCTGAT  >  NC_002947/1453033‑1453325                                                                                                                                                  |                                                                                                                                                     ctgcCGTTGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATttccattc                                                                                                                                                   <  8:40892/149‑1 (MQ=255) ctgcCGTTGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCAGGATGTTCTTGGCACCCTGCGCCGGCATttccattc                                                                                                                                                   <  1:241176/149‑1 (MQ=255)     cGTTGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAAt                                                                                                                                               >  8:180258/1‑149 (MQ=255)       ttGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCa                                                                                                                                             >  3:116528/1‑149 (MQ=255)            gATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCGTTCCAATCACCtgg                                                                                                                                        >  3:65507/1‑149 (MQ=255)            gATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCGTTCCAATCACCtgg                                                                                                                                        >  3:63568/1‑149 (MQ=255)                     caccaATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGcc                                                                                                                               >  5:117953/1‑149 (MQ=255)                      accaATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCt                                                                                                                              >  7:279973/1‑149 (MQ=255)                                  aCCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGccc                                                                                                                   <  6:88266/148‑1 (MQ=255)                                   ccccGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCg                                                                                                                 >  1:57040/1‑149 (MQ=255)                                       gCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTgg                                                                                                             <  4:116528/149‑1 (MQ=255)                                           gaggGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTgg                                                                                                                            >  2:68912/1‑130 (MQ=255)                                           gaggGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTgg                                                                                                                            <  1:68912/130‑1 (MQ=255)                                            aggGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCACGCTGGTCGAc                                                                                                        <  8:279973/149‑1 (MQ=255)                                             gggCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACg                                                                                                       <  8:99526/149‑1 (MQ=255)                                                cccGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGc                                                                                                     >  4:20997/1‑148 (MQ=255)                                                 ccGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCaa                                                                                                   >  6:333982/1‑149 (MQ=255)                                                  cGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAAc                                                                                                  <  2:57040/149‑1 (MQ=255)                                                    ttCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAAcgc                                                                                                <  3:20997/149‑1 (MQ=255)                                                          tCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGctgg                                                                                          >  4:49219/1‑149 (MQ=255)                                                               cTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGctggcctgg                                                                                     >  5:316595/1‑149 (MQ=255)                                                                  gTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCGTTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGctggcctggc                                                                                    <  4:63568/147‑1 (MQ=255)                                                                  gTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCGTTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGctggcctggc                                                                                    <  4:65507/147‑1 (MQ=255)                                                                   ttCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCtggcctggcc                                                                                   >  2:281219/1‑147 (MQ=255)                                                                     cAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑gc                                                                              >  2:196441/1‑149 (MQ=255)                                                                      aaCACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTAGACGGCAACGCCGCTGGCCTGGCCAA‑gcg                                                                             >  3:234262/1‑149 (MQ=255)                                                                       acacCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGACACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑gcg                                                                             >  1:92809/1‑148 (MQ=255)                                                                               aaGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTc                                                                        >  4:267984/1‑145 (MQ=255)                                                                               aaGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTc                                                                        <  3:267984/145‑1 (MQ=255)                                                                                        gccgccGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCa                                                           <  7:180258/149‑1 (MQ=255)                                                                                         ccgccgTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAg                                                          <  6:117953/149‑1 (MQ=255)                                                                                              gTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTAtt                                                     >  1:262817/1‑149 (MQ=255)                                                                                                cgatgacggcggCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGa                                                   >  7:101234/7‑149 (MQ=255)                                                                                                   cGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGAGAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCAtt                                                                         >  3:134792/1‑124 (MQ=255)                                                                                                   cGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCAtt                                                                         <  4:134792/124‑1 (MQ=255)                                                                                                     ccGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCCCGACCCGCCGCAGGTATTGACCTCg                                              >  6:194397/1‑149 (MQ=255)                                                                                                        gcggcCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTg                                                    <  5:48508/140‑1 (MQ=255)                                                                                                        gcggcCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTg                                                    >  6:48508/1‑140 (MQ=255)                                                                                                             ctgctgGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACct                                      <  2:92809/149‑1 (MQ=255)                                                                                                             ctgctgGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACct                                      >  6:101043/1‑149 (MQ=255)                                                                                                                   gATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTAtt                                                     >  1:342695/1‑128 (MQ=255)                                                                                                                   gATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTAtt                                                     <  2:342695/128‑1 (MQ=255)                                                                                                                   gATGTTCTTGGCACCCTGCGCCGGCATTTCCATACCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGCCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCACTATCt                                >  2:277985/1‑149 (MQ=255)                                                                                                                      gTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGt                             >  1:166565/1‑149 (MQ=255)                                                                                                                      gTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGt                             >  1:197512/1‑149 (MQ=255)                                                                                                                       ttCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTggct                            >  3:71991/1‑147 (MQ=255)                                                                                                                        tCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAg                                                                <  1:162314/112‑1 (MQ=255)                                                                                                                        tCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAg                                                                >  2:162314/1‑112 (MQ=255)                                                                                                                        tCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGc                           >  3:352900/1‑149 (MQ=255)                                                                                                                          ttGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTg                         >  4:136653/1‑149 (MQ=255)                                                                                                                            ggCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGAc                                                  >  2:37920/1‑122 (MQ=255)                                                                                                                            ggCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGAc                                                  <  1:37920/122‑1 (MQ=255)                                                                                                                            ggCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTgg                        <  5:101043/148‑1 (MQ=255)                                                                                                                            ggCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTgg                        <  4:352900/148‑1 (MQ=255)                                                                                                                            ggCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGc                       >  2:157581/1‑149 (MQ=255)                                                                                                                             gCACCCGGCGCCGGCTTTTCCATTCCAAACACTGGGTGCCCAGCTTGGGAAACCCCCGGCTGGTGTAGGGCAACCCGGGGGGCCGGGCAAACGCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGc                       <  1:281219/149‑1 (MQ=255)                                                                                                                              cACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCt                     <  2:256800/149‑1 (MQ=255)                                                                                                                                        cGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGt           >  1:349330/1‑149 (MQ=255)                                                                                                                                          gCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCa                                                           >  6:157771/1‑99 (MQ=255)                                                                                                                                          gCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCa                                                           <  5:157771/99‑1 (MQ=255)                                                                                                                                            aTTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCGTGTATCGCCGTGGcc       >  2:244061/1‑149 (MQ=255)                                                                                                                                            aTTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGcc       >  3:60551/1‑149 (MQ=255)                                                                                                                                              ttccattccaATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGCCCt     >  8:43170/1‑149 (MQ=255)                                                                                                                                                 cattccaATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGCCCTGa   <  3:136653/148‑1 (MQ=255)                                                                                                                                                 cattccaATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGCCCTGAt  >  1:286467/1‑149 (MQ=255)                                                                                                                                                  attccaATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACct                                      <  3:52078/112‑1 (MQ=255)                                                                                                                                                  attccaATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACct                                      >  4:52078/1‑112 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTGCCGTTGCTGATGCTGTCCACCAATGGCATGACCCCGCAGGAGGGCCCGTTCGCCTCCAGCTGGTTCAACACCGTCAAGGGCCTGGCCGCCGTCATCGCCGGCGGCCTGCTGGATGTTCTTGGCACCCTGCGCCGGCATTTCCATTCCAATCACCTGGTGGACAGCCTGGGCAACGCCCCGCTGGTCGACGGCAACGCCGCTGGCCTGGCCAA‑GCGCATTCACGACCAGGCGCAGGTATTGACCTCGGCTGACCTCTATCTGGTGCTGGCCTGTATCGCCGTGGCCCTGAT  >  NC_002947/1453033‑1453325

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