Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,849,916 A→T 6.8% L237Q (CTG→CAG)  PP_3394 ← 3‑hydroxy‑3‑methylglutaryl‑CoA lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,849,9160AT 6.8% 1228.8 / 51.1 368L237Q (CTG→CAG) PP_33943‑hydroxy‑3‑methylglutaryl‑CoA lyase
Reads supporting (aligned to +/- strand):  ref base A (183/159);  new base T (15/10);  total (198/170)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 369 total aligned reads displayed.

CGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCAC  >  NC_002947/3849767‑3850064
                                                                                                                                                     |                                                                                                                                                    
  cccGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAgg                                                                                                                                                    >  1:4213863/1‑150 (MQ=255)
    cGGTAGCGCTTCGCTTAGCACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGcc                                                                                                                                                  >  2:5768737/1‑150 (MQ=255)
     ggTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGccc                                                                                                                                                 >  1:2523703/1‑150 (MQ=255)
        aGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGcc                                                                                                                                              <  2:2189631/150‑1 (MQ=255)
           gcttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAgg                                                                                                                                           <  1:1836262/150‑1 (MQ=255)
             ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGa                                                                                                                                         <  1:5236250/150‑1 (MQ=255)
             ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGa                                                                                                                                         >  1:4743929/1‑150 (MQ=255)
              tcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAg                                                                                                                                        <  2:1334810/150‑1 (MQ=255)
              tcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAg                                                                                                                                        <  1:674795/150‑1 (MQ=255)
                   tAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGAGGctgc                                                                                                                                   <  2:3120139/150‑4 (MQ=255)
                   tAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATc                                                                                                                                   <  1:3782147/150‑1 (MQ=255)
                   gaGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGAGGctgc                                                                                                                                   <  1:942427/149‑4 (MQ=255)
                    aGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATcc                                                                                                                                  >  2:6514158/1‑150 (MQ=255)
                        ccTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGtt                                                                                                                              >  2:5873985/1‑150 (MQ=255)
                         cTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTc                                                                                                                             <  2:1573514/150‑1 (MQ=255)
                          tGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcc                                                                                                                            <  2:2052311/150‑1 (MQ=255)
                           gCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCGGTTccc                                                                                                                           <  2:4923889/150‑1 (MQ=255)
                           gCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTccc                                                                                                                           >  2:3883975/1‑150 (MQ=255)
                                     gATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctg                                                                                                                 >  2:1037331/1‑150 (MQ=255)
                                     gATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctg                                                                                                                 >  1:134851/1‑150 (MQ=255)
                                      aTTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgc                                                                                                                >  1:2602627/1‑150 (MQ=255)
                                      aTTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgc                                                                                                                >  1:4405368/1‑150 (MQ=255)
                                        cccgccAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCGTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgctg                                                                                                              >  2:1805062/6‑150 (MQ=255)
                                         aGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCtgctgc                                                                                                             <  2:2305031/150‑1 (MQ=255)
                                          gcgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCa                                                                                                            >  2:5399145/1‑150 (MQ=255)
                                            gcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCAcc                                                                                                          <  1:805250/150‑1 (MQ=255)
                                                gCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGcca                                                                                                      >  1:3564176/1‑150 (MQ=255)
                                                   aGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGccacca                                                                                                   >  1:4065643/1‑150 (MQ=255)
                                                     tCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAgg                                                                                                 >  1:6301419/1‑150 (MQ=255)
                                                     tCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAgg                                                                                                 >  1:3105185/1‑150 (MQ=255)
                                                     tCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAgg                                                                                                 <  1:1240917/150‑1 (MQ=255)
                                                       gATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCt                                                                                               <  1:6100618/150‑1 (MQ=255)
                                                        aTGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTg                                                                                              >  1:4844454/1‑150 (MQ=255)
                                                          gCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGtt                                                                                            <  1:2840894/150‑1 (MQ=255)
                                                           ccGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTg                                                                                           >  2:5016770/1‑150 (MQ=255)
                                                            cGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTgg                                                                                          >  1:2983746/1‑150 (MQ=255)
                                                              gtgtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGcc                                                                                        >  2:3635877/1‑150 (MQ=255)
                                                               tgtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCa                                                                                       >  2:5668771/1‑150 (MQ=255)
                                                                gtCTTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCaa                                                                                      <  2:7365668/150‑1 (MQ=255)
                                                                gtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGcggt                                                                                      <  1:6364920/150‑4 (MQ=255)
                                                                 tCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAag                                                                                     >  1:598767/1‑150 (MQ=255)
                                                                   gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggc                                                                                   >  1:2142033/1‑150 (MQ=255)
                                                                   gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggc                                                                                   >  2:7110267/1‑150 (MQ=255)
                                                                     aGCCCATGCTTTGCAGCATGACCACCAGGTCCTCGGTCACGGTATGACCAGCCGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggcca                                                                                 <  2:6092981/150‑1 (MQ=255)
                                                                     aGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggcca                                                                                 <  2:588612/150‑1 (MQ=255)
                                                                     aGCCCATGCTTTGCAGATTGTACACCATGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGTCGAACGGGCAGGCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCTATGCCCTGCTGCACCGCCACCAGGCTGTTGGCGAagtcca                                                                                 <  2:1104996/150‑1 (MQ=255)
                                                                      gCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccag                                                                                >  2:2508722/1‑150 (MQ=255)
                                                                      gCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccag                                                                                >  2:7209726/1‑150 (MQ=255)
                                                                      gCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTAACAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGACCGCCAGGGAGGAATCCAGTTCCACAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccag                                                                                >  1:6528663/1‑150 (MQ=255)
                                                                        ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                              >  1:1099485/1‑150 (MQ=255)
                                                                        ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                              <  1:7378573/150‑1 (MQ=255)
                                                                        ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                              >  1:7299322/1‑150 (MQ=255)
                                                                         cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggcc                                                                             >  2:1215137/1‑150 (MQ=255)
                                                                         cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggcc                                                                             >  2:2918575/1‑150 (MQ=255)
                                                                           tGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGcccc                                                                           >  2:5611375/1‑150 (MQ=255)
                                                                                 gcagcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGgtgt                                                                     >  1:556563/1‑150 (MQ=255)
                                                                                   agcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCa                                                                   >  2:372757/1‑150 (MQ=255)
                                                                                     caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATg                                                                 <  1:4537457/150‑1 (MQ=255)
                                                                                     caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGATGCCTCCCTGGCGGGCCTGGGAGGCTGCCCGTTcgc                                                                                                                           >  2:5295056/1‑90 (MQ=18)
                                                                                     caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGATGCCTCCCTGGCGGGCCTGGGAGGCTGCCCGTTcgc                                                                                                                           <  1:5295056/92‑3 (MQ=18)
                                                                                      aTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                <  1:5387248/150‑1 (MQ=255)
                                                                                        gAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGaa                                                              <  2:628922/150‑1 (MQ=255)
                                                                                          acacCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGt                                                            <  2:3335795/150‑1 (MQ=255)
                                                                                              cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                        >  2:2790844/1‑150 (MQ=255)
                                                                                              cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                        <  2:3266307/150‑1 (MQ=255)
                                                                                                 gTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGc                                                     >  2:750574/1‑150 (MQ=255)
                                                                                                  tCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGcc                                                    >  1:7123900/1‑150 (MQ=255)
                                                                                                    cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                  >  1:6557351/1‑150 (MQ=255)
                                                                                                     tCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTg                                                 >  2:6916404/1‑150 (MQ=255)
                                                                                                     tCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTg                                                 >  2:1264651/1‑150 (MQ=255)
                                                                                                         tCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGAc                                             >  1:393841/1‑150 (MQ=255)
                                                                                                           aCGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACgg                                           >  1:7326964/1‑150 (MQ=255)
                                                                                                              gTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCg                                        >  2:4294913/1‑150 (MQ=255)
                                                                                                               tATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGc                                       <  1:4326234/150‑1 (MQ=255)
                                                                                                                 ttACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCt                                     <  1:4150134/150‑1 (MQ=255)
                                                                                                                  tACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTg                                    >  2:2639912/1‑150 (MQ=255)
                                                                                                                    ccAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCt                                     >  1:5772314/1‑147 (MQ=255)
                                                                                                                    ccAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCa                                  >  1:1573226/1‑150 (MQ=255)
                                                                                                                    ccAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCa                                  <  1:3645356/150‑1 (MQ=255)
                                                                                                                       gACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtc                               <  1:4910658/150‑1 (MQ=255)
                                                                                                                           cgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTccc                           <  2:4075772/150‑1 (MQ=255)
                                                                                                                            gcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGGTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCg                          <  2:3012120/150‑1 (MQ=255)
                                                                                                                                        cGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAg              <  2:4669058/150‑1 (MQ=255)
                                                                                                                                         gggCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGc             <  1:6286315/150‑1 (MQ=255)
                                                                                                                                           gatgCCTCCCTGGCGGGCCTGGGAGGCTTCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCAc           >  2:3480199/4‑150 (MQ=255)
                                                                                                                                           gCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCAc           >  1:1035194/1‑150 (MQ=255)
                                                                                                                                              gCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGAGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAt        >  2:6121043/1‑150 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt       <  1:4583244/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt       >  2:4751911/1‑150 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt       <  2:4213863/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt       <  2:5981865/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGCACGGCCTGAAGCACCAtt       <  1:7000406/150‑1 (MQ=255)
                                                                                                                                                cTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg      >  1:3367974/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc     >  1:2027703/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc     >  2:4362875/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc     >  2:4540640/1‑150 (MQ=255)
                                                                                                                                                  cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGcc    >  1:2271886/1‑150 (MQ=255)
                                                                                                                                                  cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGcc    >  1:3416414/1‑150 (MQ=255)
                                                                                                                                                   ccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCa   <  1:5873985/150‑1 (MQ=255)
                                                                                                                                                    cAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCAc  >  2:6762052/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCAC  >  NC_002947/3849767‑3850064

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