Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,152	A→T	6.1%	L264Q (CTG→CAG)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,152	0	A	T	6.1%	191.7 / 2.9	66	L264Q (CTG→CAG)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (28/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 = 7.48e-01

TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAG  >  NC_002947/5181006‑5181291                                                                                                                                                     |                                                                                                                                            tATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagcag                                                                                                                                            >  1:314407/1‑148 (MQ=255)     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                                                                                                                  >  6:395496/1‑117 (MQ=255)                                                          gCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                                  <  5:395496/117‑1 (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                                                          <  6:163561/148‑1 (MQ=255)                                                                                   gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                          <  1:225685/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                                                                                                                           <  3:160153/68‑1 (MQ=255)                                                                                                  ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAatat                                                                                                                           >  4:160153/1‑68 (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           >  7:93185/1‑148 (MQ=255)                                                                                                                                  tttCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacaccac           >  6:326655/1‑148 (MQ=255)                                                                                                                                           acGAG‑CATCATCACGGCAGCCTAAATATAGAATTACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTTTAGTCGGCAAAGGTCCTCACGTAGTAGGCATATCCGAGGGAGGCGTAGATCCACACCACAACCGCCAg  <  2:77702/148‑1 (MQ=255)                                                                                                                                            cGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAg  >  6:42239/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                            TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAG  >  NC_002947/5181006‑5181291

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