Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,791,105	C→T	6.4%	G183G (GGC→GGT)	PP_3353 →	sulfatase‑modifying protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,791,105	0	C	T	6.4%	233.5 / 3.3	77	G183G (GGC→GGT)	PP_3353	sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (41/31);  new base T (2/3);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGG  >  NC_002947/3790957‑3791250                                                                                                                                                       |                                                                                                                                                  cGGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGc                                                                                                                                                   <  6:10581/149‑1 (MQ=255) cGGACGTGCTCGGCCCCGGCCGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGc                                                                                                                                                   <  3:272780/149‑1 (MQ=255)      gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTcgcgg                                                                                                                                              >  8:351703/1‑148 (MQ=255)        gCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTAGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAAc                                                                                                                                            >  6:349720/1‑149 (MQ=255)          tCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACAATCCGGTGGTGCAGATCG‑CCCTGGAACATGCCCAGGCCTATGCC‑CGCTGCGCAGGCCGCGAACt                                                                                                                                           >  4:110185/1‑148 (MQ=255)             gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGccc                                                                                                                                       <  5:132872/149‑1 (MQ=255)             gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGcc                                                                                                                                        >  8:14858/1‑148 (MQ=255)              ccccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGcc                                                                                                                                        >  2:110879/1‑147 (MQ=255)                ccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGc                                                                                                                                    >  8:299355/1‑149 (MQ=255)                ccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGACTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGc                                                                                                                                    <  8:175264/149‑1 (MQ=255)                  ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                  <  7:146807/149‑2 (MQ=255)                  ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                  <  7:351703/149‑2 (MQ=255)                         aGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAACATGCCCAGGCCTATGCCTCGCTGGGCAGTTCGCCGCCTGCCCAGCgggcatag                                                                                                                            <  5:23427/149‑8 (MQ=38)                          gTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                          >  8:143926/1‑149 (MQ=255)                             cGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                        <  1:229055/148‑1 (MQ=255)                             cGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                        <  7:299355/148‑1 (MQ=255)                                gCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGaa                                                                                                                     >  5:149869/1‑148 (MQ=255)                                 cccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATa                                                                                                                   >  8:201963/1‑149 (MQ=255)                                   cGGCGCCAACTGGCGACACCCGCCAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGATCTGCCCAGCGAGGCGCAGCTGGAATACg                                                                                                                 <  1:369313/149‑1 (MQ=255)                                           aCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                         <  1:330240/149‑1 (MQ=255)                                               gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCa                                                                                                                                      >  7:180119/1‑116 (MQ=255)                                               gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCa                                                                                                                                      <  8:180119/116‑1 (MQ=255)                                               gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                      <  2:187862/148‑1 (MQ=255)                                                 gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCggcggc                                                                                                   <  7:245745/149‑1 (MQ=255)                                                    acCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCt                                                                                                 >  4:10446/1‑148 (MQ=255)                                                    acCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTg                                                                                                >  5:339467/1‑149 (MQ=255)                                                      ccGCAAGGGCCGGGCAGCAGCCTGGCTGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                              >  1:218837/1‑149 (MQ=255)                                                      ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGTGCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGa                                                                                               >  1:108546/1‑149 (MQ=255)                                                      ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                              >  5:83217/1‑149 (MQ=255)                                                            gggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGc                                                                                        <  7:143926/149‑1 (MQ=255)                                                             ggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                       >  5:254918/1‑149 (MQ=255)                                                                      gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGc                                                                              >  6:192585/1‑149 (MQ=255)                                                                        agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTg                                                                            <  6:116168/149‑1 (MQ=255)                                                                        agcCTGGCCGGGCTAGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTg                                                                            <  5:349720/149‑1 (MQ=255)                                                                         gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                           <  4:167784/149‑1 (MQ=255)                                                                         gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                           <  6:149869/149‑1 (MQ=255)                                                                          ccTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                          <  7:201963/149‑1 (MQ=255)                                                                          ccTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                          >  6:351465/1‑149 (MQ=255)                                                                           cTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgggg                                                                         >  8:339903/1‑149 (MQ=255)                                                                              gCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTac                                                                      >  3:192199/1‑149 (MQ=255)                                                                                cGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTacca                                                                    <  2:108546/149‑1 (MQ=255)                                                                                  ggCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                  <  3:10446/149‑1 (MQ=255)                                                                                    cTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                  >  6:148964/1‑147 (MQ=255)                                                                                       gACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGag                                                               >  4:136304/1‑147 (MQ=255)                                                                                       gACAACCATCCGGTGGAGCAGATCG‑CCCTGGAAGATGCCCACGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGTAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagca                                                             >  8:87149/1‑149 (MQ=255)                                                                                         cAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcagc                                                           <  1:110879/149‑1 (MQ=255)                                                                                            ccATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGccc                                                         >  1:193328/1‑148 (MQ=255)                                                                                            ccATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGccc                                                         >  5:19604/1‑148 (MQ=255)                                                                                            ccATCCGGTCGTGCAGATCG‑CCCTGGAAGATGCCCATGCATATACC‑CGCTCGGCAGGCCCCCAACTGCACAGCGAGCCGCATCTGGAAAATGCCATACGCGTAGGCCTGACCGAAGACGACATAAGCTGGGGTACCACCGAGAAGccc                                                         >  6:200901/1‑148 (MQ=255)                                                                                             cATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGccc                                                         >  6:364550/1‑147 (MQ=255)                                                                                               tCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑GGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                     <  3:110185/149‑1 (MQ=255)                                                                                               tCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                     >  7:294485/1‑149 (MQ=255)                                                                                                       tgCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTc                                             <  6:339467/149‑1 (MQ=255)                                                                                                        gCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCa                                            >  3:202704/1‑149 (MQ=255)                                                                                                        gCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCa                                            >  5:303465/1‑149 (MQ=255)                                                                                                        gCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCa                                            >  7:130189/1‑149 (MQ=255)                                                                                                         cAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                           <  5:364550/149‑1 (MQ=255)                                                                                                          aGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATg                                          <  6:83217/149‑1 (MQ=255)                                                                                                          aGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATg                                          <  5:351465/149‑1 (MQ=255)                                                                                                            aTCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                        <  7:14858/149‑1 (MQ=255)                                                                                                                 cccTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCaa                                     >  6:139160/1‑148 (MQ=255)                                                                                                                  ccTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATa                                   <  7:339903/149‑1 (MQ=255)                                                                                                                              ccAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                        >  5:334408/1‑148 (MQ=255)                                                                                                                              ccAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                        >  3:307980/1‑148 (MQ=255)                                                                                                                              ccAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                        >  4:279480/1‑148 (MQ=255)                                                                                                                               cAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                        >  2:240007/1‑147 (MQ=255)                                                                                                                               cAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                      <  8:130189/149‑1 (MQ=255)                                                                                                                                 ggCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTTGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCGGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGt                    <  6:270924/149‑1 (MQ=255)                                                                                                                                  gCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCACTACCTGGCAGTGTCAGtt                   >  8:142128/1‑149 (MQ=255)                                                                                                                                    cTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAATGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                 <  5:139160/149‑1 (MQ=255)                                                                                                                                     tATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                >  4:85006/1‑149 (MQ=255)                                                                                                                                       tGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                <  6:19740/147‑1 (MQ=255)                                                                                                                                           tcGCTGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGCCCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc           >  8:100364/2‑149 (MQ=255)                                                                                                                                               tGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCaa       >  1:359635/1‑149 (MQ=255)                                                                                                                                                  gCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGACGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgc    >  4:213679/1‑149 (MQ=255)                                                                                                                                                   cAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcg   <  3:279480/149‑1 (MQ=255)                                                                                                                                                   cAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcg   <  3:136304/149‑1 (MQ=255)                                                                                                                                                    aGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcgg  >  5:133544/1‑149 (MQ=255)                                                                                                                                                    aGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTAACGCAATgcgg  >  2:314072/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CGGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCG‑CCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGG  >  NC_002947/3790957‑3791250

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