Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,849,916 A→T 5.7% 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 5.7% 1364.4 / 45.9 405L237Q (CTG→CAG) PP_33943‑hydroxy‑3‑methylglutaryl‑CoA lyase
Reads supporting (aligned to +/- strand):  ref base A (203/179);  new base T (11/12);  total (214/191)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 411 total aligned reads displayed.

CGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACC  >  NC_002947/3849767‑3850065
                                                                                                                                                     |                                                                                                                                                     
cGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCa                                                                                                                                                       <  2:4723665/150‑1 (MQ=255)
cGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCa                                                                                                                                                       <  2:3815620/150‑1 (MQ=255)
  cccGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAgg                                                                                                                                                     <  1:2382574/150‑1 (MQ=255)
     gggaGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGccc                                                                                                                                                  <  2:4204173/147‑1 (MQ=255)
     ggTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGccc                                                                                                                                                  <  2:2326598/150‑1 (MQ=255)
     ggTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGAGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGAAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGccc                                                                                                                                                  >  1:2522956/1‑150 (MQ=255)
        aGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGcc                                                                                                                                               <  2:2406211/150‑1 (MQ=255)
            cttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAggg                                                                                                                                           <  1:6967112/150‑1 (MQ=255)
            cttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAggg                                                                                                                                           >  1:5817377/1‑150 (MQ=255)
             ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGa                                                                                                                                          <  1:1025386/150‑1 (MQ=255)
             ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGa                                                                                                                                          >  1:7766886/1‑150 (MQ=255)
             ttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGa                                                                                                                                          <  2:1445336/150‑1 (MQ=255)
                gctgAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGc                                                                                                                                       <  1:7808415/150‑1 (MQ=255)
                 cttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGAGGct                                                                                                                                      <  2:8396150/150‑2 (MQ=255)
                  ttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCAt                                                                                                                                     >  1:2324086/1‑150 (MQ=255)
                     gTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCa                                                                                                                                  <  2:218607/150‑1 (MQ=255)
                     gTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCa                                                                                                                                  <  1:5575583/150‑1 (MQ=255)
                     gTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCa                                                                                                                                  <  1:2882142/150‑1 (MQ=255)
                      tACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAg                                                                                                                                 >  1:3537097/1‑150 (MQ=255)
                       aCCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGt                                                                                                                                >  1:8042305/1‑150 (MQ=255)
                          tGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcc                                                                                                                             >  1:7610570/1‑150 (MQ=255)
                          tGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTcc                                                                                                                             >  1:6544203/1‑150 (MQ=255)
                              gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCaa                                                                                                                         <  2:8282000/150‑1 (MQ=255)
                              gggTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCaa                                                                                                                         >  1:3296540/1‑150 (MQ=255)
                                 tttCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGc                                                                                                                      >  1:7161567/1‑150 (MQ=255)
                                 tttCGATTAGCGCAGGCAAGTCGATGCCAGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGc                                                                                                                      >  1:2729796/1‑150 (MQ=255)
                                     gATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctg                                                                                                                  >  2:7588718/1‑150 (MQ=255)
                                      aTTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgc                                                                                                                 >  1:6862966/1‑150 (MQ=255)
                                       ttAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgct                                                                                                                >  2:4046121/1‑150 (MQ=255)
                                       ttAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctgct                                                                                                                >  1:309578/1‑150 (MQ=255)
                                         aGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCtgctgc                                                                                                              <  2:5963402/150‑1 (MQ=255)
                                          gcgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCa                                                                                                             >  1:2123246/1‑150 (MQ=255)
                                           cgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCAc                                                                                                            >  1:6043300/1‑150 (MQ=255)
                                            gcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCAcc                                                                                                           >  2:3365137/1‑150 (MQ=255)
                                                  aaGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGccacc                                                                                                     >  2:7869044/1‑150 (MQ=255)
                                                    gTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAg                                                                                                   <  1:441310/150‑1 (MQ=255)
                                                       gATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGcg                                                                                                >  2:4488938/1‑149 (MQ=255)
                                                       gATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCt                                                                                                >  2:5633057/1‑150 (MQ=255)
                                                         tGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGt                                                                                              <  2:7101022/150‑1 (MQ=255)
                                                         tGCCGGAGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGt                                                                                              >  1:4488947/1‑150 (MQ=255)
                                                          gCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGtt                                                                                             <  2:5243356/150‑1 (MQ=255)
                                                           ccGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTg                                                                                            <  1:1582018/150‑1 (MQ=255)
                                                            cGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTgg                                                                                           >  2:3132534/1‑150 (MQ=255)
                                                              gtgtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGcc                                                                                         <  1:1842067/150‑1 (MQ=255)
                                                              gtgtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGcc                                                                                         >  1:1143641/1‑150 (MQ=255)
                                                                gtCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCaa                                                                                       <  1:8803877/150‑1 (MQ=255)
                                                                 tCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAag                                                                                      <  1:3611093/150‑1 (MQ=255)
                                                                  cGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAagg                                                                                     <  1:6192887/150‑1 (MQ=255)
                                                                        ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCTtggccaaca                                                                               <  2:197179/150‑4 (MQ=255)
                                                                         cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggcc                                                                              >  2:2499604/1‑150 (MQ=255)
                                                                         cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAaggccaggcc                                                                              >  2:1693732/1‑150 (MQ=255)
                                                                             cTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCgg                                                                          <  2:4330691/150‑1 (MQ=255)
                                                                              tttGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCggg                                                                         >  2:3650967/1‑150 (MQ=255)
                                                                                 gcagcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGgtgt                                                                      <  2:5396571/150‑1 (MQ=255)
                                                                                      aTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                 <  1:7531650/150‑1 (MQ=255)
                                                                                      aTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                 <  1:5692827/150‑1 (MQ=255)
                                                                                      aTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATgg                                                                 <  2:8352351/150‑1 (MQ=255)
                                                                                        gAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGaa                                                               <  1:3049899/150‑1 (MQ=255)
                                                                                          acacCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGt                                                             >  2:4819935/1‑150 (MQ=255)
                                                                                          acacCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGt                                                             <  2:7438826/150‑1 (MQ=255)
                                                                                            acCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGACGTgc                                                           >  1:585395/1‑150 (MQ=255)
                                                                                             ccAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcg                                                          <  2:1207312/150‑1 (MQ=255)
                                                                                              cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGTAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                         >  2:4888707/1‑150 (MQ=255)
                                                                                              cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcgc                                                         >  2:1744387/1‑150 (MQ=255)
                                                                                                  tCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGGTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGcc                                                     <  1:1986920/150‑1 (MQ=255)
                                                                                                    cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                   <  2:3512474/150‑1 (MQ=255)
                                                                                                     tCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTg                                                  <  2:1330549/150‑1 (MQ=255)
                                                                                                        gTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTgaga                                               <  1:4586714/150‑1 (MQ=255)
                                                                                                           aCGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACgg                                            >  2:4498429/1‑150 (MQ=255)
                                                                                                             ggTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTc                                          >  2:8538516/1‑150 (MQ=255)
                                                                                                               tATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGACGTGCGCGGCCTTGAGACGGTCGc                                        >  2:1863768/1‑150 (MQ=255)
                                                                                                               tATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGc                                        >  1:5195598/1‑150 (MQ=255)
                                                                                                                 ttACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCt                                      <  2:6351070/150‑1 (MQ=255)
                                                                                                                   aCCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGc                                    >  1:6376090/1‑150 (MQ=255)
                                                                                                                     cAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCaa                                  >  2:6142727/1‑150 (MQ=255)
                                                                                                                      aGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAt                                 <  2:7266850/150‑1 (MQ=255)
                                                                                                                       gACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtc                                >  1:8192073/1‑150 (MQ=255)
                                                                                                                       gACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtc                                >  2:4949417/1‑150 (MQ=255)
                                                                                                                          gcgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTcc                             >  2:2743997/1‑150 (MQ=255)
                                                                                                                           cgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTccc                            <  2:666488/150‑1 (MQ=255)
                                                                                                                                  ggCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGc                     <  2:5038523/150‑1 (MQ=255)
                                                                                                                                       aCGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGaa                >  2:2859838/1‑150 (MQ=255)
                                                                                                                                        cGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAg               <  2:6992626/150‑1 (MQ=255)
                                                                                                                                        cGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAg               >  2:3971719/1‑150 (MQ=255)
                                                                                                                                             aGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCa          >  1:373244/1‑150 (MQ=255)
                                                                                                                                              gCCTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCACTCTCCCGTACGGCCTGAAGCACCAt         >  2:4627116/1‑150 (MQ=255)
                                                                                                                                               ccTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt        >  1:500213/1‑150 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt        <  2:1839668/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt        <  2:6460564/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt        <  1:4046121/150‑1 (MQ=255)
                                                                                                                                               ccTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt        <  1:5612738/150‑1 (MQ=255)
                                                                                                                                                cTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg       >  2:7971124/1‑150 (MQ=255)
                                                                                                                                                 tCCCTGGCGGGCCTGGGAGGCTGCCCGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc      >  2:2986256/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc      >  1:3165879/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc      >  2:7693173/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc      >  1:1715871/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc      >  1:8399479/1‑150 (MQ=255)
                                                                                                                                                  cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGTCTGAAGCACCATTGcc     >  1:663702/1‑150 (MQ=255)
                                                                                                                                                   ccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCa    >  2:872117/1‑150 (MQ=255)
                                                                                                                                                    cAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCAc   >  1:7438769/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCCGGTAGCGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTGCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACC  >  NC_002947/3849767‑3850065

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: