Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,948,216	C→G	7.2%	R494R (CGG→CGC)	PP_1749 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,948,216	0	C	G	7.2%	172.0 / 4.2	58	R494R (CGG→CGC)	PP_1749	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGA  >  NC_002947/1948073‑1948358                                                                                                                                                |                                                                                                                                               ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                                                                          >  2:115023/1‑150 (MQ=255) ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                                                                          <  1:7520/150‑1 (MQ=255)    gCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGcc                                                                                                                                       <  1:1065644/150‑1 (MQ=255)      aTCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAg                                                                                                                                     <  1:674009/150‑1 (MQ=255)        cACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGCTGGCCGTGCGGCGCACGGCCaacc                                                                                                                                   <  1:1572536/150‑4 (MQ=255)          cGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGgcgc                                                                                                                                 <  2:864711/150‑1 (MQ=255)             aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                                                                                                              >  1:417808/1‑150 (MQ=255)             aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGAACGGCCAGGCGCTgt                                                                                                                              <  2:1142995/150‑2 (MQ=255)               cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                                                                          >  2:270032/1‑136 (MQ=255)               cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                                                                          <  1:270032/136‑1 (MQ=255)               cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGGGTGCAGGTCCGCCGTGCGCCGCACGGCCAGGCGCTggcc                                                                                                                            >  1:660533/1‑150 (MQ=255)                     gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGccggccgg                                                                                                                      <  1:1071775/150‑1 (MQ=255)                     gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGccggccgg                                                                                                                      <  1:1540970/150‑1 (MQ=255)                                   ccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCAt                                                                                                        >  2:1579510/1‑150 (MQ=255)                                   ccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCAt                                                                                                        >  1:1022735/1‑150 (MQ=255)                                    cGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATa                                                                                                       <  2:417808/150‑1 (MQ=255)                                           tCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACg                                                                                                >  2:318923/1‑150 (MQ=255)                                           tCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACg                                                                                                >  1:1281196/1‑150 (MQ=255)                                            ccAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTAACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGc                                                                                               >  1:698461/1‑150 (MQ=255)                                                agCCCGTGCCGCGCGCACGGCAGGAGGGGCCAGCCCCGGGTGCGGGCGTTGCGTGCCGTCCTCAAGGGTGCTGACGGCATGGAGGGTGTTAGCGCGTTGGAGGGCGAGGCGTTTGCCGGGGGGGTTGACGTCGTGATTGCCAACGcccac                                                                                           <  2:1421304/149‑3 (MQ=255)                                                            gcgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGt                                                                               <  1:1244133/150‑1 (MQ=255)                                                                aCGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTgcgc                                                                           <  2:1316872/150‑1 (MQ=255)                                                                  ggCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCa                                                                                                         <  2:751991/118‑1 (MQ=255)                                                                  ggCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCa                                                                                                         >  1:751991/1‑118 (MQ=255)                                                                    cagcaTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCgg                                                                       >  1:78778/1‑150 (MQ=255)                                                                           ggCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc                                                                >  2:362738/1‑150 (MQ=255)                                                                                   ccGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGc                                                                                                                                        >  2:1518196/1‑70 (MQ=255)                                                                                   ccGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGc                                                                                                                                        <  1:1518196/70‑1 (MQ=255)                                                                                            ggCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGc                                               >  2:1066930/1‑150 (MQ=255)                                                                                               gTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCgg                                                >  2:131761/1‑146 (MQ=255)                                                                                               gTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCgg                                                <  1:131761/146‑1 (MQ=255)                                                                                               gTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGctct                                            >  2:186609/1‑150 (MQ=255)                                                                                                      gACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCAcc                                     >  1:1074698/1‑150 (MQ=255)                                                                                                         gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCg                                                 >  1:1354939/1‑135 (MQ=255)                                                                                                         gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCg                                                 <  2:1354939/135‑1 (MQ=255)                                                                                                         gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGt                                  <  2:1037317/150‑1 (MQ=255)                                                                                                          tcctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTg                                 <  1:1448094/150‑1 (MQ=255)                                                                                                           cctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTgg                                <  1:20016/150‑1 (MQ=255)                                                                                                            ctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATa                                                                                                       >  1:1521516/1‑78 (MQ=255)                                                                                                            ctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATa                                                                                                       <  2:1521516/78‑1 (MQ=255)                                                                                                            ctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTGGGTTTCGGCGTCCAGCGGAATGCGGCTCTCGCGACCGGTGGc                               >  1:1306577/1‑150 (MQ=255)                                                                                                             tccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCa                              >  2:1631845/1‑150 (MQ=255)                                                                                                               cAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCg                                                 >  2:1637998/1‑129 (MQ=255)                                                                                                               cAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCg                                                 <  1:1637998/129‑1 (MQ=255)                                                                                                               cAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGc                            <  1:362738/150‑1 (MQ=255)                                                                                                                 gggTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCt                          <  2:611040/150‑1 (MQ=255)                                                                                                                    tGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       <  2:1074698/150‑1 (MQ=255)                                                                                                                      ccgccgGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGt                     <  1:1579510/150‑1 (MQ=255)                                                                                                                       cgccgGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTc                    >  1:1204225/1‑150 (MQ=255)                                                                                                                          cgGTGTGCCGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGAAACACGTCGTCATAGCCAACGCCCGCTGCGAGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGAc                 <  1:115023/150‑1 (MQ=255)                                                                                                                          cgGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGAc                 >  2:20300/1‑150 (MQ=255)                                                                                                                           ggTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACg                <  2:1022735/150‑1 (MQ=255)                                                                                                                                 cAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       >  2:1128749/1‑137 (MQ=255)                                                                                                                                 cAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       <  1:1128749/137‑1 (MQ=255)                                                                                                                                     ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       <  2:1113472/132‑1 (MQ=255)                                                                                                                                     ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       >  1:1113472/2‑133 (MQ=255)                                                                                                                                     ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTAGCGTCGACGGCTCTGCGCt      >  1:640169/2‑150 (MQ=255)                                                                                                                                         gccagcgcCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCAc                                      >  2:1404808/5‑114 (MQ=255)                                                                                                                                         gccagcgcCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCAc                                      <  1:1404808/110‑1 (MQ=255)                                                                                                                                         ccGGGCGCGGCACGGCCAGGCGCTGGCCGGCCGGCAACACGGCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGa  <  2:660533/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               CCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGA  >  NC_002947/1948073‑1948358

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