Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,574,229	G→C	5.9%	G301R (GGC→CGC)	glgX →	glycogen debranching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,574,229	0	G	C	5.9%	218.8 / 3.5	68	G301R (GGC→CGC)	glgX	glycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

ACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGC  >  NC_002947/4574083‑4574372                                                                                                                                                    |                                                                                                                                                 aCCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCgg                                                                                                                                                 >  7:174744/1‑148 (MQ=255)        cTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACa                                                                                                                                         <  6:151719/149‑1 (MQ=255)                      gCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGtactac                                                                                                                          <  5:425312/149‑1 (MQ=255)                           aTCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCt                                                                                                                     >  2:28997/1‑149 (MQ=255)                            tCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCt                                                                                                                     >  6:418411/1‑148 (MQ=255)                            tCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCt                                                                                                                     <  2:341059/148‑1 (MQ=255)                              gCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGa                                                                                                                   <  8:174744/148‑1 (MQ=255)                              gCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGAt                                                                                                                  <  5:183883/149‑1 (MQ=255)                               ccGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATg                                                                                                                 <  7:345254/149‑1 (MQ=255)                               ccGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATg                                                                                                                 >  8:194290/1‑149 (MQ=255)                                 gAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGcc                                                                                                               >  5:57001/1‑149 (MQ=255)                                 gAGTTCAAGGAGATGGGCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGcc                                                                                                               <  5:37520/149‑1 (MQ=255)                                    ttCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGga                                                                                                            <  5:418411/149‑1 (MQ=255)                                       aaGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGc                                                                                                                                     >  8:310985/1‑122 (MQ=255)                                       aaGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGc                                                                                                                                     <  7:310985/122‑1 (MQ=255)                                       aaGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGgacga                                                                                                         >  5:375060/1‑149 (MQ=255)                                       aaGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACagg                                                                                                                                       <  2:315073/120‑3 (MQ=255)                                       aaGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACagg                                                                                                                                       >  1:315073/1‑118 (MQ=255)                                        aGGAGATGGTCGCGCACCTGCACGACGCGGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCTgacga                                                                                                         >  1:232985/1‑149 (MQ=255)                                        aGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGCCACTGTCGATGCGAGGCATCGACAACGCCT‑CGTACTACCGCCAGCTGCCggtgccg                                                                                                        >  4:400692/1‑144 (MQ=255)                                        aGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCTgacgac                                                                                                        <  2:232985/149‑1 (MQ=255)                                        aGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGgacgac                                                                                                        <  8:189752/149‑1 (MQ=255)                                        aGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGAAGCCGgacga                                                                                                         >  7:189752/1‑148 (MQ=255)                                         ggAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGgacgac                                                                                                        >  6:54037/1‑148 (MQ=255)                                         ggAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACa                                                                                                       <  6:375060/149‑1 (MQ=255)                                          gagaTGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACaa                                                                                                      >  6:301784/1‑149 (MQ=255)                                               ggTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAgcgc                                                                                                  >  7:27883/1‑148 (MQ=255)                                                gTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGcta                                                                                                >  5:280564/1‑149 (MQ=255)                                                      cACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACAt                                                                                          <  6:280564/149‑1 (MQ=255)                                                           gCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCt                                                                                                                                   <  5:171299/104‑1 (MQ=255)                                                           gCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCt                                                                                                                                   >  6:171299/1‑104 (MQ=255)                                                            cACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCAATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGa                                                                                    >  1:390846/1‑149 (MQ=255)                                                             acgacgC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                               >  2:395387/1‑137 (MQ=255)                                                             acgacgC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                               <  1:395387/137‑1 (MQ=255)                                                              cgacgC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACt                                                                                  <  6:57001/149‑1 (MQ=255)                                                                acgC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTcc                                                                                <  1:28997/149‑1 (MQ=255)                                                                acgC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTcc                                                                                <  7:194290/149‑1 (MQ=255)                                                                              ggTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAaca                                                                   >  7:16999/1‑149 (MQ=255)                                                                               gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCCGTACTACCGCCTGATGCCGGACGACAAGCGct                                                                                                 <  1:55200/119‑1 (MQ=255)                                                                               gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGct                                                                                                 >  2:55200/1‑118 (MQ=255)                                                                               gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAacac                                                                  >  2:1792/1‑149 (MQ=255)                                                                                         aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACAcc                                                                 <  4:363849/140‑1 (MQ=255)                                                                                         aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACAcc                                                                 >  3:363849/1‑140 (MQ=255)                                                                                          cGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                       <  5:54037/149‑1 (MQ=255)                                                                                          cGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                       <  5:301784/149‑1 (MQ=255)                                                                                           gTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                       <  7:261360/148‑1 (MQ=255)                                                                                            tGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCAATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGc                                                     <  2:390846/149‑1 (MQ=255)                                                                                               tCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGAc                                                           >  7:344681/1‑140 (MQ=255)                                                                                               tCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGAc                                                           <  8:344681/140‑1 (MQ=255)                                                                                                cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAcc                                                 >  4:46149/1‑149 (MQ=255)                                                                                                 tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                >  2:1297/1‑149 (MQ=255)                                                                                                   cAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCt                                              <  1:1792/149‑1 (MQ=255)                                                                                                    aaCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccggc                                                                       <  5:64868/123‑1 (MQ=255)                                                                                                    aaCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccggc                                                                       >  6:64868/1‑123 (MQ=255)                                                                                                          acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgct                                       >  1:60630/1‑149 (MQ=255)                                                                                                           ccgccgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                      <  1:1297/149‑1 (MQ=255)                                                                                                             gccgAAGGGAACGAGGGCGGCCCGACGCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCa                                    <  8:27883/149‑1 (MQ=255)                                                                                                                      aaCGAGCGCTGCCCGACCCTGTCGATGCGCGGAATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTTCAGTTGGTCa                            <  1:258881/148‑1 (MQ=255)                                                                                                                      aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCCACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAc                           >  3:246644/1‑149 (MQ=255)                                                                                                                      aaCGAGCGCGGCCCGACCCTGTCGATGAGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACAGGAAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCa                            >  2:258881/1‑148 (MQ=255)                                                                                                                            cgcgGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtc                     >  6:285350/1‑149 (MQ=255)                                                                                                                              cgGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgc                   >  8:60942/1‑149 (MQ=255)                                                                                                                                   cGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgc              >  3:120344/1‑149 (MQ=255)                                                                                                                                       gttGTCGATGCCGCGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACt          >  1:221609/3‑149 (MQ=255)                                                                                                                                       ccTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACt          <  7:55690/149‑1 (MQ=255)                                                                                                                                           tCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg        <  3:183955/147‑1 (MQ=255)                                                                                                                                            cGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCa                            <  8:321633/126‑1 (MQ=255)                                                                                                                                            cGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCa                            >  7:321633/1‑126 (MQ=255)                                                                                                                                               tGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGc  <  2:60630/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGC‑GGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCT‑CGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGC  >  NC_002947/4574083‑4574372

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