Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,791,148	C→T	8.6%	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,148	0	C	T	8.6%	199.5 / 8.0	70	G183G (GGC→GGT) ‡	PP_3353	sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (3/3);  total (34/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 = 7.65e-01

GGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGG  >  NC_002947/3791001‑3791293                                                                                                                                                     |                                                                                                                                                  ggACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGACGAGCTCGCCAACGATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGAC‑CGCTGGGCAGGCcg                                                                                                                                                 >  1:371725/1‑150 (MQ=255)       gCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACt                                                                                                                                           >  2:1215347/1‑150 (MQ=255)            gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACGGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCa                                                                                                                                      <  1:996646/150‑1 (MQ=255)             ccccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                     >  1:404872/1‑150 (MQ=255)              cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                    <  1:255821/150‑1 (MQ=255)                 ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                 >  1:270384/1‑150 (MQ=255)                   cTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGgc                                                                                                                               <  2:404872/150‑1 (MQ=255)                      gCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCa                                                                                                                            >  2:754798/1‑150 (MQ=255)                         gTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCt                                                                                                                         >  2:471498/1‑150 (MQ=255)                          ttCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                        <  2:270384/150‑1 (MQ=255)                                  cGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGc                                                                                                                <  1:1344687/150‑1 (MQ=255)                                        cAACTGGCGACACCCGCAAGGGCCGGGGGGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcg                                                                                                          >  1:43185/1‑150 (MQ=255)                                            tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                      >  2:237197/1‑150 (MQ=255)                                              gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                    <  1:969654/118‑1 (MQ=255)                                              gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                    >  2:969654/1‑118 (MQ=255)                                              gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGcggcgg                                                                                                    >  1:1197349/1‑150 (MQ=255)                                                gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                      <  1:411905/146‑1 (MQ=255)                                                gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                      >  2:411905/1‑146 (MQ=255)                                                 acacCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCt                                                                                                 <  1:754798/150‑1 (MQ=255)                                                     ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                             >  1:540111/1‑150 (MQ=255)                                                        cAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGAt                                                                                          <  2:1073156/150‑1 (MQ=255)                                                               cGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACt                                                                                   >  2:741312/1‑150 (MQ=255)                                                                 ttgtccagcCCGGCCAGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                 >  2:1512653/6‑150 (MQ=255)                                                                  gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                        >  1:338922/1‑110 (MQ=255)                                                                  gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                        <  2:338922/110‑1 (MQ=255)                                                                  gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGa                                                                                     >  1:1124798/1‑145 (MQ=255)                                                                  gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGa                                                                                     <  2:1124798/145‑1 (MQ=255)                                                                       agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                           <  1:242012/150‑1 (MQ=255)                                                                       agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGAGGTGGCGGGAATACGCCATGCGCGGCGGCCTGACCGATGCGGACTTCAGCTgg                                                                           >  1:919506/1‑150 (MQ=255)                                                                       agcCCGGCCAGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                           <  1:1512653/150‑1 (MQ=255)                                                                        gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                          <  1:678587/150‑1 (MQ=255)                                                                        gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                          <  1:1215347/150‑1 (MQ=255)                                                                                  gCTGGACAGCCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                >  1:1346836/1‑150 (MQ=255)                                                                                  gCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                <  2:1197349/150‑1 (MQ=255)                                                                                  gCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                >  2:57308/1‑150 (MQ=255)                                                                                    tGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTAAACC‑CCCTCGGCAAGCCGCGAAACGTCCGGCGGGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagc                                                              >  1:1434073/1‑150 (MQ=255)                                                                                     ggACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCTCGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                  <  2:883356/82‑2 (MQ=21)                                                                                     ggACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCTCGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                  >  1:883356/1‑81 (MQ=21)                                                                                      gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgc                                                                                                           <  2:1360518/103‑1 (MQ=255)                                                                                      gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgc                                                                                                           >  1:1360518/1‑103 (MQ=255)                                                                                           ccATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATa                                                                                                                   >  2:909626/1‑90 (MQ=255)                                                                                           ccATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATa                                                                                                                   <  1:909626/90‑1 (MQ=255)                                                                                            cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                      <  1:629370/150‑1 (MQ=255)                                                                                            cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                      >  1:464601/1‑150 (MQ=255)                                                                                            cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCACGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                      <  1:1535787/150‑1 (MQ=255)                                                                                                  gtggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAg                                                >  1:348020/1‑150 (MQ=255)                                                                                                      tgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCa                                            >  1:703446/1‑150 (MQ=255)                                                                                                          gATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                        >  1:251756/1‑150 (MQ=255)                                                                                                           aTCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGcc                                       >  2:341974/1‑150 (MQ=255)                                                                                                                 cTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                 >  1:854957/1‑150 (MQ=255)                                                                                                                  tGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                     >  1:268856/1‑129 (MQ=255)                                                                                                                  tGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                     <  2:268856/129‑1 (MQ=255)                                                                                                                   ggAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagc                                                              <  2:1588180/119‑1 (MQ=255)                                                                                                                   ggAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagc                                                              >  1:1588180/1‑119 (MQ=255)                                                                                                                   ggAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTg                               <  2:540111/150‑1 (MQ=255)                                                                                                                    gAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACATgc                              <  1:471498/150‑2 (MQ=255)                                                                                                                      aGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCa                            <  1:741312/150‑1 (MQ=255)                                                                                                                      aGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCa                            <  1:341974/150‑1 (MQ=255)                                                                                                                          gCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                        <  1:1446922/150‑1 (MQ=255)                                                                                                                           cccAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTc                       >  1:207712/1‑150 (MQ=255)                                                                                                                              aGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCaa                                     >  2:1022624/1‑133 (MQ=255)                                                                                                                              aGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCaa                                     <  1:1022624/133‑1 (MQ=255)                                                                                                                                 ccTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                 <  2:854957/150‑1 (MQ=255)                                                                                                                                  cTATGCC‑CTCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                <  2:251756/150‑1 (MQ=255)                                                                                                                                  cTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGGGCCGGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCCAGCAGCCCCACGGCAAGCTCATGGCCAATATCTTGCAGGGTCAGTTccc                <  2:919506/150‑1 (MQ=255)                                                                                                                                  cTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                <  1:237197/150‑1 (MQ=255)                                                                                                                                   tATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGCTGCCGACTTCAGCTGGGGTTCGACGGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCt               <  1:1434781/150‑1 (MQ=255)                                                                                                                                       cc‑cGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc           <  1:57308/150‑1 (MQ=255)                                                                                                                                           gCTGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGCGGCAGTTCCCTTATCGCa        >  2:1032993/1‑150 (MQ=255)                                                                                                                                                 cAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcgg  <  2:207712/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCC‑CGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGG  >  NC_002947/3791001‑3791293

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