Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,216	G→A	12.5%	intergenic (+56/+8)	PP_0777 → / ← PP_0778	glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	895,216	0	G	A	12.5%	166.4 / 11.8	64	intergenic (+56/+8)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (22/34);  new base A (3/5);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

GGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGG  >  NC_002947/895068‑895359                                                                                                                                                      |                                                                                                                                                 ggCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGgcggcg                                                                                                                                                  >  3:143682/1‑149 (MQ=255)  gCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGgcgc                                                                                                                                                 >  3:9380/1‑149 (MQ=255)     aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                              <  4:143682/149‑1 (MQ=255)      ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGa                                                                                                                                             <  7:20960/149‑1 (MQ=255)       gATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                            <  1:129431/149‑1 (MQ=39)           gCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                        <  8:120432/149‑1 (MQ=255)            ccGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                        >  4:115360/1‑148 (MQ=255)             cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                        >  3:102416/1‑147 (MQ=255)              gggTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTGCAGCGCCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTcaca                                                                                                                                      <  1:68043/149‑1 (MQ=25)              gggTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAg                                                                                                                                     <  5:42127/149‑1 (MQ=38)              gggTGCTGGCGCGGTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAg                                                                                                                                     <  3:197505/149‑1 (MQ=255)                   cTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCt                                                                                                                                <  6:238705/149‑1 (MQ=255)                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                           <  6:264486/147‑1 (MQ=255)                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                           >  5:264486/1‑147 (MQ=255)                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACa                                                                                                                         >  4:61451/1‑149 (MQ=255)                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACa                                                                                                                         >  1:246747/1‑149 (MQ=255)                           tttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAt                                                                                                                        >  4:31872/1‑149 (MQ=255)                            ttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAtt                                                                                                                       >  4:215799/1‑149 (MQ=255)                                   cgcgCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAg                                                                                                                 >  4:156525/1‑148 (MQ=255)                                        aCGGCGCCGGATGACCCAGCGGTCGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACgac                                                                                                           >  7:229740/1‑147 (MQ=255)                                          ggCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGcc                                                                                                           >  1:213277/1‑147 (MQ=255)                                            cgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGc                                                                                                       <  1:278107/149‑1 (MQ=255)                                             gcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTTATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCCGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGcc                                                                                                      >  6:326711/1‑149 (MQ=255)                                                 gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTcaca                                                                                                                                      >  5:41596/1‑113 (MQ=37)                                                 gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTcaca                                                                                                                                      <  6:41596/113‑1 (MQ=37)                                                      cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGGACCa                                                                                             <  6:78203/149‑1 (MQ=255)                                                      cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAc                                                                                               >  5:95839/1‑147 (MQ=255)                                                       ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                            <  3:31872/149‑1 (MQ=255)                                                       ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                            <  3:112906/149‑1 (MQ=255)                                                              tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                                     <  2:246747/149‑1 (MQ=255)                                                                    aGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgc                                                                               >  5:76043/1‑149 (MQ=255)                                                                     ggCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcg                                                                              <  3:215799/149‑1 (MQ=255)                                                                       ccATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcg                                                                              >  7:222707/1‑147 (MQ=255)                                                                       ccATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcg                                                                              <  8:222707/147‑1 (MQ=255)                                                                          tCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                         <  3:156525/149‑1 (MQ=255)                                                                          tCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                         <  5:134427/149‑1 (MQ=255)                                                                            gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTc                                                                        <  6:95839/148‑1 (MQ=255)                                                                            gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCa                                                                       <  4:57117/149‑1 (MQ=255)                                                                              gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                         >  3:303625/1‑145 (MQ=255)                                                                              gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGa                                                                     <  4:303625/149‑1 (MQ=255)                                                                                aaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGa                                                                     >  4:102469/1‑147 (MQ=255)                                                                                    gatGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcc                                                                 <  7:47719/145‑1 (MQ=255)                                                                                          cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAg                                                                                                                                     >  5:168750/1‑73 (MQ=17)                                                                                          cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAg                                                                                                                                     <  6:168750/73‑1 (MQ=17)                                                                                           tGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑Ct                                                       >  1:296164/1‑149 (MQ=255)                                                                                             acgagTGCCTGACGTGATGGTA‑TGGCACCTGTGAGCTCGAGGGCCGCCCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑Ct                                                       <  8:203259/145‑1 (MQ=39)                                                                                              agagTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGtt                                                    <  2:213277/149‑1 (MQ=255)                                                                                              agagTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCCGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTTCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGtt                                                    <  5:326711/149‑1 (MQ=255)                                                                                               gagTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTc                                                   >  7:29713/1‑149 (MQ=255)                                                                                                         cGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGc                                          <  1:96212/148‑1 (MQ=255)                                                                                                            gATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCg                                        >  4:263750/1‑147 (MQ=255)                                                                                                            gATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCg                                        <  3:263750/147‑1 (MQ=255)                                                                                                             aTGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGcc                                     >  5:310829/1‑149 (MQ=255)                                                                                                                      gCAGCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTgg                             <  3:230841/149‑1 (MQ=255)                                                                                                                       cACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCTGCATGCGCTTTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTggc                            <  8:229740/149‑1 (MQ=255)                                                                                                                           tgtgAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGGTGC‑CAAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                       >  1:129126/1‑149 (MQ=25)                                                                                                                             tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCt                        <  4:9380/147‑1 (MQ=255)                                                                                                                             tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                       <  4:102416/148‑1 (MQ=255)                                                                                                                               agaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                    >  5:222056/1‑149 (MQ=255)                                                                                                                                gaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACCTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                    <  1:51647/149‑1 (MQ=255)                                                                                                                                 aTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt                   <  6:310829/148‑1 (MQ=255)                                                                                                                                    gAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTg               >  3:78764/1‑149 (MQ=255)                                                                                                                                      gTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTc             <  3:189056/149‑1 (MQ=255)                                                                                                                                         ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTg                                                                                        <  6:252536/72‑1 (MQ=255)                                                                                                                                         ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc           <  8:29713/148‑1 (MQ=255)                                                                                                                                           gcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCa        <  4:78764/149‑1 (MQ=255)                                                                                                                                                 gcggcgCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTAAGATACCCATCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAgg  >  6:149159/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGG  >  NC_002947/895068‑895359

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