Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,575,038	A→G	6.8%	Q570Q (CAA→CAG)	glgX →	glycogen debranching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,575,038	0	A	G	6.8%	169.8 / 2.6	59	Q570Q (CAA→CAG)	glgX	glycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TACTTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGC  >  NC_002947/4574892‑4575185                                                                                                                                                   |                                                                                                                                                     tactTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAaga                                                                                                                                                    <  5:52254/149‑1 (MQ=255)  actTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAagag                                                                                                                                                   >  5:34319/1‑149 (MQ=255)   ctTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAagacc                                                                                                                                                  >  8:254089/1‑147 (MQ=255)   ctTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgc                                                                                                                                                  >  6:211461/1‑149 (MQ=255)    ttCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgct                                                                                                                                                 >  4:311708/1‑149 (MQ=255)       gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgctgc                                                                                                                                               <  5:296103/148‑1 (MQ=255)       gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgctgc                                                                                                                                               <  1:309257/148‑1 (MQ=255)        ccACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGctgctg                                                                                                                                             <  1:219134/149‑1 (MQ=255)        ccACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGctgctg                                                                                                                                             <  5:301831/149‑1 (MQ=255)          aCCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGc                                                                                                                                           >  4:168505/1‑149 (MQ=255)             ttGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCtt                                                                                                                                        >  8:66000/1‑149 (MQ=255)               gctgctGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTc                                                                                                                                       <  2:21372/148‑1 (MQ=255)                ctgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCATCCGGACCCAGCACGGCAACAACAACGCTTACTGCCAGGACATCGAGATCGGCTGGGGGAACTGGGAGCTGGACGAGGAAGGCCAGCAGCTCTTGGCCTTcc                                                                                                                                      <  4:221473/148‑2 (MQ=255)                ctgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAGCAGCTCTTGGCCTTcc                                                                                                                                      <  5:118633/148‑2 (MQ=255)                     ggCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAgc                                                                                                                                <  1:315462/149‑1 (MQ=255)                            ggCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGTCAACAACATCGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑c                                                                                                                        <  7:8508/149‑1 (MQ=255)                               aCGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑cgcg                                                                                                                     >  7:74544/1‑149 (MQ=255)                                 gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcc                                                                                                                   >  6:246974/1‑149 (MQ=255)                                 gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcc                                                                                                                   >  2:289032/1‑149 (MQ=255)                                  ccAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcct                                                                                                                  <  3:311708/149‑1 (MQ=255)                                  ccAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcct                                                                                                                  >  1:57198/1‑149 (MQ=255)                                   cAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcctg                                                                                                                 <  8:116094/149‑1 (MQ=255)                                           gTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTgg                                                                                                          <  2:32317/148‑1 (MQ=255)                                           gTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGc                                                                                                         <  3:168505/149‑1 (MQ=255)                                            ttGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcctgcgc                                                                                                              >  1:243802/1‑143 (MQ=255)                                            ttGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑Cgcgcctgcgc                                                                                                              <  2:243802/143‑1 (MQ=255)                                                cGGCGACGACTTCTGCCGCAGCCTGCACGGCAACAACTTCGCTTAGTGCGAGGACAGCGAGTTCGGCTGGGTGAAGTGGGAGGTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGGCTGCGCCTGGCCTAcc                                                                                                    <  4:169741/149‑1 (MQ=255)                                                     acgaGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGtg                                                                                               <  7:66000/149‑1 (MQ=255)                                                      cgaGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGtgc                                                                                              >  6:316860/1‑149 (MQ=255)                                                              gCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGcgccgc                                                                                      >  3:145985/1‑149 (MQ=255)                                                                 gCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTcg                                                                                   <  2:57198/149‑1 (MQ=255)                                                                  cACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTcgc                                                                                  <  5:316860/149‑1 (MQ=255)                                                                            ggCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggt                                                                        <  8:74544/149‑1 (MQ=255)                                                                            ggCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggt                                                                        >  8:77888/1‑149 (MQ=255)                                                                              caacaacaaCGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtgg                                                                      >  4:42291/1‑149 (MQ=255)                                                                                  aacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGa                                                                  >  4:171801/1‑149 (MQ=255)                                                                                     aacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTa                                                               <  3:171801/149‑1 (MQ=255)                                                                                         cTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGa                                                                  <  7:98348/142‑1 (MQ=255)                                                                                         cTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGa                                                                  >  8:98348/1‑142 (MQ=255)                                                                                         cTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAAt                                                           <  5:246974/149‑1 (MQ=255)                                                                                            aCTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATg                                                          <  1:289032/147‑1 (MQ=255)                                                                                            aCTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATg                                                          <  7:254089/147‑1 (MQ=255)                                                                                            aCTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                         <  7:77888/148‑1 (MQ=255)                                                                                             cTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAg                                                        <  2:117423/148‑1 (MQ=255)                                                                                              tGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAgg                                                       >  6:93163/1‑148 (MQ=255)                                                                                                  aGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACTAGGAAGGCCAAGAGCTGCTGGCCTTGGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAAtt                                                  <  6:158767/149‑1 (MQ=255)                                                                                                        gCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGt                                             >  8:46388/1‑148 (MQ=255)                                                                                                        gCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGt                                             >  8:41840/1‑148 (MQ=255)                                                                                                         cGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGG‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTc                                            >  1:42766/1‑148 (MQ=255)                                                                                                         cGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTc                                            >  4:297036/1‑148 (MQ=255)                                                                                                                 gCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAg                                                        <  7:40136/128‑1 (MQ=255)                                                                                                                 gCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAg                                                        >  8:40136/1‑128 (MQ=255)                                                                                                                     ggTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGAc                                 <  4:145985/147‑1 (MQ=255)                                                                                                                      gTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctg                              >  3:149387/1‑149 (MQ=255)                                                                                                                         aaCTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctggct                           >  6:266336/1‑149 (MQ=255)                                                                                                                                    ggACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGAGGGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGAt                 >  1:230869/1‑149 (MQ=255)                                                                                                                                        gacgAAGGCCAGCAGCTCTTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAg            >  3:19080/4‑149 (MQ=255)                                                                                                                                          cgAAGGCCAGCAGCTCTTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGc           >  2:310226/2‑148 (MQ=255)                                                                                                                                           gAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCga         <  5:211461/149‑1 (MQ=255)                                                                                                                                            aaGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGG‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGCCATCGGGTGccagccag        <  2:42766/149‑1 (MQ=255)                                                                                                                                             aGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCgaga       <  3:42291/149‑1 (MQ=255)                                                                                                                                             aGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGACACCCGATGGCAGc           >  5:178343/1‑145 (MQ=255)                                                                                                                                                  aaGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGc  <  6:34319/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TACTTCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGA‑CGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGC  >  NC_002947/4574892‑4575185

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