Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	88,653	T→A	6.3%	Q222L (CAG→CTG)	betC ←	choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	88,653	0	T	A	6.3%	231.3 / 2.6	74	Q222L (CAG→CTG)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (29/40);  new base A (2/3);  total (31/43)
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.00e+00

CAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTAC  >  NC_002947/88508‑88799                                                                                                                                                  |                                                                                                                                                    cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                                     <  1:409596/146‑1 (MQ=255) cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTgg                                                                                                                                                   <  5:8816/148‑1 (MQ=255) cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTgg                                                                                                                                                   <  4:60086/148‑1 (MQ=255) cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTg                                                                                                                                                    <  3:75067/147‑1 (MQ=255) cagcaATTGCCCGATGTAGTCATCGATATAGCTGCACGCGCCGAAGTAGGCCCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGATCCAAAGCTCGAAGACAAACAGCAGGCCCTGTGAATGCGGGTCCAGTTCTGCCAGGc                                                                                                                                                  >  6:250665/1‑149 (MQ=255)  agcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGCCGAGTTCAGCCAGGc                                                                                                                                                  >  2:385539/1‑148 (MQ=255)    caATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTg                                                                                                                                                >  8:158902/1‑148 (MQ=255)          gCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCgg                                                                                                                                         >  7:377412/1‑149 (MQ=255)          gCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCg                                                                                                                                          <  5:8986/148‑1 (MQ=255)           ccGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCAGCCAGGCAGAACTCgac                                                                                                                                        <  5:404682/149‑3 (MQ=255)           ccGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCAGCCAGGCAGAACTCgac                                                                                                                                        <  8:390362/149‑3 (MQ=255)            cGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGCCGAGTTCAGCCAGGCAGAACTCgacc                                                                                                                                       <  1:385539/149‑4 (MQ=255)            cGATGTTGTCATCAATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑Ca                                                                                                                                      >  4:125165/1‑149 (MQ=255)                gttgtCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACggg                                                                                                                                  >  5:81299/1‑149 (MQ=255)                 ttgtCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACggg                                                                                                                                  <  7:381050/148‑1 (MQ=255)                      aTCGATTTAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATTGCGAACCTTATCCAGAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            <  7:152353/149‑1 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                           >  8:67486/1‑135 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                           <  7:67486/135‑1 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            <  7:221284/149‑1 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            >  2:96090/1‑149 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            <  1:401705/149‑1 (MQ=255)                      aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATg                                                                                                                             <  2:80016/148‑1 (MQ=255)                        cGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            >  2:383203/1‑147 (MQ=255)                        cGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            <  1:383203/147‑1 (MQ=255)                            ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGatat                                                                                                                      >  4:362795/1‑149 (MQ=255)                               gCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATcc                                                                                                                    <  3:362795/148‑1 (MQ=255)                                    aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCt                                                                                                              >  3:392080/1‑149 (MQ=255)                                      gcgcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTACCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTc                                                                                                            <  8:99626/149‑1 (MQ=255)                                        gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            <  3:74816/131‑1 (MQ=255)                                        gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATgg                                                                                                                            >  4:74816/1‑131 (MQ=255)                                        gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCg                                                                                                           >  1:200722/1‑148 (MQ=255)                                         ccGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                         <  8:377412/149‑1 (MQ=255)                                               tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTc                                                                                                   <  3:125165/149‑1 (MQ=255)                                                 ggCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGATCTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                                  <  1:286461/148‑1 (MQ=255)                                                    gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAg                                                                                               >  3:272821/1‑148 (MQ=255)                                                    gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGt                                                                                              <  5:13456/149‑1 (MQ=255)                                                     cggcggGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                             <  6:200729/149‑1 (MQ=255)                                                     cggcggGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                             <  2:200722/149‑1 (MQ=255)                                                            cATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGCCCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTACGTAGCGGTCCCAGTAGCgagtg                                                                                      >  6:261096/1‑145 (MQ=255)                                                             aTCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                                      <  1:96090/149‑1 (MQ=255)                                                              tCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTgg                                                                                     >  7:234937/1‑148 (MQ=255)                                                                gcgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCaa                                                                                  <  6:414671/149‑1 (MQ=255)                                                                gcgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCaa                                                                                  >  5:414671/1‑149 (MQ=255)                                                                  gAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                <  4:392080/149‑1 (MQ=255)                                                                  gAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                <  6:81299/149‑1 (MQ=255)                                                                    aCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                              <  3:199680/149‑1 (MQ=255)                                                                       ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgtg                                                                             >  1:24013/1‑147 (MQ=255)                                                                            cacaGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgggg                                                                       <  6:143534/148‑1 (MQ=255)                                                                                     tGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGc                                                                                                                                                  <  4:338586/65‑1 (MQ=255)                                                                                     tGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGc                                                                                                                                                  >  3:338586/1‑65 (MQ=255)                                                                                     tGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggt                                                             <  7:162053/149‑1 (MQ=255)                                                                                       gCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                           <  7:158902/149‑1 (MQ=255)                                                                                         ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGc                                                                 <  8:240411/141‑1 (MQ=255)                                                                                         ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGc                                                                 >  7:240411/1‑141 (MQ=255)                                                                                          tGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgggt                                                         >  1:186773/1‑148 (MQ=255)                                                                                            ttCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCAt                                                      >  8:258587/1‑149 (MQ=255)                                                                                            ttCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCAt                                                      >  5:155010/1‑149 (MQ=255)                                                                                              ccAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATCCGGGTGGGTCATc                                                     >  5:33484/1‑148 (MQ=255)                                                                                                aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGa                                                   >  1:120425/1‑148 (MQ=255)                                                                                                aaaCGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACATTCGTAGCGGTCCCAGTAGCGTTTGGCATTGGTGTAGTGGTCATGCGCGAGGGTCATCGa                                                   >  4:95558/1‑148 (MQ=255)                                                                                                     tCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGt                                             <  2:21580/149‑1 (MQ=255)                                                                                                       gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCca                                                                                                                   <  6:232763/77‑1 (MQ=255)                                                                                                       gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGt                                             >  4:200879/1‑147 (MQ=255)                                                                                                       gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                           <  2:120425/149‑1 (MQ=255)                                                                                                                cagcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGTAAAAt                                   <  2:24013/148‑1 (MQ=255)                                                                                                                   cagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggcc                               <  8:234937/149‑1 (MQ=255)                                                                                                                          tgtgAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATc                        >  2:13383/1‑149 (MQ=255)                                                                                                                                 gCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcg                 >  3:372101/1‑149 (MQ=255)                                                                                                                                  cGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc                <  2:186773/149‑1 (MQ=255)                                                                                                                                   gggTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgcg               <  5:261096/149‑1 (MQ=255)                                                                                                                                      tCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGca             >  5:322763/1‑148 (MQ=255)                                                                                                                                          gTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCTCACATgg        >  7:164441/1‑149 (MQ=255)                                                                                                                                                ccTGGCTGAACTCTG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggt                                                             <  2:193958/90‑1 (MQ=255)                                                                                                                                                ccTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCTCATGGTCGTAc  <  7:258587/149‑1 (MQ=255)                                                                                                                                                ccTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGt    <  4:372101/147‑1 (MQ=255)                                                                                                                                                  tGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc                <  4:298494/133‑1 (MQ=255)                                                                                                                                                  tGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc                >  3:298494/1‑133 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGG‑CACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTAC  >  NC_002947/88508‑88799

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