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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,575,0380AG 7.2% 197.2 / 4.8 70Q570Q (CAA→CAG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (27/37);  new base G (2/3);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGC  >  NC_002947/4574896‑4575185
                                                                                                                                               |                                                                                                                                                   
tCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAgctg                                                                                                                                               >  4:2066/1‑149 (MQ=255)
  gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAgctgct                                                                                                                                             <  6:299900/149‑1 (MQ=255)
  gCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAgctgct                                                                                                                                             <  2:159713/149‑1 (MQ=255)
      ccTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGcc                                                                                                                                         <  3:113736/149‑1 (MQ=255)
          gctgctGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCtt                                                                                                                                       <  4:178995/147‑1 (MQ=255)
          gctgctGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTc                                                                                                                                      <  6:357658/148‑1 (MQ=255)
           ctgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAGCAGCTCTTGGCCTTcc                                                                                                                                     <  6:352553/148‑2 (MQ=255)
            tgctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGTCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAGCAGCCCTTGGCCTTCCTc                                                                                                                                   <  6:282180/149‑1 (MQ=255)
              ctgGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGACGAAGGCCAGCAGCTCTTGGCCTTCCTc                                                                                                                                   <  5:374117/147‑1 (MQ=255)
                 gCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCATCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAgc                                                                                                                               <  5:84381/148‑1 (MQ=255)
                   ccAGGGCACGCCAATGATTGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGCAAGGCCAAGAGCTGCTGGACTTCGTCAAGCGcc                                                                                                                            >  8:271353/1‑149 (MQ=255)
                       ggCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGAc                                                                                                                        >  2:190918/1‑149 (MQ=255)
                       ggCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGAc                                                                                                                        >  1:24608/1‑149 (MQ=255)
                            gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcc                                                                                                                   >  6:172490/1‑149 (MQ=255)
                            gCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcc                                                                                                                   >  6:237446/1‑149 (MQ=255)
                                   aTCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCAGGGTGATCTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACgcgcctgcgcct                                                                                                            >  6:237508/1‑149 (MQ=255)
                                    tCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGcgcctgcgcctg                                                                                                           <  1:190918/149‑1 (MQ=255)
                                               gacgaGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGt                                                                                                >  4:294554/1‑149 (MQ=255)
                                                  gaGTTCAGCCGCACCCAGCCCGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGtgct                                                                                             <  6:250248/149‑1 (MQ=255)
                                                     ttCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTgcg                                                                                          >  7:150783/1‑149 (MQ=255)
                                                     ttCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTgcg                                                                                          >  6:276530/1‑149 (MQ=255)
                                                        aGCCGCACCCAGCAAGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGcgccg                                                                                       >  2:65706/1‑149 (MQ=255)
                                                          ccGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGGTCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCt                                                                                     >  2:296826/1‑149 (MQ=255)
                                                                ccAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCt                                                                               >  7:182047/1‑149 (MQ=255)
                                                                 cAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCgg                                                                                                 >  3:153656/1‑130 (MQ=255)
                                                                 cAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCgg                                                                                                 <  4:153656/130‑1 (MQ=255)
                                                                  aGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCttt                                                                             <  8:182047/149‑1 (MQ=255)
                                                                     actgCAACAAGAACGCTTACTGCCAGGAGAGCGAGAGCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTGGTCAAGCGCCTGACGCGCCTGCGCCTGGCGTACCCGGTGCTGCGCCGCTCGCGCTTTCtg                                                                          <  2:303800/146‑1 (MQ=255)
                                                                     aCGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtg                                                                          <  5:71762/149‑1 (MQ=255)
                                                                        gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                       <  3:158071/149‑1 (MQ=255)
                                                                        gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                       <  2:67109/149‑1 (MQ=255)
                                                                        gCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtg                                                                       <  2:61719/149‑1 (MQ=255)
                                                                         caacaacaaCGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCtggtgg                                                                      <  8:150783/149‑1 (MQ=255)
                                                                           acaacaacTCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGAGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                    <  8:83252/149‑1 (MQ=255)
                                                                           acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                    <  5:237446/149‑1 (MQ=255)
                                                                           acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                    <  5:172490/149‑1 (MQ=255)
                                                                           acaacaacGCTTACTGCCAGGACAGCGAGATCGGCTCGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                    <  5:359976/149‑1 (MQ=255)
                                                                            caacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTggtgggt                                                                    >  6:76630/1‑148 (MQ=255)
                                                                             aacaacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGa                                                                  >  7:26388/1‑149 (MQ=255)
                                                                               caacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACCGTGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGAc                                                                 >  1:84042/1‑149 (MQ=255)
                                                                               caacGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGAc                                                                 >  3:156271/1‑148 (MQ=255)
                                                                                  cGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTAc                                                              >  5:19089/1‑148 (MQ=255)
                                                                                       aCTGCCAGGACAGCGAGATCGTCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                         <  3:197603/148‑1 (MQ=255)
                                                                                        cTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAg                                                        <  2:36615/148‑1 (MQ=255)
                                                                                           ccAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTACGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCaa                                                    >  7:189266/1‑149 (MQ=255)
                                                                                              ggACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTg                                                 <  4:156271/149‑1 (MQ=255)
                                                                                              ggACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTg                                                 <  3:294554/149‑1 (MQ=255)
                                                                                              ggACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGAGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCAGCTCGCGCTTTCTGGTGGGTGACTACAATTAGGCAATTg                                                 >  4:106155/1‑149 (MQ=255)
                                                                                               gACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTgg                                                <  8:26388/149‑1 (MQ=255)
                                                                                                aCAGCGAGATCGGCGGGGTGAACGGGGACC‑TGGACGAGGACGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTggg                                               <  6:418/149‑1 (MQ=255)
                                                                                                    cGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTc                                            >  7:5273/1‑148 (MQ=255)
                                                                                                      agaTCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAg                                         <  5:76630/149‑1 (MQ=255)
                                                                                                           ggCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGt                                    <  3:2066/149‑1 (MQ=255)
                                                                                                            gCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGt                                    >  2:239713/1‑148 (MQ=255)
                                                                                                             cTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                         >  6:146785/1‑126 (MQ=255)
                                                                                                             cTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGa                                                         <  5:146785/126‑1 (MQ=255)
                                                                                                              tGGGTGAACTGGGACC‑TGGACGAGGAAGGCCACGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGCGTACATTGAGGCTATTGGGGGCAAGGAGGTgcc                                 >  6:42584/1‑147 (MQ=255)
                                                                                                              tGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCGTCGTCAAGCGCCTGACGCGCCGGCGCCTGGCCTCCCCGGTGCTGCGCCGCTCGCGCGTTGGGGGGGGTGACTGCAAGGAGGCAATGGGGGGCAAGGACGTGa                                  >  6:37176/1‑148 (MQ=255)
                                                                                                                ggTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGAc                                 <  1:296826/147‑1 (MQ=255)
                                                                                                                     aCTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctggctg                          <  3:84426/149‑1 (MQ=255)
                                                                                                                      cTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACctggctgg                         <  2:84042/149‑1 (MQ=255)
                                                                                                                             c‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGa                  <  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)
                                                                                                                                               |                                                                                                                                                   
TCGCCACCTTGCTGCTGGCCCAGGGCACGCCAATGATCGTTGCCGGCGACGAGTTCAGCCGCACCCAGCACGGCAACAACAACGCTTACTGCCAGGACAGCGAGATCGGCTGGGTGAACTGGGACC‑TGGACGAGGAAGGCCAAGAGCTGCTGGCCTTCGTCAAGCGCCTGACGCGCCTGCGCCTGGCCTACCCGGTGCTGCGCCGCTCGCGCTTTCTGGTGGGTGACTACAATGAGGCAATTGGGGTCAAGGACGTGACCTGGCTGGCACCCGATGGCAGCGAGATGAGC  >  NC_002947/4574896‑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: ‑