Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,538,381	A→G	6.3%	E153E (GAA→GAG)	PP_4871 →	lysine decarboxylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,538,381	0	A	G	6.3%	179.1 / 2.2	63	E153E (GAA→GAG)	PP_4871	lysine decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGG  >  NC_002947/5538233‑5538525                                                                                                                                                     |                                                                                                                                                  ttgccttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCgaa                                                                                                                                                   <  8:14324/149‑1 (MQ=255)    ccttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCgaagaa                                                                                                                                                <  5:226149/149‑1 (MQ=255)      ttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAgg                                                                                                                                              >  8:120518/1‑149 (MQ=255)          cGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCACTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTc                                                                                                                                          <  4:302043/149‑1 (MQ=255)            ggTGGGCTGGGTAGGCTGGAGGTACTGTTCGAAGTGTGGAGCTGGGTGCAGCTGGGCTTGCACGCCAAGCCGCTCGGGCTACTGGATGTGTACGTTTTCTACGTCTAGCTGGGCGTGTTCCTCGACCATATCGTCGAAGAAGGTTTCGt                                                                                                                                        <  6:301805/149‑1 (MQ=255)            ggTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGt                                                                                                                                        <  5:132590/149‑1 (MQ=255)                 gCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGAACTGGGGGTAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGcc                                                                                                                                   >  4:182803/1‑149 (MQ=255)                     ggTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGc                                                                                                                                >  3:125085/1‑148 (MQ=255)                            tGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCATCGAAGAAGGTTTCGTGCGCCCGCAACACCgg                                                                                                                        >  5:311863/1‑149 (MQ=255)                             ggaggaACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCggg                                                                                                                       >  3:28758/1‑149 (MQ=255)                              gaggaACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGc                                                                                                                      >  2:33134/1‑149 (MQ=255)                                ggaACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACCCCGGGCGa                                                                                                                    >  5:361153/1‑149 (MQ=255)                                          gAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTtgctgc                                                                                                           <  3:228911/148‑1 (MQ=255)                                          gAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTtgctgc                                                                                                           >  4:228911/1‑148 (MQ=255)                                          gAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTgctgc‑t                                                                                                         >  4:342256/1‑149 (MQ=255)                                              tgtgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tgggg                                                                                                     >  2:287643/1‑149 (MQ=255)                                                tgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tgggg                                                                                                     >  5:367637/1‑147 (MQ=255)                                                tgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGgca                                                                                                   <  6:367637/149‑1 (MQ=255)                                                  gACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGgcagc                                                                                                 <  7:120518/149‑1 (MQ=255)                                                   aCCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGgcagca                                                                                                >  5:102081/1‑149 (MQ=255)                                                          ggCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGAc                                                                                         <  4:125085/149‑1 (MQ=255)                                                          ggCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGAc                                                                                         <  8:256067/149‑1 (MQ=255)                                                               cTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAgct                                                                                    <  1:286556/149‑1 (MQ=255)                                                                tGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAgctg                                                                                   >  7:51526/1‑149 (MQ=255)                                                                    cTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGctgctg                                                                                >  4:364103/1‑148 (MQ=255)                                                                     tATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tggg                                                                                                      <  6:246789/125‑1 (MQ=255)                                                                     tATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tggg                                                                                                      >  5:246789/1‑125 (MQ=255)                                                                       tCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGaa                                                                             <  1:329971/148‑1 (MQ=255)                                                                       tCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAg                                                                            >  2:352369/1‑149 (MQ=255)                                                                          cGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGTAAGGGa                                                                         >  6:378428/1‑149 (MQ=255)                                                                           gCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAgg                                                                           >  4:114887/1‑146 (MQ=255)                                                                           gCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAgg                                                                           <  3:114887/146‑1 (MQ=255)                                                                            ccAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCATCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATg                                                                       <  6:311863/149‑1 (MQ=255)                                                                              aaGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGa                                                                     <  7:180012/149‑1 (MQ=255)                                                                                         gCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATAGTCGAAGAAGGTTTAGTGCGCCCGCAACACCGGGCGATGTTTCTGC‑TTGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTgg                                                          >  1:191907/1‑149 (MQ=255)                                                                                           tACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGgc                                                                                                    >  5:148717/1‑105 (MQ=255)                                                                                           tACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGgc                                                                                                    <  6:148717/105‑1 (MQ=255)                                                                                             cTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGcc                                                      >  2:155362/1‑149 (MQ=255)                                                                                                     gAACGGTTTCTACGCCAAGCGGGGCGGGTTCCTCGACCATATCGGCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTTCGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTcgcg                                               <  1:287643/149‑1 (MQ=255)                                                                                                      aaCGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCgcgc                                              <  3:182803/148‑1 (MQ=255)                                                                                                      aaCGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCGGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGcc                                             >  5:217064/1‑149 (MQ=255)                                                                                                       aCGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCt                                            <  3:364103/149‑1 (MQ=255)                                                                                                        cGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgt                                                            <  8:214330/132‑1 (MQ=255)                                                                                                        cGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTtgt                                                            >  7:214330/1‑132 (MQ=255)                                                                                                        cGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCt                                            >  6:105302/1‑148 (MQ=255)                                                                                                           tttCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGt                                        <  1:352369/149‑1 (MQ=255)                                                                                                           tttCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGGCGCGCCTAAGt                                        <  6:102081/149‑1 (MQ=255)                                                                                                             tCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCGGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTg                                       <  6:217064/148‑1 (MQ=255)                                                                                                              cTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTggg                                     >  1:37354/1‑149 (MQ=255)                                                                                                                     aaGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tgg                                                                                                       <  5:365503/76‑1 (MQ=255)                                                                                                                     aaGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑Tgg                                                                                                       >  6:365503/1‑76 (MQ=255)                                                                                                                       gCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAag                             >  4:185767/1‑148 (MQ=255)                                                                                                                                gTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTAGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgact                   >  6:297749/1‑149 (MQ=255)                                                                                                                                cttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGt                                        <  1:341682/127‑1 (MQ=255)                                                                                                                                cttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGt                                        >  2:341682/2‑128 (MQ=255)                                                                                                                                 ttcttcGACGATATGGTCGAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCctgact                   >  7:159412/1‑148 (MQ=255)                                                                                                                                    cTCGACGATATCGTCGTAGAAGGTTGCGTTCGCCCGCAACACCGGGCTATGTTGCGGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGGCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGaa                <  6:361153/148‑1 (MQ=255)                                                                                                                                    cTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGaa                <  5:159493/148‑1 (MQ=255)                                                                                                                                     tCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGc              >  5:269782/1‑149 (MQ=255)                                                                                                                                           atatCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGAc        <  3:342256/149‑1 (MQ=255)                                                                                                                                                gtcgAGGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACtttt    <  8:159412/148‑1 (MQ=255)                                                                                                                                                gtcgAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTg   >  4:299415/1‑149 (MQ=255)                                                                                                                                                 tcgAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGTAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTgg  <  5:378428/149‑1 (MQ=255)                                                                                                                                                 tcgAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTgg  <  6:269782/149‑1 (MQ=255)                                                                                                                                                 tcgAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTgg  <  2:191907/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGC‑TGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGG  >  NC_002947/5538233‑5538525

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