Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,219	G→C	12.5%	intergenic (+59/+5)	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,219	0	G	C	12.5%	172.8 / 14.5	64	intergenic (+59/+5)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base C (3/5);  total (26/38)
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.19e-01

AGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGG  >  NC_002947/895072‑895359                                                                                                                                                     |                                                                                                                                              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                                                                                                                                      <  6:41596/113‑1 (MQ=37)                                             gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTcaca                                                                                                                                      >  5:41596/1‑113 (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                                                                                                                                     <  6:168750/73‑1 (MQ=17)                                                                                      cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAg                                                                                                                                     >  5:168750/1‑73 (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                                        <  3:263750/147‑1 (MQ=255)                                                                                                        gATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCg                                        >  4:263750/1‑147 (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)                                                                                                                                                     |                                                                                                                                              AGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAA‑CTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGG  >  NC_002947/895072‑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: ‑