Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,551,118	A→G	8.0%	Q76R (CAG→CGG) ‡	PP_3135 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,551,118	0	A	G	8.0%	161.1 / 2.4	58	Q76R (CAG→CGG) ‡	PP_3135	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base G (2/3);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGG  >  NC_002947/3550969‑3551265                                                                                                                                                      |                                                                                                                                                    gcgcCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCa                                                                                                                                                     <  2:128555/150‑1 (MQ=255)  cgcCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAg                                                                                                                                                    <  1:987862/150‑1 (MQ=255)   gcCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCGGCCCAGGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAgt                                                                                                                                                   <  2:89334/150‑1 (MQ=255)        gCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGttgttg                                                                                                                                              <  2:1036468/150‑1 (MQ=255)         ccGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGTCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             <  1:1064807/150‑1 (MQ=255)          cGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTGgc                                                                                                                                            <  2:978651/150‑1 (MQ=255)          cGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgc                                                                                                                                            >  1:1085172/1‑150 (MQ=255)                    cgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAg                                                                                                                                  <  2:846046/150‑1 (MQ=255)                       cTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGc                                                                                                                               >  1:582006/1‑150 (MQ=255)                          gCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAg                                                                                                                            >  1:310106/1‑150 (MQ=255)                           ccAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGt                                                                                                                           <  2:582006/150‑1 (MQ=255)                            cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                          <  1:496614/150‑1 (MQ=255)                            cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                          >  1:597553/1‑150 (MQ=255)                            cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                          >  1:692884/1‑150 (MQ=255)                                       tggtCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTtgctgc                                                                                                               <  2:597553/150‑1 (MQ=255)                                        ggtCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGAAAGTATGCCTTgctgct                                                                                                              <  2:244029/150‑1 (MQ=255)                                           cGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCaa                                                                                                           >  1:488913/1‑150 (MQ=255)                                            gACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAAt                                                                                                          >  1:807739/1‑150 (MQ=255)                                              cAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATAc                                                                                                        >  1:498069/1‑150 (MQ=255)                                               aaTGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             >  1:238308/1‑112 (MQ=255)                                               aaTGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             <  2:238308/112‑1 (MQ=255)                                                                cTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAg                                                                                                                                                    <  2:519446/88‑1 (MQ=255)                                                                cTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgc                                                                                      >  1:752553/1‑150 (MQ=255)                                                                cTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgc                                                                                      >  2:599047/1‑150 (MQ=255)                                                                cTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGAGCCAACAACCAg                                                                                                                                                    >  1:519446/1‑88 (MQ=255)                                                                 tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                     >  1:541287/1‑150 (MQ=255)                                                                        aGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGc                                                                              <  2:310106/150‑1 (MQ=255)                                                                        aGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGc                                                                              >  2:774601/1‑150 (MQ=255)                                                                           gcATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGc                                                                           >  1:993383/1‑150 (MQ=255)                                                                           gcATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGc                                                                           <  2:692884/150‑1 (MQ=255)                                                                                 cccAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCaaaaaa                                                                     >  2:600902/1‑150 (MQ=255)                                                                                     ggTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 >  2:961066/1‑150 (MQ=255)                                                                                     ggTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 >  2:915730/1‑150 (MQ=255)                                                                                     gatgggCCCCGAGGCACGGCGGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 <  2:1085172/145‑1 (MQ=255)                                                                                      gtgtGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTc                                                                >  2:323754/1‑150 (MQ=255)                                                                                             ccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACa                                                         >  1:664520/1‑150 (MQ=255)                                                                                              cGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACAt                                                        <  1:774601/150‑1 (MQ=255)                                                                                                        ctgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAc                                              <  1:961066/150‑1 (MQ=255)                                                                                                            tgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCgg                                          >  1:363133/1‑150 (MQ=255)                                                                                                                 aacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCGAAAAACACTCGACTACATGGACGCCCACCCGGAGTgc                                     >  1:396732/1‑150 (MQ=255)                                                                                                                    aaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcggc                                  >  1:15602/1‑150 (MQ=255)                                                                                                                     aGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcggcg                                 >  1:1085472/1‑150 (MQ=255)                                                                                                                           aCGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTgg                           <  2:664520/150‑1 (MQ=255)                                                                                                                               gggTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAACAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGt                       <  1:604386/150‑1 (MQ=255)                                                                                                                                     cGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCt                 <  2:498069/150‑1 (MQ=255)                                                                                                                                         cGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             <  1:599047/150‑1 (MQ=255)                                                                                                                                         cGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             <  1:915730/150‑1 (MQ=255)                                                                                                                                            gCCAACAACTGGTCGTTGGCGCACCTCAAGGGAAAGGATGCCTTGCTGCGCAATACCGAGGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGt          >  2:1030968/1‑150 (MQ=255)                                                                                                                                             taaacaacCAGTGGTTGGCGGACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTc         <  2:1085472/148‑1 (MQ=255)                                                                                                                                              caacaaCCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcg        >  1:807150/1‑150 (MQ=255)                                                                                                                                               aacaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcgc       >  2:436432/1‑150 (MQ=255)                                                                                                                                                acaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcgcg      >  1:1088577/1‑150 (MQ=255)                                                                                                                                                  aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt    >  2:1106887/1‑150 (MQ=255)                                                                                                                                                  aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt    <  1:323754/150‑1 (MQ=255)                                                                                                                                                  aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt    >  2:1129561/1‑150 (MQ=255)                                                                                                                                                  aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGACAACTGCGCACGACTACATGGGCGCCCACCAGGACTGCGGCGTCCTGGGGGTTCGCCTGGAGGGTCGCGAt    >  1:860878/1‑150 (MQ=255)                                                                                                                                                   acCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATg   <  1:1106887/150‑1 (MQ=255)                                                                                                                                                    ccAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATgg  <  2:921000/150‑1 (MQ=255)                                                                                                                                                    ccAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATgg  >  2:945666/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGG  >  NC_002947/3550969‑3551265

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