Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,146	A→T	6.1%	L266H (CTC→CAC)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,146	0	A	T	6.1%	195.6 / 2.8	66	L266H (CTC→CAC)	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.22e-01

ACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAG  >  NC_002947/5181002‑5181291                                                                                                                                                  |                                                                                                                                                   aCGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTCACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGACGAAGTCTTTCACGCCCGGGTTCTTGCATTCGCCTGACATGTGCTCGAACACCGCATTCATTTCGGC‑GCCGag‑cag                                                                                                                                               >  4:70406/1‑149 (MQ=255)  cGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagc                                                                                                                                              <  3:33545/149‑1 (MQ=255)  cGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGag‑cagc                                                                                                                                              <  7:220966/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:375381/149‑1 (MQ=255)                                                             ggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGcc                                                                                   <  3:108720/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                                                          <  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           >  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)                                                                                                                                                  |                                                                                                                                                   ACGCTATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAG‑CAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAG  >  NC_002947/5181002‑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: ‑