Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,141:1	+T	5.8%	coding (802/957 nt)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,141	1	.	T	5.8%	281.6 / 11.3	69	coding (802/957 nt)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base . (30/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

TTCGGCTCGACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCAC  >  NC_002947/5180993‑5181282                                                                                                                                                      |                                                                                                                                               ttCGGCTCGACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGgc                                                                                                                                                 >  3:210154/1‑149 (MQ=255)     gCTCGACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGgctgcc                                                                                                                                             >  1:233051/1‑145 (MQ=255)          aCGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTCACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGACGAAGTCTTTCACGCCCGGGTTCTTGCATTCGCCTGACATGTGCTCGAACACCGCATTCATTTCGGC‑GCCGag‑cag                                                                                                                                      >  4:70406/1‑149 (MQ=255)           cGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagc                                                                                                                                     <  7:220966/149‑1 (MQ=255)           cGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagc                                                                                                                                     <  3:33545/149‑1 (MQ=255)            gCTATGAGGCTTGCGCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagc                                                                                                                                     <  3:70406/148‑1 (MQ=255)              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:130139/1‑149 (MQ=255)                                                                                                          ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                      >  7:51207/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)                                                                                                                                                      |                                                                                                                                               TTCGGCTCGACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCAC  >  NC_002947/5180993‑5181282

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