Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	87,956	G→C	8.8%	R454R (CGC→CGG)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	87,956	0	G	C	8.8%	176.2 / 6.4	57	R454R (CGC→CGG)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base C (2/3);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87812‑88102                                                                                                                                                 |                                                                                                                                                   ggcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCggcggcggc                                                                                                                                                >  8:125015/1‑149 (MQ=255)         gtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                        >  3:325023/1‑149 (MQ=255)          tgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         <  2:181645/147‑1 (MQ=255)          tgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCgccggctg                                                                                                                                        <  5:98417/148‑1 (MQ=255)            gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                      <  2:175842/148‑1 (MQ=255)            gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                      >  1:175842/1‑148 (MQ=255)              tAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAg                                                                                                                                   >  1:183423/1‑149 (MQ=255)                aCGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTCCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAg                                                                                                                                   <  5:192509/147‑1 (MQ=255)                   tGCCTGGCGTTCGAGTGCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGGTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACct                                                                                                                              <  6:349449/149‑1 (MQ=255)                      cTTGCGTTCGAGGTGATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                  <  8:305865/126‑1 (MQ=255)                      cTTGCGTTCGAGGTGATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                  >  7:305865/1‑126 (MQ=255)                       ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCTCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                           <  6:12232/132‑1 (MQ=255)                       ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCTCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                           >  5:12232/1‑132 (MQ=255)                           gTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgcgc                                                                                                                      <  2:183423/149‑1 (MQ=255)                           gTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgcg                                                                                                                       >  6:212385/1‑148 (MQ=255)                                  gTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGACCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                               >  7:69575/1‑149 (MQ=255)                                   tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                  <  4:93386/113‑1 (MQ=255)                                   tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCcggcggcg                                                                                                                                                  >  3:93386/1‑113 (MQ=255)                                   tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  4:325023/149‑1 (MQ=255)                                   tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGGCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  8:69575/149‑1 (MQ=255)                                           aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  6:224371/1‑148 (MQ=255)                                           aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  8:64034/1‑148 (MQ=255)                                              gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         <  7:256766/111‑1 (MQ=255)                                              gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         >  8:256766/1‑111 (MQ=255)                                                 gtgGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCt                                                                                                >  1:350532/1‑149 (MQ=255)                                                 gtgGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCAGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCt                                                                                                >  2:153844/1‑149 (MQ=255)                                                    gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  5:20423/1‑132 (MQ=255)                                                    gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  6:20423/132‑1 (MQ=255)                                                    gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                             <  7:125015/149‑1 (MQ=255)                                                           aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  5:182532/1‑125 (MQ=255)                                                           aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  6:182532/125‑1 (MQ=255)                                                                   ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                               >  3:199174/1‑148 (MQ=255)                                                                     ggctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACaaa                                                                             <  2:350532/148‑1 (MQ=255)                                                                      gctggcGTCCACCAGTGGTTGGTGGTCCCAGCTCTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   >  6:249644/1‑141 (MQ=255)                                                                      gctggcGTCCACCAGTGGTTGGTGGTCCCAGCTCTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   <  5:249644/141‑1 (MQ=255)                                                                       ctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGc                                                                          <  4:199174/149‑1 (MQ=255)                                                                          gcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTg                                                                       >  2:221163/1‑149 (MQ=255)                                                                                accaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcag                                                                 >  2:73962/1‑149 (MQ=255)                                                                                       gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                          <  5:224371/149‑1 (MQ=255)                                                                                               tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                  >  6:119499/1‑149 (MQ=255)                                                                                               tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                  <  5:212385/149‑1 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGTGCTGGCTGGCCAGTACCTGATGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCCCCTGGTGGTCCGGGc                                                >  1:156438/1‑149 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGc                                                >  7:235058/1‑149 (MQ=255)                                                                                                      ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                           <  6:7810/133‑1 (MQ=255)                                                                                                      ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                           >  5:7810/1‑133 (MQ=255)                                                                                                      ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCg                                           >  5:338491/1‑149 (MQ=255)                                                                                                           aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      <  5:367540/149‑1 (MQ=255)                                                                                                           aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      >  7:228663/1‑149 (MQ=255)                                                                                                                        ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCCTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCg                         >  5:240044/1‑149 (MQ=255)                                                                                                                           gcgcTTGGGCCACGAGGCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCTGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCTTGCTGCCGGTGAGGTTCTCCCGCTCGTGc                      <  5:367978/149‑1 (MQ=255)                                                                                                                           gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                      <  8:228663/149‑1 (MQ=255)                                                                                                                             gcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                    >  7:123998/1‑149 (MQ=255)                                                                                                                                   gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg              >  4:217959/1‑149 (MQ=255)                                                                                                                                   gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg              >  2:146261/1‑149 (MQ=255)                                                                                                                                      gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAAtac                                                                          >  2:220596/2‑84 (MQ=38)                                                                                                                                      gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAAtac                                                                          <  1:220596/85‑3 (MQ=38)                                                                                                                                        cagccagcGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       >  6:336553/1‑119 (MQ=255)                                                                                                                                        cagccagcGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       <  5:336553/119‑1 (MQ=255)                                                                                                                                        caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      <  7:5948/120‑1 (MQ=255)                                                                                                                                        caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      >  8:5948/1‑120 (MQ=255)                                                                                                                                        caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCa         <  8:235058/149‑1 (MQ=255)                                                                                                                                           ccGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGt      <  1:73962/149‑1 (MQ=255)                                                                                                                                             ggcggcggcGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa    <  1:221163/149‑1 (MQ=255)                                                                                                                                              gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa    >  2:107401/1‑148 (MQ=255)                                                                                                                                               cggcggcgCTGGCAGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa  <  1:153844/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87812‑88102

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