Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,435,225 A→T 6.2% L249Q (CTG→CAG)  PP_1256 ← alpha‑ketoglutarate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,435,2250AT 6.2% 297.9 / 8.1 96L249Q (CTG→CAG) PP_1256alpha‑ketoglutarate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base T (3/3);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CGGCGCCCTGGGTCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTGC  >  NC_002947/1435076‑1435368
                                                                                                                                                     |                                                                                                                                                
cGGCGCCCTGGGTCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGcaca                                                                                                                                                  >  2:1868177/1‑150 (MQ=255)
           gTCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgccatgcca                                                                                                                                       <  1:1671505/150‑1 (MQ=255)
            tCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCAc                                                                                                                                      >  2:363538/1‑150 (MQ=255)
            tCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGATGATCGAACGCATAACGGCATACACAGGGATGGGCTCCGGGCGCGCCTGGACGGCTAGGCACAGttccatgccgc                                                                                                                                      >  2:451389/1‑148 (MQ=255)
              aGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCg                                                                                                                                    >  1:2152353/1‑150 (MQ=255)
               gCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGc                                                                                                                                   >  2:647246/1‑150 (MQ=255)
               gCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGc                                                                                                                                   >  1:2251053/1‑150 (MQ=255)
               gCGAGGCGAAGAACACTTGTCACCGTGCTCAGGCTCGGGCTTGCAGGGGGGGCCCTAACAGGAACACCGGGTTGAGCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGc                                                                                                                                   >  2:1233812/1‑150 (MQ=255)
                cGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGTGGCATGGCACTGtgccaagccgcccag                                                                                                                                  <  1:458046/150‑7 (MQ=255)
                cGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCt                                                                                                                                  >  2:1439926/1‑150 (MQ=255)
                         gaAACCTCGTGCCAGTGCTTCGGCACGGGCTTGAAGGGCGGCGTCGAACAGGAAAACCGGGTTGCTAGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGAACGCCTGGGCCGCTTGGCACAGTCCCATGCCACCGCCGCGCGAGcc                                                                                                                         >  2:650042/1‑150 (MQ=255)
                          aaaCCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCAAAGTGCCATGCCACCGCTGCGCGAGCCg                                                                                                                        <  2:1248744/150‑1 (MQ=255)
                           aaCCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCgg                                                                                                                       <  1:1439926/150‑1 (MQ=255)
                            aCCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGt                                                                                                                      <  1:647246/150‑1 (MQ=255)
                              cTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGa                                                                                                                    >  2:1431679/1‑150 (MQ=255)
                               ttGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGaa                                                                                                                   >  2:742590/1‑150 (MQ=255)
                                    ccAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCaa                                                                                                              <  2:2152353/150‑1 (MQ=255)
                                    ccAGTGCTTCGCCACGCGCTTGCAGCGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCaa                                                                                                              <  2:1157644/150‑1 (MQ=255)
                                     cAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGGGAAACCAAc                                                                                                             >  1:1701554/1‑150 (MQ=255)
                                           ttCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCTATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCtt                                                                                                       <  2:1395540/150‑1 (MQ=255)
                                           ttCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCtt                                                                                                       <  1:1868177/150‑1 (MQ=255)
                                                  cGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCACCGGCCTTGAtgcat                                                                                                <  2:1299405/150‑5 (MQ=255)
                                                  cGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGt                                                                                                <  2:1830922/150‑1 (MQ=255)
                                                  cGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGt                                                                                                <  2:1444512/150‑1 (MQ=255)
                                                  cGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGt                                                                                                <  1:2173257/150‑1 (MQ=255)
                                                   gggCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgcca                                                                                                                                            <  2:2279015/105‑1 (MQ=255)
                                                   gggCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgcca                                                                                                                                            >  1:2279015/1‑105 (MQ=255)
                                                       ttGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTc                                                                                          >  1:1471821/1‑150 (MQ=255)
                                                           aGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTg                                                                                      >  2:1375135/1‑150 (MQ=255)
                                                             ggcggcGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCt                                                                                       <  2:27446/147‑1 (MQ=255)
                                                             ggcggcGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCt                                                                                       >  1:27446/1‑147 (MQ=255)
                                                               cggcgTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCa                                                                                  <  1:1375135/150‑1 (MQ=255)
                                                                 gcgTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCaga                                                                                <  1:363538/150‑1 (MQ=255)
                                                                 gcgTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCaga                                                                                <  2:2087834/150‑1 (MQ=255)
                                                                      gAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgccat                                                                                                                                           <  1:616396/87‑1 (MQ=255)
                                                                      gAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgccat                                                                                                                                           >  2:616396/1‑87 (MQ=255)
                                                                      gAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCATGGCACTGtgccaagc                                                                                                                                         >  2:1293454/1‑86 (MQ=255)
                                                                      gAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCATGGCACTGtgccaagc                                                                                                                                         <  1:1293454/89‑1 (MQ=255)
                                                                      acaCAGTAACACCGGGTTGCTCGAACTCTCTTCGGCATACATCGGGATTGGCTCGGGTCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGGTGCGGGAGCCGTTGAAACCATCGGCATTGATACGTGGGTCTCTGTCCAGAGCAAt                                                                           <  2:505892/148‑1 (MQ=255)
                                                                       aaCAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATg                                                                          >  1:2395412/1‑150 (MQ=255)
                                                                       aaCAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATg                                                                          >  2:886078/1‑150 (MQ=255)
                                                                        aCAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGTCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGc                                                                         >  1:2301130/1‑150 (MQ=255)
                                                                        aCAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGc                                                                         >  2:1696563/1‑150 (MQ=255)
                                                                        aCAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGGGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGc                                                                         >  2:1312693/1‑150 (MQ=255)
                                                                            gAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAAt                                                                           >  1:29521/1‑144 (MQ=255)
                                                                            gAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAAt                                                                           <  2:29521/144‑1 (MQ=255)
                                                                            gAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCAc                                                                     >  1:1377697/1‑150 (MQ=255)
                                                                                acCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCaga                                                                                <  1:77896/135‑1 (MQ=255)
                                                                                acCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCaga                                                                                >  2:77896/1‑135 (MQ=255)
                                                                                    ggTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCAc                                                                     <  2:1372419/142‑1 (MQ=255)
                                                                                    ggTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCAc                                                                     >  1:1372419/1‑142 (MQ=255)
                                                                                        gATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCGGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGa                                                         <  2:349148/150‑1 (MQ=255)
                                                                                             aaCTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgcc                                                                                                                                             <  1:499596/62‑1 (MQ=255)
                                                                                             aaCTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGtgcc                                                                                                                                             >  2:499596/1‑62 (MQ=255)
                                                                                                  aTTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCg                                                                                                                        <  1:1038961/78‑1 (MQ=255)
                                                                                                  aTTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCg                                                                                                                        >  2:1038961/1‑78 (MQ=255)
                                                                                                  aTTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAagc                                               >  1:2259338/1‑150 (MQ=255)
                                                                                                     tCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAagcagc                                            >  1:889377/1‑150 (MQ=255)
                                                                                                         cATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCgaga                                        >  2:1264085/1‑150 (MQ=255)
                                                                                                          aTACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCATAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGaa                                       <  1:451389/150‑1 (MQ=255)
                                                                                                          aTACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGaa                                       >  2:2158642/1‑150 (MQ=255)
                                                                                                              acCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCAt                                                    >  1:1408978/1‑133 (MQ=255)
                                                                                                              acCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCAt                                                    <  2:1408978/133‑1 (MQ=255)
                                                                                                              acCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAacaca                                   >  1:415394/1‑150 (MQ=255)
                                                                                                              acCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAacaca                                   >  2:2027848/1‑150 (MQ=255)
                                                                                                               ccGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAacacac                                  <  1:1233812/150‑1 (MQ=255)
                                                                                                                cGGGAGGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACAcc                                 <  1:650042/150‑1 (MQ=255)
                                                                                                                  ggATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTcgcg                                                               >  1:198640/1‑118 (MQ=255)
                                                                                                                  ggATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTcgcg                                                               <  2:198640/118‑1 (MQ=255)
                                                                                                                  ggATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCt                               >  2:1991125/1‑150 (MQ=255)
                                                                                                                   gATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACTGCCTTGCTTTGAAGCTGTGCTTGGAGAGTAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTc                              >  1:891799/1‑150 (MQ=255)
                                                                                                                   gATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTc                              <  2:2251053/150‑1 (MQ=255)
                                                                                                                    aTGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTcc                             >  1:113940/1‑150 (MQ=255)
                                                                                                                         cTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAacaca                                   <  1:580567/139‑1 (MQ=255)
                                                                                                                         cTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAacaca                                   >  2:580567/1‑139 (MQ=255)
                                                                                                                         cTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTaa                        >  2:132849/1‑150 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATTCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAg                       >  1:1042470/1‑150 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCt                                                                                                        <  2:1078586/70‑1 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCt                                                                                                        >  1:1078586/1‑70 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATa                                                   <  2:2265584/122‑1 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATa                                                   >  1:2265584/1‑122 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAg                       >  2:1674936/1‑150 (MQ=255)
                                                                                                                          tCCGGGCGCGCCTGGCCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAg                       >  2:169890/1‑150 (MQ=255)
                                                                                                                              ggCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGc                   <  1:1431679/150‑1 (MQ=255)
                                                                                                                                     cTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGcactgctt            <  2:889377/150‑1 (MQ=255)
                                                                                                                                      tGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGcactgcttc           >  2:2221382/1‑150 (MQ=255)
                                                                                                                                              ttGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGcc                                                      <  2:1032668/99‑1 (MQ=255)
                                                                                                                                              ttGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGcc                                                      >  1:1032668/1‑99 (MQ=255)
                                                                                                                                              ttGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTg   >  1:628512/1‑150 (MQ=255)
                                                                                                                                              ttGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTg   >  1:1450779/1‑150 (MQ=255)
                                                                                                                                              atGGCACTGTGCCAAGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAagcag                                             >  2:1076485/2‑108 (MQ=255)
                                                                                                                                              atGGCACTGTGCCAAGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAagcag                                             <  1:1076485/107‑1 (MQ=255)
                                                                                                                                              atGGCACTGTGCCAAGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTg   >  1:249545/2‑150 (MQ=255)
                                                                                                                                               tGGCACAGTGCCATGCCACCGCTGAGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGACGccctacgtcactgggttg  <  2:1731006/150‑6 (MQ=255)
                                                                                                                                               tGGCACAGGGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCGGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTgt  <  1:468898/150‑2 (MQ=255)
                                                                                                                                                 gCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCtt                                                                                                       <  2:1512892/48‑1 (MQ=255)
                                                                                                                                                 gCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCtt                                                                                                       >  1:1512892/1‑48 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGGCGCCCTGGGTCAGCGAGGCGACGAAACCTTGTGCCAGTGCTTCGGCACGGGCTTGCAGGGCGGCGTCGAACAGGAACACCGGGTTGATCGAACTCATTTCGGCATACACCGGGATGGGCTCCGGGCGCGCCTGGGCGGCTTGGCACAGTGCCATGCCACCGCTGCGCGAGCCGGTGAAACCAACGGCCTTGATACGTGGGTCGCTGACCAGAGCAATGCCCACTTCGCGGCCGGAGCCATAAAGCAGCGAGAACACACCCTCCGGTAAGCCGCACTGCTTCACTGCTTGC  >  NC_002947/1435076‑1435368

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: