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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	87,959	0	G	C	7.0%	268.1 / 9.1	85	R453R (CGC→CGG)	betC	choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (45/34);  new base C (4/2);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGA  >  NC_002947/87812‑88104                                                                                                                                                    |                                                                                                                                                     ggcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGgcggcggcg                                                                                                                                                    >  2:996156/1‑150 (MQ=255)    gggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCAGGAGGCGGCgctg                                                                                                                                                 >  2:101823/1‑150 (MQ=255)     ggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                <  1:867548/150‑1 (MQ=255)        tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                <  2:430808/147‑1 (MQ=255)        tgtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                                >  1:430808/1‑147 (MQ=255)            gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCCGACTGCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCa                                                                                                                                      >  1:1071392/1‑150 (MQ=255)              tAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGt                                                                                                                                    >  2:388039/1‑150 (MQ=255)                      cTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgctg                                                                                                                            >  2:411309/1‑150 (MQ=255)                       ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTgg                                                                                                                           >  1:310358/1‑150 (MQ=255)                        tGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTATTGGCTGGCATCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGgc                                                                                                                          >  1:1038711/1‑150 (MQ=255)                                 ggTCATCGAGATCGAGGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 <  2:249702/150‑1 (MQ=255)                                     aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgc                                                                                                                              >  1:679814/1‑133 (MQ=255)                                     aTCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACctgc                                                                                                                              <  2:679814/133‑1 (MQ=255)                                      tCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATg                                                                                                            >  1:814478/1‑150 (MQ=255)                                      tCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATg                                                                                                            >  2:411766/1‑150 (MQ=255)                                       cGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCACTACCTGCTGGCGCAGGCTGGTGATGt                                                                                                           <  2:15883/150‑1 (MQ=255)                                              gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAgcg                                                                                                    >  2:1281674/1‑150 (MQ=255)                                                   ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:1320878/150‑1 (MQ=255)                                                   ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:1209866/150‑1 (MQ=255)                                                   ggTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                               <  1:411766/150‑1 (MQ=255)                                                    gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 <  1:884128/131‑1 (MQ=255)                                                    gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                                 >  2:884128/1‑131 (MQ=255)                                                         gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                       >  2:1154447/1‑136 (MQ=255)                                                         gcATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCa                                                                                                       <  1:1154447/136‑1 (MQ=255)                                                          cATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTc                                                                                        <  2:846704/150‑1 (MQ=255)                                                          cATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTc                                                                                        <  1:682178/150‑1 (MQ=255)                                                              tACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         <  1:551449/97‑1 (MQ=255)                                                              tACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         >  2:551449/1‑97 (MQ=255)                                                               aCTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCGTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGa                                                                                   >  1:608516/1‑150 (MQ=255)                                                                 tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTTAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  1:625750/1‑150 (MQ=255)                                                                 tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  1:898696/1‑150 (MQ=255)                                                                 tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                                 >  2:216860/1‑150 (MQ=255)                                                                 tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         <  1:806926/94‑1 (MQ=38)                                                                 tGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCGccggct‑‑‑gg                                                                                                                                         >  2:806926/1‑94 (MQ=38)                                                                   ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggct‑‑‑ggc                                                                                                                                        <  1:761152/93‑1 (MQ=255)                                                                   ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggct‑‑‑ggc                                                                                                                                        >  2:761152/1‑93 (MQ=255)                                                                        tggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          >  1:138474/1‑150 (MQ=255)                                                                        tggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          <  1:996156/150‑1 (MQ=255)                                                                        aggcGTGCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCt                                                                          >  2:1065482/2‑150 (MQ=255)                                                                                accaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcagc                                                                  >  2:405431/1‑150 (MQ=255)                                                                                 ccaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcagca                                                                 >  1:1220419/1‑150 (MQ=255)                                                                                  caGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCAc                                                                >  2:1092827/1‑150 (MQ=255)                                                                                    gTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                              <  2:17521/102‑1 (MQ=255)                                                                                    gTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGa                                                                                                              >  1:17521/1‑102 (MQ=255)                                                                                     tggttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                             <  2:608516/150‑1 (MQ=255)                                                                                       gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggt                                                           >  2:99451/1‑150 (MQ=255)                                                                                        ttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtggtg                                                          <  2:138474/150‑1 (MQ=255)                                                                                              gtCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                    <  1:216860/150‑1 (MQ=255)                                                                                               tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTAATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:946141/1‑150 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTa                                                                                                                                   <  2:1388833/68‑1 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTa                                                                                                                                   >  1:1388833/1‑68 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                      >  1:120161/1‑113 (MQ=255)                                                                                                 ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCAt                                                                                      <  2:120161/113‑1 (MQ=255)                                                                                                  cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:882845/1‑147 (MQ=255)                                                                                                  cAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   <  2:882845/147‑1 (MQ=255)                                                                                                   aGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCGGTTGCGCTGCATCGACAAAAGCCTGCAGCAGCACCTGGTGGTCCGGGCTgt                                               >  2:654163/1‑149 (MQ=255)                                                                                                      ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCgg                                            >  2:355166/1‑150 (MQ=255)                                                                                                        ttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTg                                          >  1:45853/1‑150 (MQ=255)                                                                                                         tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                         <  1:411309/150‑1 (MQ=255)                                                                                                         tCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                         <  2:1038711/150‑1 (MQ=255)                                                                                                          cAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCAAGCGCTATTGCGCTTCATCGACAAAGACCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                        <  1:388039/150‑1 (MQ=255)                                                                                                           aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAgg                                       >  2:1348175/1‑150 (MQ=255)                                                                                                                 ccGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTcc                                 <  2:1241011/150‑1 (MQ=255)                                                                                                                  cGATGGCCAGCGCTTCGAAAACAAGCCGGCGGCGGCGCTGTCT‑‑‑GTCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCCGCGATGTTACGCTTCATCGACAAATGTCGCGATCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTccc                                >  2:1020983/1‑150 (MQ=255)                                                                                                                   gATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCg                               >  1:707732/1‑150 (MQ=255)                                                                                                                    aTGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGc                              <  2:814478/150‑1 (MQ=255)                                                                                                                        ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                            <  2:894137/116‑1 (MQ=255)                                                                                                                        ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                            >  1:894137/1‑116 (MQ=255)                                                                                                                        ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGt                          >  2:25721/1‑150 (MQ=255)                                                                                                                         cAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTg                         >  2:1130791/1‑150 (MQ=255)                                                                                                                          aGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                        <  1:1224803/150‑1 (MQ=255)                                                                                                                          aGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                        <  2:898696/150‑1 (MQ=255)                                                                                                                           gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                    >  1:755243/1‑137 (MQ=255)                                                                                                                           gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTtc                                    <  2:755243/137‑1 (MQ=255)                                                                                                                                   gCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCgg               <  2:1071392/150‑1 (MQ=255)                                                                                                                                     ggccagccagCGCCGCCGCCGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCt             >  2:276604/3‑150 (MQ=255)                                                                                                                                     caccaGCCGGCGGCGGCGCTGGCTGGCGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg                >  2:1258255/1‑150 (MQ=255)                                                                                                                                     caccaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCt             >  1:53860/1‑150 (MQ=255)                                                                                                                                     caccaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCACGCTCGTGCGGGGCGCGGCt             >  2:1327261/1‑150 (MQ=255)                                                                                                                                      accaGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTc            >  2:915190/1‑150 (MQ=255)                                                                                                                                        cagccagcGCCGCCGCCGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAg          >  1:130311/1‑150 (MQ=255)                                                                                                                                         aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   >  1:1263987/1‑108 (MQ=255)                                                                                                                                         aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCggg                                                   <  2:1263987/108‑1 (MQ=255)                                                                                                                                         aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         <  2:583917/150‑1 (MQ=255)                                                                                                                                         aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         <  1:1281674/150‑1 (MQ=255)                                                                                                                                         aGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAgg         >  1:1261863/1‑150 (MQ=255)                                                                                                                                            cggcggcggCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa      >  1:1150388/1‑150 (MQ=255)                                                                                                                                              gcggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa    <  1:654163/150‑1 (MQ=255)                                                                                                                                                ggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATaga  <  2:310358/150‑1 (MQ=255)                                                                                                                                                ggcggcgCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATaga  <  1:355166/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCT‑‑‑GGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGA  >  NC_002947/87812‑88104

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