Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,695,065	G→A	5.6%	R330Q (CGG→CAG)	wspC →	biofilm formation methyltransferase WspC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,695,065	0	G	A	5.6%	212.1 / 2.8	71	R330Q (CGG→CAG)	wspC	biofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base A (2/2);  total (33/38)
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.19e-01

TTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCG  >  NC_002947/1694920‑1695212                                                                                                                                                  |                                                                                                                                                    ttCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTg                                                                                                                                                 >  2:509607/1‑150 (MQ=255)   cACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGc                                                                                                                                               <  1:638204/150‑1 (MQ=255)        cgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATg                                                                                                                                          <  1:509607/150‑1 (MQ=255)         gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGc                                                                                                                                         >  1:250764/1‑150 (MQ=255)         atctcccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTCGCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGc                                                                                                                                         <  1:15460/146‑1 (MQ=255)           tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                       >  2:1104304/1‑150 (MQ=255)           tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                       <  1:1284605/150‑1 (MQ=255)              ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGCGGGCGCCACCGCGTTCCTGTCCTGCATTTCCCGGCGGGCCAATGccagct                                                                                                                                    >  2:15460/1‑146 (MQ=255)              ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                    >  1:1547842/1‑150 (MQ=255)                    cGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCg                                                                                                                              >  2:1234513/1‑150 (MQ=255)                             cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCTGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                     >  2:1094987/1‑150 (MQ=255)                              gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                    >  2:122874/1‑150 (MQ=255)                                aGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGACAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                  <  1:323588/150‑1 (MQ=255)                                 gCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAg                                                                                                                 <  2:174579/150‑1 (MQ=255)                                     aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                          <  2:1288240/137‑1 (MQ=255)                                     aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                          >  1:1288240/1‑137 (MQ=255)                                     aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                             <  1:3410/150‑1 (MQ=255)                                     aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATCGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                             >  2:301269/1‑150 (MQ=255)                                      gtgtGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCTGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                            <  1:1094987/150‑1 (MQ=255)                                      gtgtGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                            <  1:164684/150‑1 (MQ=255)                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                   <  1:67652/116‑1 (MQ=255)                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                   >  2:67652/1‑116 (MQ=255)                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 <  1:122874/150‑1 (MQ=255)                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 <  1:1234513/150‑1 (MQ=255)                                                         gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                         <  2:250764/150‑1 (MQ=255)                                                         gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                         <  2:767570/150‑1 (MQ=255)                                                                   ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                      <  2:176532/127‑1 (MQ=255)                                                                   ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                      >  1:176532/1‑127 (MQ=255)                                                                    tgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGt                                                                              <  2:1547842/150‑1 (MQ=255)                                                                      ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCg                                                                            >  2:1093657/1‑150 (MQ=255)                                                                        ggTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCAAGTCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGc                                                                          >  1:198389/1‑150 (MQ=255)                                                                            gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAg                                                                      <  2:494246/150‑1 (MQ=255)                                                                            gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCAACGAGGAGTCTAGTCTTCGACTTGCCGGCATTGCCCGGCTGTTATATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAg                                                                      <  1:774253/150‑1 (MQ=255)                                                                              ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCGCAGgt                                                                    <  1:1093657/150‑1 (MQ=255)                                                                                 caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                             >  2:1232756/1‑106 (MQ=255)                                                                                 caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                             <  1:1232756/106‑1 (MQ=255)                                                                                    gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                      <  2:1630045/110‑1 (MQ=255)                                                                                    gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                      >  1:1630045/1‑110 (MQ=255)                                                                                    ccccgcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                              >  1:939774/5‑150 (MQ=255)                                                                                     ccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGc                                                                                                                                                   <  1:672433/64‑1 (MQ=255)                                                                                     ccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGc                                                                                                                                                   >  2:672433/1‑64 (MQ=255)                                                                                      cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACTCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTactact                                                            >  2:993795/1‑150 (MQ=255)                                                                                       agCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActg                                                           >  2:25743/1‑150 (MQ=255)                                                                                          gcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggct                                                        <  2:1554170/150‑1 (MQ=255)                                                                                           cgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctg                                                       <  1:25743/150‑1 (MQ=255)                                                                                               cccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                                   <  1:626737/150‑1 (MQ=255)                                                                                                ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                     <  2:1276045/83‑1 (MQ=255)                                                                                                ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCGGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                     >  1:1276045/1‑83 (MQ=255)                                                                                                      gCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTAAGCTGATACCGTCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                            <  1:1104304/150‑1 (MQ=255)                                                                                                       cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                           <  2:939774/150‑1 (MQ=255)                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGt                                                                              <  1:1026972/110‑1 (MQ=255)                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGt                                                                              >  2:1026972/1‑110 (MQ=255)                                                                                                              aCGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCg                                    >  2:1051445/1‑150 (MQ=255)                                                                                                              aCGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCg                                    >  1:1079037/1‑150 (MQ=255)                                                                                                                  gggggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                                <  2:1079037/150‑1 (MQ=255)                                                                                                                   ggggCCAACGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAg                                                                                                                 >  1:1524929/1‑68 (MQ=255)                                                                                                                   ggggCCAACGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAg                                                                                                                 <  2:1524929/68‑1 (MQ=255)                                                                                                                     ggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                             >  2:662496/1‑150 (MQ=255)                                                                                                                     ggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                             <  1:662496/150‑1 (MQ=255)                                                                                                                        cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                 >  1:401832/1‑127 (MQ=255)                                                                                                                        cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                 <  2:401832/127‑1 (MQ=255)                                                                                                                                   tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagca                    <  1:305370/145‑1 (MQ=255)                                                                                                                                   tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagca                    >  2:305370/1‑145 (MQ=255)                                                                                                                                      ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                            >  1:519392/1‑118 (MQ=255)                                                                                                                                      ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                            <  2:519392/118‑1 (MQ=255)                                                                                                                                       cGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 >  2:197321/1‑64 (MQ=21)                                                                                                                                       cGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 <  1:197321/64‑1 (MQ=21)                                                                                                                                       cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg           >  1:1093359/1‑150 (MQ=255)                                                                                                                                         gCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                             <  1:1115171/114‑1 (MQ=255)                                                                                                                                         gCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                             >  2:1115171/1‑114 (MQ=255)                                                                                                                                                ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg  >  2:1400592/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCG  >  NC_002947/1694920‑1695212

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