Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,584,214	T→G	8.9%	E457A (GAG→GCG)	glgE ←	alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,584,214	0	T	G	8.9%	165.1 / 6.1	56	E457A (GAG→GCG)	glgE	alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base G (2/3);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTGCGC  >  NC_002947/4584066‑4584360                                                                                                                                                     |                                                                                                                                                    ccACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAgctgct                                                                                                                                                     >  2:149585/1‑149 (MQ=255)  cacaAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTc                                                                                                                                                    <  2:665/149‑1 (MQ=255)         cTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGaatatt                                                                                                                                             >  5:230268/1‑145 (MQ=255)          tGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGaatattt                                                                                                                                            >  5:139898/1‑144 (MQ=255)             gCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgc                                                                                                                                         <  8:278113/149‑1 (MQ=255)               ccATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                                                                                        >  3:24668/1‑148 (MQ=255)                cATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                                                                                        >  4:133604/1‑147 (MQ=255)                 aTGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAg                                                                                                                                     <  3:277588/149‑1 (MQ=255)                    gTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑Tc                                                                                                                                 >  3:271510/1‑149 (MQ=255)                     tCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑Tc                                                                                                                                 <  5:65040/148‑1 (MQ=255)                       gCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGC‑Tcaat                                                                                                                              <  6:230268/149‑4 (MQ=255)                        ccAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGc                                                                                                                             <  3:133604/149‑1 (MQ=255)                         cAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGC‑Tcaatca                                                                                                                            <  6:139898/149‑6 (MQ=255)                          aGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGC‑Tcaatcag                                                                                                                           <  5:262980/149‑7 (MQ=255)                                 cACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTg                                                                                                                    >  3:281741/1‑149 (MQ=255)                                   cGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTgcg                                                                                                                  <  2:111726/149‑1 (MQ=255)                                         aGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAgg                                                                                                             <  3:274784/148‑1 (MQ=255)                                          ggAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTg                                                                                                           <  3:291624/149‑1 (MQ=255)                                             aaCCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGaaa                                                                                                        <  4:156693/149‑1 (MQ=255)                                                   cgcgGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGt                                                                                                   <  4:271510/148‑1 (MQ=255)                                                      ggCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAg                                                                                                              <  5:56034/134‑1 (MQ=255)                                                      ggCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAg                                                                                                              >  6:56034/1‑134 (MQ=255)                                                        ccTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGct                                                                                             >  6:220597/1‑149 (MQ=255)                                                         cTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAAt                                                                                                       <  2:101304/138‑1 (MQ=255)                                                         cTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCCGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAAt                                                                                                       >  1:101304/1‑138 (MQ=255)                                                           gAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtg                                                                                          >  6:41199/1‑149 (MQ=255)                                                            aGGTGTGCAGGAAAAACGGGTTGATGTCCTGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtgt                                                                                         >  5:5883/1‑149 (MQ=255)                                                            aGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtgt                                                                                         <  5:41199/149‑1 (MQ=255)                                                              gtgtGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGc                                                                                       <  3:306474/149‑1 (MQ=255)                                                                   cAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGc                                                                                  <  4:281741/149‑1 (MQ=255)                                                                          atcGGGTTGATGTCCGGTGTATTGACGATGAAGTTGGGCCGGTAGCTCTGGCTCCACGGCGGGTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCtt                                                                           <  4:204192/147‑1 (MQ=255)                                                                          aaCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCtt                                                                           <  5:220597/149‑1 (MQ=255)                                                                           aCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtgt                                                                                         <  4:302743/134‑1 (MQ=255)                                                                           aCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtgt                                                                                         >  3:302743/1‑134 (MQ=255)                                                                           aCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTg                                                                          >  7:161402/1‑149 (MQ=255)                                                                             gggTTGATGTCCTGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCTAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGcc                                                                        <  6:5883/149‑1 (MQ=255)                                                                              ggTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCGCCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGcc                                                                        <  1:253944/149‑1 (MQ=255)                                                                               gTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGcc                                                                        >  7:139591/1‑147 (MQ=255)                                                                                   aTGTCCGGTGTATTGAAGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTAGAAATATTCGAGCAGC‑TCCTGCTTGTGGTTGCGCCATGTGAAATAGGTATAGCTCTGAGCGTAGCCGACCTTACCCAggag                                                                  >  8:112855/1‑147 (MQ=255)                                                                                          gtgtATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCcat                                                            >  5:50148/1‑148 (MQ=255)                                                                                                    gaagaaGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCt                                                 >  7:43541/1‑149 (MQ=255)                                                                                                    gaagaaGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCt                                                 >  2:212281/1‑149 (MQ=255)                                                                                                    gaagaaGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCt                                                 >  2:19427/1‑149 (MQ=255)                                                                                                       gaagTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTgg                                              >  2:21297/1‑149 (MQ=255)                                                                                                          gTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTa                                            >  5:113923/1‑148 (MQ=255)                                                                                                                ccGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCtt                                     <  4:24668/149‑1 (MQ=255)                                                                                                                 cGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTAGCGCAGA‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCATGCACAGGCGCGCCATCATCGCAGGCTTGGTAAAGGATTc                                    >  3:58778/1‑149 (MQ=255)                                                                                                                      gCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGa                               <  8:161402/149‑1 (MQ=255)                                                                                                                      gCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGa                               <  1:149585/149‑1 (MQ=255)                                                                                                                      gCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGa                               <  1:19427/149‑1 (MQ=255)                                                                                                                       cACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCg                                   <  2:247633/144‑1 (MQ=255)                                                                                                                       cACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCg                                   >  1:247633/1‑144 (MQ=255)                                                                                                                        aCTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGTTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAgg                             <  7:48504/149‑1 (MQ=255)                                                                                                                               ccACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGa                      >  6:238484/1‑149 (MQ=255)                                                                                                                                cACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGAc                     <  6:3996/149‑1 (MQ=255)                                                                                                                                cACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGAc                     >  2:24789/1‑149 (MQ=255)                                                                                                                                 aCGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACg                    >  6:118230/1‑149 (MQ=255)                                                                                                                                                 ctgctCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTgc    <  8:43541/149‑1 (MQ=255)                                                                                                                                                   gctCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTgcg   >  4:196237/1‑148 (MQ=255)                                                                                                                                                    ctCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTgcgc  <  5:118230/148‑1 (MQ=255)                                                                                                                                                    ctCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTgcgc  <  1:172071/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGC‑TCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCTGGCTGCGC  >  NC_002947/4584066‑4584360

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