Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,833	G→A	7.3%	intergenic (‑72/+59)	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,833	0	G	A	7.3%	195.6 / 3.9	68	intergenic (‑72/+59)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base A (3/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598692‑1598981                                                                                                                                               |                                                                                                                                                     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                                                                                                                       >  2:119511/1‑148 (MQ=255)                          cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                       >  1:330151/1‑148 (MQ=255)                          cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaa                                                                                                                       >  4:78272/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                                                                                                                                          <  8:129116/121‑1 (MQ=255)                                  tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAg                                                                                                                                          >  7:129116/1‑121 (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                                                                                                                                                     <  7:321021/89‑1 (MQ=255)                                                       gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTgc                                                                                                                                                     >  8:321021/1‑89 (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                                                                       <  6:320245/141‑1 (MQ=255)                                                                                 tCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                       >  5:320245/1‑141 (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                                                                                                      <  6:318396/83‑1 (MQ=255)                                                                                                            tctTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                      >  5:318396/1‑83 (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                                >  2:367620/1‑134 (MQ=37)                                                                                                                               gcCTACAGGTGC‑TGCAGAGTACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                                <  1:367620/134‑1 (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             <  6:351607/141‑1 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCa             >  5:351607/1‑141 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGACTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAcc                                        >  5:103277/1‑114 (MQ=255)                                                                                                                                            tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAGGTGGCCTTTTTCATTCAAAGCAATTCAAAGCTCTCCAGCTGCACAGCCTGGGAATCCAGCCCTACATGGACATTGAACTCCCCTGGTTCGGCAACCCGcag     >  8:204840/1‑147 (MQ=255)                                                                                                                                               gcagcaCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCaa  >  4:113924/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     TGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTAC‑TCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598692‑1598981

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