Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,849,925 A→T 6.4% L234Q (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,9250AT 6.4% 1103.3 / 43.1 329L234Q (CTG→CAG) PP_33943‑hydroxy‑3‑methylglutaryl‑CoA lyase
Reads supporting (aligned to +/- strand):  ref base A (155/152);  new base T (15/6);  total (171/158)
Fisher's exact test for biased strand distribution p-value = 7.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 333 total aligned reads displayed.

GCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGGCCCGGA  >  NC_002947/3849776‑3850074
                                                                                                                                                     |                                                                                                                                                      
gcgcTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCa                                                                                                                                                        <  1:5051253/150‑1 (MQ=255)
 cgcttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAg                                                                                                                                                       <  2:6604970/150‑1 (MQ=255)
    ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGa                                                                                                                                                    >  1:341416/1‑150 (MQ=255)
        cttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCa                                                                                                                                                >  2:1280810/1‑150 (MQ=255)
         ttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCAt                                                                                                                                               >  1:393332/1‑150 (MQ=255)
         ttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCAt                                                                                                                                               >  1:4534981/1‑150 (MQ=255)
         ttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCAt                                                                                                                                               <  2:4487316/150‑1 (MQ=255)
         ttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGATGCCTCCCTGGCGGGCCTGGGAGGctg                                                                                                                                               >  2:2511505/1‑148 (MQ=255)
           aGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATcc                                                                                                                                             <  1:1904202/150‑1 (MQ=255)
            gTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGATGCCTCCCTGGCGGGCCTGGGAGGctgccc                                                                                                                                            <  1:6075896/150‑6 (MQ=255)
             tACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAg                                                                                                                                           <  1:1280810/150‑1 (MQ=255)
             tACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAg                                                                                                                                           >  2:5144932/1‑150 (MQ=255)
                 tGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcc                                                                                                                                      >  2:2741557/1‑150 (MQ=255)
                 tGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcc                                                                                                                                      >  2:6456424/1‑150 (MQ=255)
                  gCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACTGGATGCCTCCCTGGCGGGCCTGGGAGGCTGCCCGTTcgc                                                                                                                                      <  1:5568589/150‑3 (MQ=37)
                   ccGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcccc                                                                                                                                    >  2:220773/1‑150 (MQ=255)
                   ccGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcccc                                                                                                                                    <  1:465261/150‑1 (MQ=255)
                     gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAATGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCaa                                                                                                                                  >  2:5602752/1‑150 (MQ=255)
                     gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCaa                                                                                                                                  >  1:6755559/1‑150 (MQ=255)
                     gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCaa                                                                                                                                  >  1:4780714/1‑150 (MQ=255)
                     gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCACCCAGTTCCCCaa                                                                                                                                  >  1:1725300/1‑150 (MQ=255)
                       gTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATg                                                                                                                                <  1:2261548/150‑1 (MQ=255)
                        tttCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGc                                                                                                                               >  1:3309518/1‑150 (MQ=255)
                        tttCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGc                                                                                                                               >  1:6549938/1‑150 (MQ=255)
                             aTTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgc                                                                                                                          >  2:1558457/1‑150 (MQ=255)
                                  cgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCAc                                                                                                                     >  1:5439453/1‑150 (MQ=255)
                                  cgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCAc                                                                                                                     >  1:1153613/1‑150 (MQ=255)
                                      ggCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGcc                                                                                                                 <  2:585803/150‑1 (MQ=255)
                                         aaGTCGATGCCGGGGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGGACCGccacc                                                                                                              <  2:3561760/150‑1 (MQ=255)
                                            tCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAgg                                                                                                           >  1:232451/1‑150 (MQ=255)
                                                 gCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGtt                                                                                                      <  2:2590576/150‑1 (MQ=255)
                                                  ccGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTg                                                                                                     <  1:3879106/150‑1 (MQ=255)
                                                  ccGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTg                                                                                                     <  2:2389669/150‑1 (MQ=255)
                                                   cGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTgg                                                                                                    >  1:3516768/1‑150 (MQ=255)
                                                     gtgtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGcc                                                                                                  >  2:2911858/1‑150 (MQ=255)
                                                       gtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCaa                                                                                                <  1:5829115/150‑1 (MQ=255)
                                                          gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggc                                                                                             >  2:5301338/1‑150 (MQ=255)
                                                           tAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggcc                                                                                            >  1:6257334/1‑150 (MQ=255)
                                                            aGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggcca                                                                                           <  1:2007238/150‑1 (MQ=255)
                                                            aGCCCATGCTTTGCAACATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggcca                                                                                           <  1:2707069/150‑1 (MQ=255)
                                                               ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                                        >  2:5339329/1‑150 (MQ=255)
                                                               ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                                        >  1:4195368/1‑150 (MQ=255)
                                                                    cTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCgg                                                                                   <  2:4336707/150‑1 (MQ=255)
                                                                    cTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCgg                                                                                   <  2:4833867/150‑1 (MQ=255)
                                                                        gcagcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGAgtgt                                                                               >  1:5224144/1‑150 (MQ=255)
                                                                           gcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCAt                                                                            <  1:2958185/150‑1 (MQ=255)
                                                                            caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCAt                                                                            <  1:2051481/149‑1 (MQ=255)
                                                                             aTGAACACCCGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                          <  2:3483411/150‑1 (MQ=255)
                                                                             aTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                          <  1:6125835/150‑1 (MQ=255)
                                                                               gAACACCAGGTCCTCGGTCACGGTCTTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCCAGGCCAGGCCCCGGGTGTCATGGaa                                                                        >  2:6058920/1‑150 (MQ=255)
                                                                               gAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGaa                                                                        <  2:1682535/150‑1 (MQ=255)
                                                                               gAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGaa                                                                        <  1:2603712/150‑1 (MQ=255)
                                                                               gAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGaa                                                                        <  1:2810886/150‑1 (MQ=255)
                                                                                   acCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgc                                                                    >  1:998705/1‑150 (MQ=255)
                                                                                    ccAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcg                                                                   >  1:4448975/1‑150 (MQ=255)
                                                                                    ccAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcg                                                                   >  1:2564337/1‑150 (MQ=255)
                                                                                     cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                                  <  1:1646126/150‑1 (MQ=255)
                                                                                     cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                                  <  1:1480367/150‑1 (MQ=255)
                                                                                     cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                                  >  1:3575212/1‑150 (MQ=255)
                                                                                     cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                                  >  1:3717514/1‑150 (MQ=255)
                                                                                        gTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGc                                                               <  2:1623565/150‑1 (MQ=255)
                                                                                           cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                            <  1:4072383/150‑1 (MQ=255)
                                                                                           cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                            <  2:3044540/150‑1 (MQ=255)
                                                                                           cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                            <  1:2467523/150‑1 (MQ=255)
                                                                                             cGGTCCCGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGGTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTga                                                          >  2:725879/1‑150 (MQ=255)
                                                                                              ggTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCCAGGCCCGGCCCCGGGTGTCAGGGAAGTGCGCGGCCTTgag                                                         >  1:1397654/1‑150 (MQ=255)
                                                                                               gTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTgaga                                                        >  2:628429/1‑150 (MQ=255)
                                                                                               gTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTgaga                                                        <  1:6477780/150‑1 (MQ=255)
                                                                                                 cACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACg                                                      >  2:6539121/1‑150 (MQ=255)
                                                                                                 cACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACg                                                      >  2:3077360/1‑150 (MQ=255)
                                                                                                  aCGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACgg                                                     <  1:609318/150‑1 (MQ=255)
                                                                                                  aCGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACgg                                                     >  1:3752698/1‑150 (MQ=255)
                                                                                                     gTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCg                                                  >  1:1100711/1‑150 (MQ=255)
                                                                                                        ttACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCt                                               <  1:3038431/150‑1 (MQ=255)
                                                                                                        ttACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCt                                               <  1:4234089/150‑1 (MQ=255)
                                                                                                           ccAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCa                                            <  2:497884/150‑1 (MQ=255)
                                                                                                                cgcgCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtctc                                       <  2:3909657/150‑1 (MQ=255)
                                                                                                                 gcgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTcc                                      <  2:2085694/150‑1 (MQ=255)
                                                                                                                      aGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTAc                                 >  2:4038960/1‑150 (MQ=255)
                                                                                                                        gggCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCACTCTCCCGTACgg                               >  2:4232783/1‑150 (MQ=255)
                                                                                                                          gCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCTATGCCCTGCTGGTGCGCCACCAGGCTGTTGGGCAAGGCCAGGCCCCGGGTGTCATGGAATTGCGCGGCCTTGAGACGGTCGCCTGCTATCTCCCGTACGGcc                             <  2:5011310/150‑1 (MQ=255)
                                                                                                                           cGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCt                            <  2:1899564/150‑1 (MQ=255)
                                                                                                                            gAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGGCGCCTGCAATCTCCCGTACGGCATg                           >  1:315653/1‑150 (MQ=255)
                                                                                                                             aaCGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGa                          >  2:2815021/1‑150 (MQ=255)
                                                                                                                              aCGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGaa                         <  1:1340992/150‑1 (MQ=255)
                                                                                                                              aCGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGaa                         <  1:949518/150‑1 (MQ=255)
                                                                                                                                   atgCCTCCCTGGCGGGCCTGGGAGGCTGCCCGTTCGCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCAc                     >  2:1479105/3‑150 (MQ=38)
                                                                                                                                     gCCTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAt                  >  1:755123/1‑150 (MQ=255)
                                                                                                                                       cTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg                >  1:1057933/1‑150 (MQ=255)
                                                                                                                                       cTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg                <  1:498463/150‑1 (MQ=255)
                                                                                                                                        tCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc               >  2:5217964/1‑150 (MQ=255)
                                                                                                                                        tCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc               >  2:5927601/1‑150 (MQ=255)
                                                                                                                                        tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc               >  1:3811901/1‑150 (MQ=255)
                                                                                                                                        tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc               >  1:5654396/1‑150 (MQ=255)
                                                                                                                                         cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGcc              >  2:3692183/1‑150 (MQ=255)
                                                                                                                                         cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGcc              >  2:272698/1‑150 (MQ=255)
                                                                                                                                              gCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTg         >  1:1798639/1‑150 (MQ=255)
                                                                                                                                               cccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTgg        <  1:3057249/150‑1 (MQ=255)
                                                                                                                                                ccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGGc       >  1:5974442/1‑150 (MQ=255)
                                                                                                                                                ccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGGc       >  2:1350989/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGGCCCGGA  >  NC_002947/3849776‑3850074

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: