Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,831	C→G	6.1%	intergenic (‑70/+61)	dctB ← / ← bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,598,831	0	C	G	6.1%	197.0 / 4.2	66	intergenic (‑70/+61)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base G (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

TGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTG  >  NC_002947/1598692‑1598978                                                                                                                                             |                                                                                                                                                    tGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCAcc                                                                                                                                            <  3:392501/149‑1 (MQ=255)       cGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAg                                                                                                                                       <  6:321329/148‑1 (MQ=255)       cGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAg                                                                                                                                       <  3:407310/148‑1 (MQ=255)            cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGgcgag                                                                                                                                  >  5:344776/1‑146 (MQ=255)             ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGt                                                                                                                                >  7:317750/1‑149 (MQ=255)              cAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAATGAAAAGGGCATCGGGGGGGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGtt                                                                                                                               >  1:292417/1‑149 (MQ=255)                gCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGtttt                                                                                                                             <  8:10231/149‑1 (MQ=255)                  aCGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGC‑TGCAGAGTACCTGTAGGCGCGGGTTTAc                                                                                                                           <  2:21970/149‑1 (MQ=25)                      tCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGGCCCTGTAGGCGCGGTTTTACCcgc                                                                                                                       >  5:241552/1‑149 (MQ=255)                         gcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGa                                                                                                                     >  2:325544/1‑148 (MQ=255)                          cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                    >  4:78272/1‑148 (MQ=255)                          cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                    >  1:330151/1‑148 (MQ=255)                          cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                    >  2:119511/1‑148 (MQ=255)                           agGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                    >  2:328279/1‑147 (MQ=255)                            gggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAg                                                                                                                 <  2:292417/149‑1 (MQ=255)                            gggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAg                                                                                                                 <  2:330151/149‑1 (MQ=255)                             ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAgg                                                                                                                <  4:133393/149‑1 (MQ=255)                                  tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAg                                                                                                                                       >  7:129116/1‑121 (MQ=255)                                  tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAg                                                                                                                                       <  8:129116/121‑1 (MQ=255)                                    aGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAAc                                                                                                         <  8:317750/149‑1 (MQ=255)                                           gACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcggc                                                                                                  <  8:28441/149‑1 (MQ=255)                                                  ggCGCGAAAGGAAAAGGGCATCGGGGGCGGAGCCGGTGGTGGGATGTCGGGGCCTGCCTCTACGCGGGTAAACCCGCGCCTACAGGTCC‑GCGGCAGCCCCGGTCGCCGCGGGGTTACCCGCGAAAAGGCCGGACCAGGCGGCCttttt                                                                                            >  5:404614/1‑148 (MQ=255)                                                   gcgcGTAAGGAAATGGGCCTGGGGGGCTGATCCGGTGGTGGGAGGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                          <  6:404614/149‑1 (MQ=37)                                                   gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                          <  5:176831/149‑1 (MQ=255)                                                   gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                          <  6:364149/149‑1 (MQ=255)                                                     gcgAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                                        >  1:228735/1‑149 (MQ=255)                                                       gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTgc                                                                                                                                                  >  8:321021/1‑89 (MQ=255)                                                       gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTgc                                                                                                                                                  <  7:321021/89‑1 (MQ=255)                                                       gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCCCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGACCGAACAGGCGGCCTATTTAAGTc                                                                                       >  5:422457/1‑148 (MQ=38)                                                          aGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                    <  1:328279/149‑1 (MQ=255)                                                             aaaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTgc                                                                                                                                                  <  7:364486/83‑1 (MQ=255)                                                             aaaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTgc                                                                                                                                                  >  8:364486/1‑83 (MQ=255)                                                               aaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAtt                                                                              <  3:78272/149‑1 (MQ=255)                                                               aaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAA‑CCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAtt                                                                              >  6:159424/1‑148 (MQ=255)                                                                  ggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCa                                                                            >  5:304192/1‑148 (MQ=255)                                                                     aTCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGAGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGc                                                                        >  1:292518/1‑149 (MQ=255)                                                                           ggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgct                                                                  <  3:90771/149‑1 (MQ=255)                                                                               gATCCTGTGGTGGTATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgctgctg                                                              <  4:151115/149‑1 (MQ=255)                                                                                 tCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                    >  5:320245/1‑141 (MQ=255)                                                                                 tCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                    <  6:320245/141‑1 (MQ=255)                                                                                 tCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCa                                                            <  1:119511/149‑1 (MQ=255)                                                                                  tcGGTGGTTGGATGTCGGTGCCTGCCTCTTCGCGGGTAAACCGGCGCCTACAGTTAC‑TCTGCAGCACCTGTATTCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCAc                                                           <  8:213670/148‑1 (MQ=38)                                                                                    ggtggtGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGAAGCACCTGTAGGAGCGGTTTTACCAGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAATCAAAGCACTGCTGCTGCAc                                                           >  7:213670/1‑147 (MQ=38)                                                                                        gtggGATGTCGGGGCCTGCCTCTTCGCGGTTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTg                                                     <  5:255611/149‑1 (MQ=255)                                                                                        gtggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCt                                                      <  6:304192/148‑1 (MQ=255)                                                                                        gtggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTg                                                     >  8:123370/1‑149 (MQ=255)                                                                                                tCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCCCGCTTTTAACCGCGAAAAGGCCGGAAAAGGCAGCCTTTTTAAGTCAAAGCAAATCAAAGCTCTGATGCTGCACCGCCTTGGAATcc                                              >  8:146265/1‑148 (MQ=38)                                                                                                ggggggCCTGCCGCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGTATCCa                                             <  6:422457/147‑1 (MQ=255)                                                                                                  ggggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                    <  2:86139/108‑1 (MQ=255)                                                                                                  ggggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAg                                            <  7:146265/148‑1 (MQ=255)                                                                                                  ggggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTACAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                    >  1:86139/1‑108 (MQ=255)                                                                                                   gggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGcc                                          >  4:305680/1‑149 (MQ=255)                                                                                                          cctCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTg                                   <  1:325544/149‑1 (MQ=255)                                                                                                            tctTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                   >  5:318396/1‑83 (MQ=255)                                                                                                            tctTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                   <  6:318396/83‑1 (MQ=255)                                                                                                                 gctggTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTg                            <  6:241552/146‑1 (MQ=255)                                                                                                                         aCCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTcccc                    <  3:305680/149‑1 (MQ=255)                                                                                                                          cccGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt                   <  2:228735/149‑1 (MQ=255)                                                                                                                             gcgcCTACAGGTGC‑TGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGt                >  2:207236/1‑149 (MQ=37)                                                                                                                               gcCTACAGGTGC‑TGCAGAGTACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                             <  1:367620/134‑1 (MQ=37)                                                                                                                               gcCTACAGGTGC‑TGCAGAGTACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                             >  2:367620/1‑134 (MQ=37)                                                                                                                                 cTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg            >  5:213965/1‑149 (MQ=255)                                                                                                                                 cTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg            <  7:273317/149‑1 (MQ=255)                                                                                                                                   aCAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc           >  2:145875/1‑148 (MQ=255)                                                                                                                                        tAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTg                  >  6:392302/1‑136 (MQ=255)                                                                                                                                        tAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTg                  <  5:392302/136‑1 (MQ=255)                                                                                                                                            tctGCCGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTg  >  7:16511/1‑149 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAcc                                     <  6:103277/114‑1 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCa          >  5:351607/1‑141 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCa          <  6:351607/141‑1 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGACTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAcc                                     >  5:103277/1‑114 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAGGTGGCCTTTTTCATTCAAAGCAATTCAAAGCTCTCCAGCTGCACAGCCTGGGAATCCAGCCCTACATGGACATTGAACTCCCCTGGTTCGGCAACCCGcag  >  8:204840/1‑147 (MQ=255)                                                                                                                                             |                                                                                                                                                    TGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTG  >  NC_002947/1598692‑1598978

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