Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,578,136 C→A 8.4% E659* (GAG→TAG) ‡ glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,578,1360CA 8.4% 233.9 / 9.3 79E659* (GAG→TAG) ‡glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base A (4/3);  total (44/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 = 3.13e-01

CCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGG  >  NC_002947/4577989‑4578278
                                                                                                                                                   |                                                                                                                                              
ccACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTcg                                                                                                                                               >  3:66270/1‑149 (MQ=255)
 cACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCTCCACCAGCACCGGCTcgc                                                                                                                                              >  5:122290/1‑149 (MQ=255)
      ccccGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctg                                                                                                                                         <  4:159207/149‑1 (MQ=255)
          gagaTTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCAGCGGCGAGCCggtgctgg                                                                                                                                     <  4:30333/149‑2 (MQ=255)
            gaTTGCCGCCATTCGGCCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCAGCGGCGAGCCggtgctggtg                                                                                                                                   <  1:371990/149‑4 (MQ=255)
             aTTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTc                                                                                                                                   <  8:254059/148‑1 (MQ=255)
              ttGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTAGCGGGGTGAAGTTGGCCTCCACCAGCACCGGCTCGCCGCTGCTGCTCCa                                                                                                                                 <  6:122290/149‑1 (MQ=255)
                 ccGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGTTCGCCGCTGCTGCTCCAgcg                                                                                                                              >  3:148153/1‑149 (MQ=255)
                 ccGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                              <  3:19540/149‑1 (MQ=255)
                 ccGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgc                                                                                                                               <  5:206000/148‑1 (MQ=255)
                 ccGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgc                                                                                                                               <  7:37614/148‑1 (MQ=255)
                  cGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                              >  5:48856/1‑148 (MQ=255)
                     cATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGc                                                                                                                          <  4:158920/149‑1 (MQ=255)
                      aTTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGcc                                                                                                                         >  8:308909/1‑149 (MQ=255)
                       ttCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctg                                                                                                                                         >  3:242223/1‑132 (MQ=255)
                       ttCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctg                                                                                                                                         <  4:242223/132‑1 (MQ=255)
                       ttCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGc                                                                                                                          >  2:58660/1‑147 (MQ=255)
                          gACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGc                                                                                                                     <  6:48856/149‑1 (MQ=255)
                            cccGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCg                                                                                                                    >  2:297066/1‑148 (MQ=255)
                                gCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCGGCAGCTCGTGCAGCGGAGGCTGGGGTacac                                                                                                               >  2:159550/1‑149 (MQ=255)
                                    acaccgcGGCCTCCCCGTTCAGCAGCTCCTGCCAACGCTCGCCGCACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTg                                                                                                           <  8:257031/143‑1 (MQ=255)
                                     cAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTg                                                                                                           <  1:58660/148‑1 (MQ=255)
                                      aGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGt                                                                                                          >  7:234820/1‑148 (MQ=255)
                                        cTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTcgccgc                                                                                                                                           <  2:62108/113‑1 (MQ=255)
                                        cTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTcgccgc                                                                                                                                           >  1:62108/1‑113 (MQ=255)
                                               tCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgtc                                                                                                <  1:297066/149‑1 (MQ=255)
                                               tCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgtc                                                                                                <  2:372225/149‑1 (MQ=255)
                                                     ttCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCAcc                                                                                          >  5:247656/1‑149 (MQ=255)
                                                       cAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGa                                                                                        >  6:251591/1‑149 (MQ=255)
                                                         gAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGAt                                                                                       >  7:111728/1‑148 (MQ=255)
                                                                   ccAACGCTCGCCGAACTGCACGCCAATGCGGTAACCCTCACGGGGAACCGGGGTTAAGTTGGCCAACACCAGCACCGGCTCGCCGCTGCTGCTACAGCACAGCAAGGCGTACACGATGTTCTGCGCGTCGTACCCGATCAACCACTgg                                                                             >  6:178457/1‑148 (MQ=255)
                                                                      aCGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAAc                                                                         <  7:308909/149‑1 (MQ=255)
                                                                       cGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAAcc                                                                        <  5:251591/149‑1 (MQ=255)
                                                                              cGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGc                                                                 >  5:200895/1‑149 (MQ=255)
                                                                               gAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCt                                                                <  7:141368/149‑1 (MQ=255)
                                                                                 aCGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTgg                                                              <  8:111728/149‑1 (MQ=255)
                                                                                    gCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctg                                                                                                                                         <  1:85610/71‑1 (MQ=255)
                                                                                    gCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctg                                                                                                                                         >  2:85610/1‑71 (MQ=255)
                                                                                      tcgacAATTCGGTAACCCTCACGTCGTACCGGGGTGAAGTTGGCCAACACCATCACCTGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGAGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATc                                                         <  5:178554/145‑1 (MQ=255)
                                                                                        gCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCAAGCGAAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCTATAAACCACTGGAAACCCTCCGGCTGGCAATc                                                         >  6:178554/1‑147 (MQ=255)
                                                                                                ggTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGACCCGGCTAGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATAAACCACTGGAAACCCTGCGGCTGGCAATCCAGCTCGTg                                                >  1:282866/1‑148 (MQ=255)
                                                                                                      ccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGcc                                         >  3:153254/1‑149 (MQ=255)
                                                                                                      ccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGcc                                         <  4:153254/149‑1 (MQ=255)
                                                                                                         cACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGc                                      >  6:239609/1‑149 (MQ=255)
                                                                                                          aCGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCgg                                       <  3:251518/147‑1 (MQ=255)
                                                                                                              ggTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAGCCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGctctac                                 >  5:410602/1‑147 (MQ=255)
                                                                                                              ggTACCGGGGTGAAGTTGGACACCATCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTTAGCCTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTc                                 >  8:244735/1‑149 (MQ=255)
                                                                                                                tACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTcgc                               >  8:4573/1‑149 (MQ=255)
                                                                                                                  ccGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTcgcg                              <  5:239609/148‑1 (MQ=255)
                                                                                                                  ccGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCgg                             >  5:391800/1‑149 (MQ=255)
                                                                                                                   cGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCgg                             >  3:359564/1‑148 (MQ=255)
                                                                                                                     gggTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTAc                          <  5:398889/149‑1 (MQ=255)
                                                                                                                        tGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAgg                       >  1:342054/1‑149 (MQ=255)
                                                                                                                            gTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGt                   >  4:236882/1‑149 (MQ=255)
                                                                                                                               ggCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGtt                  >  7:71572/1‑147 (MQ=255)
                                                                                                                                gCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCgg                    >  3:265582/1‑144 (MQ=255)
                                                                                                                                gCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCgg                    <  4:265582/144‑1 (MQ=255)
                                                                                                                                gCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg               >  8:302725/1‑149 (MQ=255)
                                                                                                                                  caccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGt             >  2:272232/1‑149 (MQ=255)
                                                                                                                                  caccaccaGCACCGACTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGt             >  4:162590/1‑149 (MQ=255)
                                                                                                                                   tggagcagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTc            >  3:193323/6‑149 (MQ=255)
                                                                                                                                   tggagcagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTc            >  3:233987/6‑149 (MQ=255)
                                                                                                                                   accaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTc            >  1:155790/1‑149 (MQ=255)
                                                                                                                                    ccaccaGCACCGGTTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCGGGCTCTTCGCGGTACAGGCGGTTGAGGTCg           <  4:148153/149‑1 (MQ=255)
                                                                                                                                    ccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCg           >  4:298402/1‑149 (MQ=255)
                                                                                                                                       gcagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCt                                                                <  7:165227/92‑1 (MQ=255)
                                                                                                                                       gcagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCt                                                                >  8:165227/2‑93 (MQ=255)
                                                                                                                                       ccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGccc        <  4:359564/149‑1 (MQ=255)
                                                                                                                                       ccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGcc         >  1:321233/1‑148 (MQ=255)
                                                                                                                                        caGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCcca       <  2:282866/149‑1 (MQ=255)
                                                                                                                                        caGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCcca       >  8:202729/1‑149 (MQ=255)
                                                                                                                                          gCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                               <  2:122409/91‑1 (MQ=255)
                                                                                                                                          gCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                               >  1:122409/1‑91 (MQ=255)
                                                                                                                                             ccGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAgg  >  6:196889/1‑149 (MQ=255)
                                                                                                                                             ccGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAgg  >  5:279628/1‑149 (MQ=255)
                                                                                                                                              cGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAACCCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGa                >  5:95438/1‑134 (MQ=255)
                                                                                                                                              cGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGa                <  6:95438/134‑1 (MQ=255)
                                                                                                                                                gCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTcg                                >  6:47599/1‑116 (MQ=255)
                                                                                                                                                gCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTcg                                <  5:47599/116‑1 (MQ=255)
                                                                                                                                                  tCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTgg                                                              >  7:187224/1‑84 (MQ=255)
                                                                                                                                                  tCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTgg                                                              >  7:181913/1‑84 (MQ=255)
                                                                                                                                                  tCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTgg                                                              <  8:181913/84‑1 (MQ=255)
                                                                                                                                                  tCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTgg                                                              <  8:187224/84‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGG  >  NC_002947/4577989‑4578278

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: