Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,677,630	G→C	6.4%	intergenic (+23/‑108)	thrC → / → PP_5471	threonine synthase/two‑component system response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,677,630	0	G	C	6.4%	198.2 / 3.7	62	intergenic (+23/‑108)	thrC/PP_5471	threonine synthase/two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (2/2);  total (29/33)
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.35e-01

AGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGCC  >  NC_002947/1677486‑1677774                                                                                                                                                  |                                                                                                                                                 aGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTgcagcg                                                                                                                                                >  1:55543/1‑147 (MQ=255)    tACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGcc                                                                                                                                           >  6:314671/1‑149 (MQ=255)        cgcgCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCaa                                                                                                                                       <  4:293840/149‑1 (MQ=255)         gcgcACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAAc                                                                                                                                      <  8:341747/149‑1 (MQ=255)             aCCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtct                                                                                                                                  <  8:117746/149‑1 (MQ=255)             aCCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtct                                                                                                                                  <  5:176734/149‑1 (MQ=255)                tCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAg                                                                                                                               >  7:349179/1‑149 (MQ=255)                               gCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGt                                                                                                                <  2:47432/149‑1 (MQ=255)                                ccGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTc                                                                                                               <  1:196242/149‑1 (MQ=255)                                ccGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAACGCTGTACAGGGCTTAGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTc                                                                                                               >  3:40676/1‑149 (MQ=255)                                 cGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTc                                                                                                               >  7:67450/1‑148 (MQ=255)                                 cGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTc                                                                                                               >  2:335984/1‑148 (MQ=255)                                 cGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCt                                                                                                              <  4:86402/149‑1 (MQ=255)                                          cGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAAc                                                                                                     <  6:70020/149‑1 (MQ=255)                                           gTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACa                                                                                                    >  7:151810/1‑149 (MQ=255)                                            ttGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACAt                                                                                                   <  8:349179/149‑1 (MQ=255)                                            ttGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACAt                                                                                                   <  5:374222/149‑1 (MQ=255)                                              gCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATc                                                                                                  <  5:314671/148‑1 (MQ=255)                                                    tCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAtt                                                                                             >  7:308662/1‑147 (MQ=255)                                                      cTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGc                                                                                         <  2:55543/149‑1 (MQ=255)                                                           cAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                    >  1:101783/1‑149 (MQ=255)                                                            aaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCg                                                                                   <  2:101783/149‑1 (MQ=255)                                                               gACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                                >  3:325512/1‑149 (MQ=255)                                                                 ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCAGTGGCGTTTGGGGCTGCAGCGCAGCCACAACGTCTCAGCCTTTTCCCTCCTGTCAATCAGAAACATAAAATTCAGCCTCCCGGgcg                                                                               >  1:303323/1‑149 (MQ=255)                                                                 ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgcgc                                                                              >  6:211340/1‑149 (MQ=255)                                                                    gtaGGCTGTACAGTGCTTTGTCATCCAGCACGGCAACCGCGGCAAGGGGCTGTAATTGCA‑TGGCGGTTGTGGCTGCGCTGCAGCGCCaaac                                                                                                                                     <  4:51695/89‑3 (MQ=255)                                                                    gAATGCTGTACAGGGCTTTGTCAGCCAGCACCGCAACCGCGGCAAGCCGCTGAAATTGCA‑AGGCGTTTGGGGCTGCAGCGCAGCCCCAACCTCTCAGCCATTTCCCTCCAGTCTATCAGAAACATCCAATTCAGACAGCCGGGcccc                                                                             >  5:176455/1‑147 (MQ=255)                                                                    gAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGGTTGGGGCTGCGCTGCAGCCCCaaac                                                                                                                                     >  3:51695/1‑89 (MQ=38)                                                                     aaGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCaa                                                                           <  5:211340/148‑1 (MQ=255)                                                                         cTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTa                                                                                               <  6:74608/124‑1 (MQ=255)                                                                         cTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTa                                                                                               >  5:74608/1‑124 (MQ=255)                                                                             ctatctcTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAAt                                                                   <  5:214877/142‑1 (MQ=255)                                                                                gggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                               >  1:355786/1‑149 (MQ=255)                                                                                 ggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGcc                                                              >  2:282957/1‑149 (MQ=255)                                                                                 ggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGcc                                                              <  3:113131/149‑1 (MQ=255)                                                                                    tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAAt                                                                   >  6:214877/1‑141 (MQ=255)                                                                                        tCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcctc                                                       <  5:229043/149‑1 (MQ=255)                                                                                        tCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcctc                                                       >  6:229043/1‑149 (MQ=255)                                                                                         cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcctcc                                                      >  2:352424/1‑149 (MQ=255)                                                                                               gcACGGCAACCGCGGCAAGCCGCTGTAATTGCT‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAg                                                <  5:263611/149‑1 (MQ=255)                                                                                                  cGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAg                                             <  4:325512/149‑1 (MQ=255)                                                                                                      aaCCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                                         <  4:66554/149‑1 (MQ=255)                                                                                                      aaCCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                                         >  5:94891/1‑149 (MQ=255)                                                                                                           ctgCAAGCCGCGCTATTTGCT‑TGGCGGTTGGGGCTGCTGGGCAGCCCCTACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCATGCATAATTTAGCCGGTCCTCCCATCAGGAGCGGCTCATg                                    <  7:301823/147‑1 (MQ=255)                                                                                                            ggCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGc                                   <  4:40676/149‑1 (MQ=255)                                                                                                                 gCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGGGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAg                                                <  2:177866/131‑1 (MQ=255)                                                                                                                 gCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAg                                                >  1:177866/1‑131 (MQ=255)                                                                                                                        ttaTTTCA‑TGGGGGTTGGGGGTGCAGGGCAGGGCCAACGTGTCTGCCTTTTGCCTCCTGGCTATGAGTAACATCCTATTGAGGCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt                        <  2:303323/146‑1 (MQ=255)                                                                                                                          aTTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAtgtg                     <  8:151810/149‑1 (MQ=255)                                                                                                                             gCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGtt                   <  6:94891/148‑1 (MQ=255)                                                                                                                              cA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGa                 >  2:216625/1‑149 (MQ=255)                                                                                                                               a‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAt                >  5:65862/1‑149 (MQ=255)                                                                                                                               a‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAt                >  5:249245/1‑149 (MQ=255)                                                                                                                                gagacGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAt                >  4:2346/5‑149 (MQ=255)                                                                                                                                 tGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAt                >  2:323858/1‑148 (MQ=255)                                                                                                                                   gCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTc              <  1:282957/148‑1 (MQ=255)                                                                                                                                     gTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCAcc           <  1:323858/149‑1 (MQ=255)                                                                                                                                         ggggCTGCGCTGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCgg                                          >  2:4056/1‑114 (MQ=255)                                                                                                                                         ggggCTGCGCTGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCgg                                          <  1:4056/114‑1 (MQ=255)                                                                                                                                           ggCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTc     >  2:89133/1‑149 (MQ=255)                                                                                                                                            gCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCg    >  7:95495/1‑149 (MQ=255)                                                                                                                                             cTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGc   <  8:95495/149‑1 (MQ=255)                                                                                                                                              tGCAGCGCAGCCCCAACGTGTCAGCCTTTTCCCTCCTGTCTTTCTGATAGATCCTATTCAGCCTGGCGTGCGCCATGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGGTCATGCACCAGCCCTATGTGTTTATTCACCAGGGTGGcc  <  6:176455/149‑1 (MQ=255)                                                                                                                                                cagcgcagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGTCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGcc  <  4:173223/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 AGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCA‑TGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGCC  >  NC_002947/1677486‑1677774

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