Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,575,045 C→T 7.2% L573L (CTG→TTG)  glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,575,0450CT 7.2% 202.7 / 4.1 69L573L (CTG→TTG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/3);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

GCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGCGT  >  NC_002947/4574898‑4575187
                                                                                                                                                    |                                                                                                                                              
gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgctgct                                                                                                                                               <  6:299900/149‑1 (MQ=255)
gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAgctgct                                                                                                                                               <  2:159713/149‑1 (MQ=255)
    ccTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGcc                                                                                                                                           <  3:113736/149‑1 (MQ=255)
        gctgctGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCtt                                                                                                                                         <  4:178995/147‑1 (MQ=255)
        gctgctGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTc                                                                                                                                        <  6:357658/148‑1 (MQ=255)
         ctgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAGCAGCTCTTGGCCTTcc                                                                                                                                       <  6:352553/148‑2 (MQ=255)
          tgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGTCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAGCAGCCCTTGGCCTTCCTc                                                                                                                                     <  6:282180/149‑1 (MQ=255)
            ctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGACGAAGGCCAGCAGCTCTTGGCCTTCCTc                                                                                                                                     <  5:374117/147‑1 (MQ=255)
               gCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCATCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAgc                                                                                                                                 <  5:84381/148‑1 (MQ=255)
                 ccAGGGCACGCCAATGATTGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGCAAGGCCAAGAGCTGCTGGACTTCGTCAAGCGcc                                                                                                                              >  8:271353/1‑149 (MQ=255)
                     ggCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGAc                                                                                                                          >  1:24608/1‑149 (MQ=255)
                     ggCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGAc                                                                                                                          >  2:190918/1‑149 (MQ=255)
                          gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcc                                                                                                                     >  6:172490/1‑149 (MQ=255)
                          gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcc                                                                                                                     >  6:237446/1‑149 (MQ=255)
                                 aTCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCAGGGTGATCTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcctgcgcct                                                                                                              >  6:237508/1‑149 (MQ=255)
                                  tCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGcgcctgcgcctg                                                                                                             <  1:190918/149‑1 (MQ=255)
                                             gacgaGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGt                                                                                                  >  4:294554/1‑149 (MQ=255)
                                                gaGTTCAGCCGCACCCAGCCCGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGtgct                                                                                               <  6:250248/149‑1 (MQ=255)
                                                   ttCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTgcg                                                                                            >  7:150783/1‑149 (MQ=255)
                                                   ttCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTgcg                                                                                            >  6:276530/1‑149 (MQ=255)
                                                      aGCCGCACCCAGCAAGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGcgccg                                                                                         >  2:65706/1‑149 (MQ=255)
                                                        ccGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGGTCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCt                                                                                       >  2:296826/1‑149 (MQ=255)
                                                              ccAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCt                                                                                 >  7:182047/1‑149 (MQ=255)
                                                               cAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCgg                                                                                                   >  3:153656/1‑130 (MQ=255)
                                                               cAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCgg                                                                                                   <  4:153656/130‑1 (MQ=255)
                                                                aGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCttt                                                                               <  8:182047/149‑1 (MQ=255)
                                                                   actgCAACAAGAACGCTTACTGCCAGGAGAGCGAGAGCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTGGTCAAGCGCCTGACGCGCCTGCGCCTGGCGTACCCGGTGCTGCGCCGCTCGCGCTTTCtg                                                                            <  2:303800/146‑1 (MQ=255)
                                                                   aCGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtg                                                                            <  5:71762/149‑1 (MQ=255)
                                                                      gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                         <  3:158071/149‑1 (MQ=255)
                                                                      gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                         <  2:67109/149‑1 (MQ=255)
                                                                      gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                         <  2:61719/149‑1 (MQ=255)
                                                                       caacaacaaCGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtgg                                                                        <  8:150783/149‑1 (MQ=255)
                                                                         acaacaacTCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGAGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                      <  8:83252/149‑1 (MQ=255)
                                                                         acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                      <  5:172490/149‑1 (MQ=255)
                                                                         acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                      <  5:237446/149‑1 (MQ=255)
                                                                         acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTCGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                      <  5:359976/149‑1 (MQ=255)
                                                                          caacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                      >  6:76630/1‑148 (MQ=255)
                                                                           aacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGa                                                                    >  7:26388/1‑149 (MQ=255)
                                                                             caacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCGTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGAc                                                                   >  1:84042/1‑149 (MQ=255)
                                                                             caacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGAc                                                                   >  3:156271/1‑148 (MQ=255)
                                                                                cGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTAc                                                                >  5:19089/1‑148 (MQ=255)
                                                                                     aCTGCCAGGACAGCGAGATCGTCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                           <  3:197603/148‑1 (MQ=255)
                                                                                      cTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAg                                                          <  2:36615/148‑1 (MQ=255)
                                                                                         ccAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTACGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCaa                                                      >  7:189266/1‑149 (MQ=255)
                                                                                            ggACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTg                                                   <  4:156271/149‑1 (MQ=255)
                                                                                            ggACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTg                                                   <  3:294554/149‑1 (MQ=255)
                                                                                            ggACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGAGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCAGCTCGCGCTTTCTGGTGGGTGACTACAATTAGGCAATTg                                                   >  4:106155/1‑149 (MQ=255)
                                                                                             gACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTgg                                                  <  8:26388/149‑1 (MQ=255)
                                                                                              aCAGCGAGATCGGCGGGGTGAACGGGGACCTGGACGAGGACGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTggg                                                 <  6:418/149‑1 (MQ=255)
                                                                                                  cGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTc                                              >  7:5273/1‑148 (MQ=255)
                                                                                                    agaTCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAg                                           <  5:76630/149‑1 (MQ=255)
                                                                                                         ggCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGt                                      <  3:2066/149‑1 (MQ=255)
                                                                                                          gCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGt                                      >  2:239713/1‑148 (MQ=255)
                                                                                                           cTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                           >  6:146785/1‑126 (MQ=255)
                                                                                                           cTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                           <  5:146785/126‑1 (MQ=255)
                                                                                                            tGGGTGAACTGGGACCTGGACGAGGAAGGCCACGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGCGTACATTGAGGCTATTGGGGGCAAGGAGGTgcc                                   >  6:42584/1‑147 (MQ=255)
                                                                                                            tGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCGTCGTCAAGCGCCTGACGCGCCGGCGCCTGGCCTCCCCGGTGCTGCGCCGCTCGCGCGTTGGGGGGGGTGACTGCAAGGAGGCAATGGGGGGCAAGGACGTGa                                    >  6:37176/1‑148 (MQ=255)
                                                                                                              ggTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGAc                                   <  1:296826/147‑1 (MQ=255)
                                                                                                                   aCTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctggctg                            <  3:84426/149‑1 (MQ=255)
                                                                                                                    cTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctggctgg                           <  2:84042/149‑1 (MQ=255)
                                                                                                                           cTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGa                    <  5:276530/149‑1 (MQ=255)
                                                                                                                               gACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATg                  <  7:271353/148‑1 (MQ=255)
                                                                                                                                  gAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAg              <  2:24608/149‑1 (MQ=255)
                                                                                                                                  gAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAg              <  6:19089/149‑1 (MQ=255)
                                                                                                                                   acgAAGGCCAGCAGCTCTTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGc             >  2:81915/3‑149 (MQ=255)
                                                                                                                                      aaGGCCAGCAGCTCTTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGc             >  1:5739/1‑146 (MQ=255)
                                                                                                                                         gCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATg       <  1:65706/149‑1 (MQ=255)
                                                                                                                                           cAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAg     >  8:186182/1‑149 (MQ=255)
                                                                                                                                           cAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAg     >  4:355190/1‑149 (MQ=255)
                                                                                                                                            aaGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGc    >  6:297437/1‑149 (MQ=255)
                                                                                                                                              gagCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGAGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGCGt  >  4:320647/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGCGT  >  NC_002947/4574898‑4575187

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: