Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,450,211	A→G	10.3%	Q30R (CAG→CGG) ‡	ybeY →	endoribonuclease YbeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,450,211	0	A	G	10.3%	179.0 / 9.4	68	Q30R (CAG→CGG) ‡	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (4/3);  total (36/32)
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.47e-01

GAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCG  >  NC_002947/5450068‑5450356                                                                                                                                                 |                                                                                                                                                  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                                                                                                                                  <  4:420032/149‑1 (MQ=255)             cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCgg                                                                                                                                  <  2:294218/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                                                                                                                      <  6:53055/149‑1 (MQ=255)                         cGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTc                                                                                                                      >  3:124693/1‑149 (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                                <  7:288745/149‑1 (MQ=255)                                                                                                               ccTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGcc                                <  3:230382/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              >  5:193579/1‑149 (MQ=255)                                                                                                                                 aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg              >  8:305566/1‑149 (MQ=255)                                                                                                                                 aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg              >  4:374098/1‑149 (MQ=255)                                                                                                                                  aTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGc             >  1:267400/1‑149 (MQ=255)                                                                                                                                       cgctGCGC‑TGGCGCAACGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGa         >  1:55623/4‑148 (MQ=255)                                                                                                                                       cgctGCGC‑TGGCGCAACGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGa         >  7:45934/4‑148 (MQ=255)                                                                                                                                       cGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGat        <  8:21130/149‑1 (MQ=255)                                                                                                                                        gctGCGC‑TGGCGCAACGCCAATTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGAc                                                      <  7:201294/100‑1 (MQ=255)                                                                                                                                        gctGCGC‑TGGCGCAACGCCAATTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGAc                                                      >  8:201294/3‑102 (MQ=255)                                                                                                                                        gTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGata       >  2:404918/1‑149 (MQ=255)                                                                                                                                        gTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGata       >  4:162383/1‑149 (MQ=255)                                                                                                                                           gcgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATc     <  2:48285/148‑1 (MQ=255)                                                                                                                                           gcgc‑CAGCGCAGCGCCGACTCAGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATc     <  8:265496/148‑1 (MQ=255)                                                                                                                                            cgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATcc    <  4:45347/148‑1 (MQ=255)                                                                                                                                            cgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATcc    >  3:45347/1‑148 (MQ=255)                                                                                                                                             gc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCg  >  4:285559/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCG  >  NC_002947/5450068‑5450356

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