Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,695,789	C→A	6.0%	A165A (GCG→GCT)	PP_4155 ←	amine oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,695,789	0	C	A	6.0%	209.3 / 3.2	66	A165A (GCG→GCT)	PP_4155	amine oxidase
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base A (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CGATGGCGAACTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTGTGCGTGC  >  NC_002947/4695641‑4695934                                                                                                                                                     |                                                                                                                                                   cGATGGCGAACTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGcc                                                                                                                                                    >  5:110072/1‑149 (MQ=255)   aTGGCGAACTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGc                                                                                                                                                  <  2:299858/149‑1 (MQ=255)   aTGGCGAACTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGc                                                                                                                                                  <  1:198715/149‑1 (MQ=255)         aaCTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAg                                                                                                                                            >  6:62567/1‑149 (MQ=255)           cTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCg                                                                                                                                          <  4:113216/149‑1 (MQ=255)              cGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCggg                                                                                                                                        >  5:135463/1‑148 (MQ=255)              cGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCggg                                                                                                                                        >  1:178706/1‑148 (MQ=255)                tGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCaa                                                                                                                                     <  2:178706/149‑1 (MQ=255)                    aCGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTg                                                                                                                                  >  3:16525/1‑148 (MQ=255)                      gcgcgAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGaaa                                                                                                                               >  7:103501/1‑149 (MQ=255)                        gcgAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATc                                                                                                                             <  3:64022/149‑1 (MQ=255)                         cgAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCg                                                                                                                            <  6:105964/149‑1 (MQ=255)                         cgAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCg                                                                                                                            <  2:80474/149‑1 (MQ=255)                                   gTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCTTGCAGCGGCACGCggcggg                                                                                                                                     >  3:5886/1‑126 (MQ=255)                                   gTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCTTGCAGCGGCACGCggcggg                                                                                                                                     <  4:5886/130‑5 (MQ=255)                                   gTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACg                                                                                                                   >  8:259212/1‑148 (MQ=255)                                    tCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGca                                                                                                                 <  5:103979/149‑1 (MQ=255)                                      gAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGcaa                                                                                                                <  4:16525/148‑1 (MQ=255)                                             gCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCTTGCAGCGGCACGCggcgg                                                                                                                                      <  1:296237/119‑4 (MQ=255)                                             gCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCTTGCAGCGGCACGCggcgg                                                                                                                                      >  2:296237/1‑116 (MQ=255)                                               cATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCAt                                                                                                      >  1:240488/1‑149 (MQ=255)                                               cATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCAt                                                                                                      >  4:32145/1‑149 (MQ=255)                                               cATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCAt                                                                                                      >  4:35993/1‑149 (MQ=255)                                                      ccATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAg                                                                                               >  2:230294/1‑149 (MQ=255)                                                       cATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGt                                                                                              <  7:259212/149‑1 (MQ=255)                                                        atcgatcgCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAagccaagcca                                                                                                       >  5:38303/1‑139 (MQ=255)                                                        atcgatcgCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAagccaagcca                                                                                                       <  6:38303/139‑1 (MQ=255)                                                                   cAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGc                                                                                  >  3:237097/1‑148 (MQ=255)                                                                     caCAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCaa                                                                               >  2:69295/2‑149 (MQ=255)                                                                      aCAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGAAACAG‑CGGCCa                                                                                >  5:291436/1‑147 (MQ=255)                                                                       cAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGCAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATg                                                                             <  1:230294/149‑1 (MQ=255)                                                                          aGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCACCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGt                                                                          >  1:111893/1‑149 (MQ=255)                                                                           ggTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCaa                                                                                          >  6:237061/1‑133 (MQ=255)                                                                               gttgCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTtctc                                                                      <  6:110072/148‑1 (MQ=255)                                                                                ttgCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCgg                                                                    <  8:103501/149‑1 (MQ=255)                                                                                ttgCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCg                                                                     <  6:135463/148‑1 (MQ=255)                                                                                  gCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑Cg                                                                                    >  2:161848/1‑131 (MQ=255)                                                                                  gCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑Cg                                                                                    <  1:161848/131‑1 (MQ=255)                                                                                          agaCGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTg                                                          <  3:35993/149‑1 (MQ=255)                                                                                          agaCGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTg                                                          <  5:62567/149‑1 (MQ=255)                                                                                          agaCGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTg                                                          <  3:32145/149‑1 (MQ=255)                                                                                          agaCGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCACGCCGCGTGTCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGATCTCGGCAAAATGCt                                                           >  6:255687/1‑148 (MQ=255)                                                                                           gaCGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTg                                                          >  1:141339/1‑148 (MQ=255)                                                                                               ccTGGTCAAGGCATCCCAGGCGGCGTCGTTATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGcc                                                     <  2:111893/149‑1 (MQ=255)                                                                                               ccTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGc                                                      >  2:84007/1‑148 (MQ=255)                                                                                                      aaGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGc                                              >  3:56740/1‑149 (MQ=255)                                                                                                           aTCCCATGCGGCGTCGATATCCGTCAACGCCCCGCCGCGTGCCGCTTCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCAATGTTGAAGTGCAACAG‑CGGCCAATGTGTTCTCGGCAAAATTCTGAAGCCATGCCGCACCtt                                         <  6:291436/149‑1 (MQ=255)                                                                                                              ccAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCAt                                      >  7:85948/1‑149 (MQ=255)                                                                                                                  gcggcgTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGc                                  <  2:240488/149‑1 (MQ=255)                                                                                                                      cgTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTc                              <  8:85948/149‑1 (MQ=255)                                                                                                                        tCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCa                                                    <  6:95082/125‑1 (MQ=255)                                                                                                                          gATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAGCCAGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCAc                            <  1:84007/148‑1 (MQ=255)                                                                                                                             atCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTc                                        <  3:12806/132‑1 (MQ=255)                                                                                                                             atCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTc                                        >  4:12806/1‑132 (MQ=255)                                                                                                                              tCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTTTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCAcc                      <  4:237097/149‑1 (MQ=255)                                                                                                                              tCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCAcc                      >  8:123701/1‑149 (MQ=255)                                                                                                                                cGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGACAccc                     >  4:208387/1‑148 (MQ=255)                                                                                                                                cGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATAGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGACAccc                     >  8:18509/1‑148 (MQ=255)                                                                                                                                  gCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAAc                  <  4:19437/149‑1 (MQ=255)                                                                                                                                  gCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAAc                  >  1:54947/1‑149 (MQ=255)                                                                                                                                       gCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGa              >  5:7386/1‑148 (MQ=255)                                                                                                                                        ccccGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTTAAGCCATGCCGCAGCTTCATCTGCTTTCACGCCACCCAACTCGACa            <  1:69295/149‑1 (MQ=255)                                                                                                                                           cgccgcGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGCAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAgtg         <  1:71590/149‑1 (MQ=255)                                                                                                                                                cgTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTgtgcgt    <  4:56740/149‑1 (MQ=255)                                                                                                                                                cgTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGATGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTgtgcgt    >  8:46335/1‑149 (MQ=255)                                                                                                                                                 gTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTgtgcgtg   >  5:283741/1‑149 (MQ=255)                                                                                                                                                  tGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTgtgcgtgc  >  3:101348/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGATGGCGAACTGCGTGGCACGCGCGAGCGCGGGGTCGAGTGTGGCCATGTAGCCATCGATCGCGTCAGACAAAGGTTGTTGCCGATCGAGACGCCTGGTCAAGGCATCCCAGGCGGCGTCGATATCCGGCAACGCCCCGCCGCGTGCCGCTGCAAGCGGGCAACTGAAATCGGTCTCGACGCAAGCCAAGCCATGTTGAAGTGCAACAG‑CGGCCAATGGGTTCTCGGCAAAATGCTGAAGCCATGCCGCACCTTCATCTGCTTTCACGCCACCCAACTCGACAGTGTGCGTGC  >  NC_002947/4695641‑4695934

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