Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,450,217	G→A	10.3%	S32N (AGC→AAC)	ybeY →	endoribonuclease YbeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,450,217	0	G	A	10.3%	184.4 / 9.0	68	S32N (AGC→AAC)	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base A (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 = 6.51e-01

CCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGCTGCTGGG  >  NC_002947/5450072‑5450364                                                                                                                                                   |                                                                                                                                                    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                                                                                                                              <  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                                                                            <  3:104075/149‑1 (MQ=255)                                                                       aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGcac                                                                            <  5:118267/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                      >  5:193579/1‑149 (MQ=255)                                                                                                                             aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg                      >  4:374098/1‑149 (MQ=255)                                                                                                                             aaTTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTg                      >  8:305566/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               >  4:162383/1‑149 (MQ=255)                                                                                                                                    gTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGata               >  2:404918/1‑149 (MQ=255)                                                                                                                                       gcgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATc             <  2:48285/148‑1 (MQ=255)                                                                                                                                       gcgc‑CAGCGCAGCGCCGACTCAGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATc             <  8:265496/148‑1 (MQ=255)                                                                                                                                        cgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATcc            >  3:45347/1‑148 (MQ=255)                                                                                                                                        cgc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATcc            <  4:45347/148‑1 (MQ=255)                                                                                                                                         gc‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCg          >  4:285559/1‑149 (MQ=255)                                                                                                                                                 cagcgcCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGCTGCTgg   <  3:285559/149‑1 (MQ=255)                                                                                                                                                  agcgcCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGCTGCTggg  >  8:396938/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGC‑CAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGCTGCTGGG  >  NC_002947/5450072‑5450364

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