Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,155	A→T	6.1%	V263E (GTG→GAG)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,155	0	A	T	6.1%	198.0 / 2.8	66	V263E (GTG→GAG)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

AGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCG  >  NC_002947/5181010‑5181302                                                                                                                                                    |                                                                                                                                                    aGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACgg                                                                                                                                                  >  2:237263/1‑149 (MQ=255) aGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTGTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACg                                                                                                                                                   <  2:327259/148‑1 (MQ=255)     ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGc                                                                                                                                              >  7:340903/1‑149 (MQ=255)     ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCAGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGc                                                                                                                                              >  4:396241/1‑149 (MQ=255)       gCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCTATAACCGCATTAATTTCGGC‑GCCGAG‑CAGCATCACGGCac                                                                                                                                               >  8:168994/1‑145 (MQ=255)           gggTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGaaa                                                                                                                                         <  2:436785/148‑1 (MQ=255)           gggTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGaa                                                                                                                                          <  4:164471/147‑1 (MQ=255)             gTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTG‑CTGCTGCTCGGC‑GCCGAAatga                                                                                                                                      <  2:233051/149‑3 (MQ=38)             gTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAata                                                                                                                                       >  5:360922/1‑148 (MQ=255)                tGGGGGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTTATTTCGGCTGCCGTG‑CTGCTGCTCGGC‑GCCGAAatgaat                                                                                                                                    <  8:349026/148‑5 (MQ=37)                 gggTGCTCGTGACCCAGTACAGTCAAGGTATCCTTCGGCTTATCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGGACGATCACCGCATTCATTTCTGC‑GCCGAG‑CACCAGCACGGCAGCCGAAAtatacg                                                                                                                                   >  8:113013/1‑146 (MQ=255)                      cTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATac                                                                                                                              <  3:168057/148‑1 (MQ=255)                          tGACACAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCAGTCGAAGTCTTTCTCGCCCGGGTTCTTTCCTTCGGCTGACATTTGCTCGATCACCGCATTCATTTCTGA‑GCCGAG‑CAGACGCACGGCAGCCGAAATATAGAAATacagcaa                                                                                                                         >  7:139240/1‑149 (MQ=255)                            aCCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATacagcaac                                                                                                                        <  8:340903/148‑1 (MQ=255)                                 gTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAAtata                                                                                                                                     <  7:220263/130‑1 (MQ=255)                                 gTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCAc                                                                                                                   >  8:77351/1‑148 (MQ=255)                                    cAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACga                                                                                                                 <  1:90636/148‑1 (MQ=255)                                       tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                            >  1:77702/1‑149 (MQ=255)                                       tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                            >  8:46339/1‑149 (MQ=255)                                        cAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGAATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                            <  3:349201/148‑1 (MQ=255)                                            gTTTCCTTGGGCTTTTCCCCGACGAAGTCTTTCTCGCCCGGGTTCATGCCTTCGGCTGACAAGTGCTCGAACACCGCATTCATTTCGGC‑GCCGAG‑CAGCACCACGGCACCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                        >  1:67451/1‑148 (MQ=255)                                              ttCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCaa                                                                                                     >  7:171823/1‑149 (MQ=255)                                                ccctGGGCTTCTCCCCGCCGAACGCTCTCTCGCCCGCGTGCGCGCCTTCGGCTGACATGGGCGCGAGCACCGCATTCATTTCCGC‑GCCGAG‑CAGCACCACGGCAGCCGAAATAAAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                                    <  4:315058/145‑1 (MQ=255)                                                ccTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                                    <  7:77351/148‑1 (MQ=255)                                                 cTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGT‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtg                                                                                                   <  5:314028/148‑1 (MQ=255)                                                    gggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑ACCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACGGCAACAGCACGATGATTGCACCAAtgctgc                                                                                               >  3:372245/1‑149 (MQ=255)                                                     ggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGcc                                                                                              <  3:108720/149‑1 (MQ=255)                                                     ggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGcc                                                                                              <  3:375381/149‑1 (MQ=255)                                                      gCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                                             <  5:395496/117‑1 (MQ=255)                                                      gCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                                             >  6:395496/1‑117 (MQ=255)                                                            cccaGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                                        <  3:396241/148‑1 (MQ=255)                                                             cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATg                                                                                       <  8:171823/148‑1 (MQ=255)                                                                  cGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACAAATGCTGCCGTACATGGCGTTg                                                                                 >  1:311774/1‑149 (MQ=255)                                                                    aaGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTa                                                                               <  6:438871/149‑1 (MQ=255)                                                                           tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCActa                                                                                                                 <  2:113974/108‑3 (MQ=255)                                                                           tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                                          <  7:46339/147‑1 (MQ=255)                                                                           tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGc                                                                         <  6:360922/148‑1 (MQ=255)                                                                           tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTGCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACga                                                                                                                 >  1:113974/1‑108 (MQ=255)                                                                               gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAgg                                                                    >  5:439646/1‑149 (MQ=255)                                                                               gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                                     <  1:225685/148‑1 (MQ=255)                                                                               gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                                     <  6:163561/148‑1 (MQ=255)                                                                               gCCAGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCAGGGCAGCCGAATTATAGAAATACAGCAACAGCACGATGATTGCACCAATTCTGCCTTACATGGCGTTGTAGTCGGCAAAg                                                                     <  1:36383/148‑1 (MQ=255)                                                                                      ttCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCAc                                                             <  1:143456/149‑1 (MQ=255)                                                                                      ttCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCAc                                                             <  1:237263/149‑1 (MQ=255)                                                                                         ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                          >  7:51207/1‑149 (MQ=255)                                                                                         ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                          >  7:130139/1‑149 (MQ=255)                                                                                             cTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtaggag                                                      <  6:143601/149‑1 (MQ=255)                                                                                              ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAatat                                                                                                                                      >  4:160153/1‑68 (MQ=255)                                                                                              ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAatat                                                                                                                                      <  3:160153/68‑1 (MQ=255)                                                                                              ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACttagtatg                                                     >  4:35049/1‑147 (MQ=255)                                                                                                cGGCTGACTTGTGCTCGATCACCGGATTCATTTCGGC‑GGCGAG‑CAGCAGCACGGCAGCCGAAATATAGATTTTCAGCAAGAGGAGGATGATTGCACCAATGCTGGCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGc                                                    <  7:113013/148‑1 (MQ=255)                                                                                                cGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGc                                                    >  2:18996/1‑148 (MQ=255)                                                                                                   cTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTc                                                                            >  1:387788/1‑121 (MQ=255)                                                                                                   cTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTc                                                                            <  2:387788/121‑1 (MQ=255)                                                                                                          gtgCTCGATGAGCGCATTCATTTCGGC‑GCCGAT‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCATCAGCACGATGTTTGCTCCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAg                                           <  2:67451/147‑1 (MQ=255)                                                                                                                gATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCa                                   >  2:195499/1‑149 (MQ=255)                                                                                                                gATCACCGCATTCATTTCGGC‑ACCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACGGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGc                                     <  4:372245/147‑1 (MQ=255)                                                                                                                   cACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGa                                >  5:428266/1‑149 (MQ=255)                                                                                                                          ttcattTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATacagcaaca                                                                                                                       <  1:398215/55‑1 (MQ=255)                                                                                                                          ttcattTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATacagcaaca                                                                                                                       >  2:398215/1‑55 (MQ=255)                                                                                                                          tatattTCGGCTGCCGTG‑CTGCTGCTCGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                         >  2:272612/4‑149 (MQ=38)                                                                                                                              tttCGGCTGCCGTG‑CTGCTGCTCGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacaccac                      >  4:348586/1‑148 (MQ=255)                                                                                                                              tttCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacaccac                      >  6:326655/1‑148 (MQ=255)                                                                                                                              tttCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacaccac                      >  7:93185/1‑148 (MQ=255)                                                                                                                                       acGAG‑CATCATCACGGCAGCCTAAATATAGAATTACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTTTAGTCGGCAAAGGTCCTCACGTAGTAGGCATATCCGAGGGAGGCGTAGATCCACACCACAACCGCCAg             <  2:77702/148‑1 (MQ=255)                                                                                                                                        cGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAg             >  6:42239/1‑148 (MQ=255)                                                                                                                                                   caCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCg  <  6:428266/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCG  >  NC_002947/5181010‑5181302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑