Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,450,204	T→C	7.4%	L28L (TTG→CTG)	ybeY →	endoribonuclease YbeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,450,204	0	T	C	7.4%	190.2 / 3.7	67	L28L (TTG→CTG)	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATA  >  NC_002947/5450057‑5450351                                                                                                                                                    |                                                                                                                                                     aGCGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGtt                                                                                                                                                     <  5:25486/149‑1 (MQ=255)   cGTATCGTCGAAGCCTACGACCGCTTCGATGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgt                                                                                                                                                   <  1:431461/149‑2 (MQ=255)   cGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgc                                                                                                                                                   <  6:378201/149‑1 (MQ=255)   cGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgc                                                                                                                                                   >  5:322748/1‑149 (MQ=255)   cGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTgc                                                                                                                                                   >  5:364774/1‑149 (MQ=255)     tATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgcgc                                                                                                                                                 >  5:434828/1‑149 (MQ=255)     tATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgcgc                                                                                                                                                 <  3:430142/149‑1 (MQ=255)            gAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑cagcg                                                                                                                                           >  2:60954/1‑147 (MQ=255)                ccTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGC‑Tggcgcaacgc                                                                                                                                     <  6:364774/149‑1 (MQ=255)                  tACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGcc                                                                                                                                    >  1:29729/1‑148 (MQ=255)                    cGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGCCGTCCGGCAAAGATGCTTGAACTTGATCTCCCACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGC‑TGGCGCAACGccaa                                                                                                                                  <  5:322358/148‑3 (MQ=255)                     gACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACt                                                                                                                                >  6:269820/1‑149 (MQ=255)                     caCCGCTTCCAAGCCCGCCAGCCCAAGCCCGCGGCGGCCGGCAATGAGGCTTCAAATGGATATCCAACGGGCCACGCATGCTGCAACT‑CAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGC‑TGGCGCAACGccaattc                                                                                                                               <  3:304535/148‑6 (MQ=255)                       ccGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCg                                                                                                                              <  1:310652/149‑1 (MQ=255)                       ccGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGCCCACGGATGCTGCCACTCCACAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCg                                                                                                                              >  6:118267/1‑149 (MQ=255)                        cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCgg                                                                                                                             <  2:294218/149‑1 (MQ=255)                        cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCgg                                                                                                                             <  4:420032/149‑1 (MQ=255)                             cGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAAt                                                                                                                         <  7:37057/148‑1 (MQ=255)                              gAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGa                                                                                                                       <  2:402355/149‑1 (MQ=255)                                aGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGAcc                                                                                                                     <  6:434828/149‑1 (MQ=255)                                    cGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTc                                                                                                                 >  3:124693/1‑149 (MQ=255)                                    cGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTc                                                                                                                 <  6:53055/149‑1 (MQ=255)                                      ccAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCACCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGc                                                                                                               <  8:31678/149‑1 (MQ=255)                                         gCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTg                                                                                                            >  5:226016/1‑149 (MQ=255)                                          cccAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGCCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCCCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTgg                                                                                                           >  6:4901/1‑149 (MQ=255)                                              aGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTcga                                                                                                       <  7:253382/149‑1 (MQ=255)                                                 ccGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTcgacg                                                                                                     >  2:328440/1‑148 (MQ=255)                                                           cGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGc                                                                                          >  6:370412/1‑149 (MQ=255)                                                                aaGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCCGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGa                                                                                     <  1:329889/149‑1 (MQ=255)                                                                  gATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGc                                                                                   <  2:29729/149‑1 (MQ=255)                                                                    tGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCt                                                                                  <  4:124693/148‑1 (MQ=255)                                                                     gCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTg                                                                                 >  4:230382/1‑148 (MQ=255)                                                                         gAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACca                                                                            >  1:212362/1‑149 (MQ=255)                                                                            cTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCacac                                                                         <  6:322748/149‑1 (MQ=255)                                                                               gATATCCAACGGGCCATGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTa                                                                      >  1:78246/1‑149 (MQ=255)                                                                               gATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTa                                                                      >  4:104075/1‑149 (MQ=255)                                                                               gATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATAACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTa                                                                      >  7:21130/1‑149 (MQ=255)                                                                                  atCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCg                                                                   >  2:256590/1‑149 (MQ=255)                                                                                      aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGACGCCGAAGGCCGCGAGCTGAACCACACCTACCGGca                                                                >  6:123641/1‑148 (MQ=255)                                                                                      aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGcac                                                               <  5:118267/149‑1 (MQ=255)                                                                                      aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGcac                                                               <  3:104075/149‑1 (MQ=255)                                                                                        cGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACaa                                                             <  2:58871/149‑1 (MQ=255)                                                                                         gggCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAg                                                            >  5:235222/1‑149 (MQ=255)                                                                                           gCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGa                                                          >  3:404397/1‑149 (MQ=255)                                                                                               cGGACGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTAc                                                      >  4:286612/1‑149 (MQ=255)                                                                                                 gATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTAcgc                                                    >  8:288745/1‑149 (MQ=255)                                                                                                    gctgcCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGAc                                                 >  3:179435/1‑149 (MQ=255)                                                                                                       gcCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCAGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCaa                                              >  7:265496/1‑149 (MQ=255)                                                                                                         cACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCGCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAAt                                             >  1:102181/1‑149 (MQ=255)                                                                                                         cACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAAtg                                            <  1:328440/149‑1 (MQ=255)                                                                                                               agaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGtgttgt                                      >  1:48285/1‑149 (MQ=255)                                                                                                                gaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTc                                     <  6:226016/149‑1 (MQ=255)                                                                                                                 tctaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGtgttgt                                      <  7:326765/144‑1 (MQ=255)                                                                                                                     cgccgcCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGtgttgt                                      >  8:326765/1‑143 (MQ=255)                                                                                                                          ccTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGcc                           <  3:230382/149‑1 (MQ=255)                                                                                                                          ccTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGcc                           <  7:288745/149‑1 (MQ=255)                                                                                                                           cTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCg                          <  5:4901/149‑1 (MQ=255)                                                                                                                               cgccgcTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGAGTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGAtgt                      <  5:370412/149‑1 (MQ=255)                                                                                                                                       gTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGa                <  4:179435/147‑1 (MQ=255)                                                                                                                                       gTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATg              <  2:102181/149‑1 (MQ=255)                                                                                                                                          cgagtcggCGCTGCGC‑TGGCGCAACGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCCCAAGGACTACGCGACCAAAGAGTTGTCGTTCCCGGCCGCTGTGCCCGATGAc            >  8:209364/7‑148 (MQ=255)                                                                                                                                            aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg         >  4:374098/1‑149 (MQ=255)                                                                                                                                            aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg         >  8:305566/1‑149 (MQ=255)                                                                                                                                            aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg         >  5:193579/1‑149 (MQ=255)                                                                                                                                             aTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGc        >  1:267400/1‑149 (MQ=255)                                                                                                                                                  cGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGat   <  8:21130/149‑1 (MQ=255)                                                                                                                                                   gTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGata  >  2:404918/1‑149 (MQ=255)                                                                                                                                                   gTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGata  >  4:162383/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGCGTATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATA  >  NC_002947/5450057‑5450351

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