Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,948,173	C→G	9.5%	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,173	0	C	G	9.5%	188.7 / 10.3	63	R494R (CGG→CGC)	PP_1749	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (4/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

GCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGAT  >  NC_002947/1948025‑1948316                                                                                                                                                     |                                                                                                                                                 gCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgcgc                                                                                                                                                  >  4:27477/1‑149 (MQ=255)  cTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgcc                                                                                                                                                 <  1:176547/149‑1 (MQ=255)     cccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCa                                                                                                                                              >  3:285672/1‑149 (MQ=255)      ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                                                                                                             >  5:262816/1‑149 (MQ=255)      ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                                                                                                             <  6:262816/149‑1 (MQ=255)      ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                                                                                                             >  4:236689/1‑149 (MQ=255)         gCATCACGCACCATGAAGTCCAGCCCCACCACCGGACTTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGGGCCGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGc                                                                                                                                          <  7:307739/149‑1 (MQ=255)                     aTGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTggcc                                                                                                                              <  2:133109/149‑1 (MQ=255)                        aaGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTggccggc                                                                                                                           <  8:220984/149‑1 (MQ=255)                         aGTCCAGCCCCTCCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGGAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTggccggcc                                                                                                                          <  4:227653/149‑1 (MQ=255)                         aGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTggccggcc                                                                                                                          <  3:236689/149‑1 (MQ=255)                              aGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                                                                                     <  8:224139/149‑1 (MQ=255)                                 cccACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAac‑ac                                                                                                                 >  8:203545/1‑149 (MQ=255)                                   caccacCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCACC‑Acgt                                                                                                               >  6:92439/1‑149 (MQ=255)                                    accaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑Acgtc                                                                                                              <  4:285672/149‑1 (MQ=255)                                    accaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑Acgtc                                                                                                              <  6:228636/149‑1 (MQ=255)                                       accGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑Acgtcgt                                                                                                            >  3:118728/1‑148 (MQ=255)                                          ggAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCa                                                                                                          >  5:245071/1‑147 (MQ=255)                                          ggAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATa                                                                                                        <  1:207665/149‑1 (MQ=255)                                           gAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAg                                                                                                       >  8:322726/1‑149 (MQ=255)                                           gAATTTCAATTGCCGGTGCCGCGCGCACGGCATCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAg                                                                                                       <  6:252577/149‑1 (MQ=255)                                             aTTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGcc                                                                                                     <  7:322726/149‑1 (MQ=255)                                                  cAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGc                                                                                                >  8:3006/1‑149 (MQ=255)                                                         cGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTgc                                                                                         >  1:131279/1‑149 (MQ=255)                                                          gTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTgcg                                                                                        >  4:79961/1‑149 (MQ=255)                                                            gCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTgcgcg                                                                                      >  6:284384/1‑149 (MQ=255)                                                              cgcgcgCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCCCCCCCCGGCCACGCCGGGGCCCGCCGGGAAA‑CCGGGGGCATACCCAACCCCCCCTgccccc                                                                                     >  1:110914/1‑144 (MQ=255)                                                                  cgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGt                                                                                <  6:245071/149‑1 (MQ=255)                                                                     aCGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACTACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTgc                                                                              <  1:65411/149‑1 (MQ=255)                                                                            caTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGtt                                                                      <  5:284384/149‑1 (MQ=255)                                                                             aTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGttt                                                                     >  6:74842/1‑149 (MQ=255)                                                                                   cAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTgc                                                                              <  5:361970/134‑1 (MQ=255)                                                                                   cAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTgc                                                                              >  6:361970/1‑134 (MQ=255)                                                                                   cAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTcgtcgt                                                               >  7:149165/1‑149 (MQ=255)                                                                                   cAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTcgtcgt                                                               >  4:221658/1‑149 (MQ=255)                                                                                    aGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCATGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCgtcgtc                                                              >  4:339927/1‑149 (MQ=255)                                                                                         cGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCg                                                         >  8:293420/1‑149 (MQ=255)                                                                                                         gTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtc                                                                    <  8:329834/122‑1 (MQ=255)                                                                                                         gTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtc                                                                    >  7:329834/1‑122 (MQ=255)                                                                                                          tGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCATGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCa                                        <  3:339927/149‑1 (MQ=255)                                                                                                          tGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCa                                        >  5:279232/1‑149 (MQ=255)                                                                                                              gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGct                                               <  2:156001/138‑1 (MQ=255)                                                                                                              gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGct                                               >  1:156001/1‑138 (MQ=255)                                                                                                               tcctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGt                                   <  4:118728/149‑1 (MQ=255)                                                                                                                 ctccAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTgg                                 >  8:143171/1‑149 (MQ=255)                                                                                                                    cAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAg                              >  1:173146/1‑149 (MQ=255)                                                                                                                        gTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTg                          >  6:145715/1‑149 (MQ=255)                                                                                                                        gTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTg                          >  8:31608/1‑149 (MQ=255)                                                                                                                          gccgccGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                        <  3:192093/149‑1 (MQ=255)                                                                                                                                    tgCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGct              <  7:143171/149‑1 (MQ=255)                                                                                                                                      cAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGctct            <  7:3006/149‑1 (MQ=255)                                                                                                                                        ggTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGTCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTGCATCGGTATGCTGCTCTCGACACCGGTGGCAGCCTGGCGTCGACGGATCTgc          <  2:118456/149‑1 (MQ=255)                                                                                                                                        ggTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTg                          >  3:6813/1‑133 (MQ=255)                                                                                                                                        ggTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTg                          <  4:6813/133‑1 (MQ=255)                                                                                                                                         cctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCt                           <  3:97907/129‑1 (MQ=255)                                                                                                                                         cctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCt                           >  4:97907/3‑131 (MQ=255)                                                                                                                                         cctGGCCGTGCGGAGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcg         >  6:90970/3‑149 (MQ=255)                                                                                                                                          ttGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGGGTCGACGGCTCTgcgc        <  6:279232/149‑1 (MQ=255)                                                                                                                                          ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGc                                                   >  5:311386/2‑106 (MQ=255)                                                                                                                                          ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGc                                                   <  6:311386/105‑1 (MQ=255)                                                                                                                                          ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcgc        >  3:105887/2‑149 (MQ=255)                                                                                                                                             gCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTc     <  8:149165/149‑1 (MQ=255)                                                                                                                                              ccGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCg    <  7:293420/149‑1 (MQ=255)                                                                                                                                               cGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGa   <  3:221658/149‑1 (MQ=255)                                                                                                                                                gTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCCCCGGTGGCGGCCTGGCGTCGACGGCTCTGCGCTTCGAt  >  8:225824/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAC‑ACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGAT  >  NC_002947/1948025‑1948316

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